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ABSTRACT 


This  report  presents  an  approximate  theory  of  the  growth  of  a  point 
source  explosion  in  an  infinite  solid  or  liquid  medium  in  the  range  of  shock 
pressure  below  that  required  to  induce  a  physical  change  of  phase  of  the  ma¬ 
terial,  and  above  pressures  at  which  the  wave  propagation  can  be  considered 
acoustic.  The  analysis  is  based  on  the  characterization  of  the  shock  state 
locus  for  the  material  by  a  general  form  containing  three  constants  related 
to  the  compressibility  of  the  material.  Both  porous  and  non-porous  solids 
are  considered. 

Peak  shock  pressure  and  arrival  time  data  obtained  by  machine 
calculation  as  functions  of  shock  radius  are  presented  graphically  and 
numerically  in  dimensionless  form  for  several  values  of  the  compressibility 
constants  involved.  Methods  of  deducing  values  of  the  constants  from 
compressibility  data  for  particular  materials  are  presented.  The  results 
are  applicable  to  estimates  of  shock  attenuation  in  a  variety  of  blast  absorbs 
ing  or  shielding  materials^ 
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POINT  SOURCE  EXPLOSION  IN  A  SQUID 


L  INTRODUCTION 

Recently  interest  has  arisen  in  methods  of  reducing  the  mechanical 
or  blast  effects  of  a  sudden,  explosion-like  release  of  energy.  This  interest 
has  grown  as  a  result  of  the  need  for  defensive  or  protective  measures 
against  the  high  energy  yield  explosive  weapons  of  modern  warfare,  and 
the  need  for  the  confinement  and  direction  of  large  quantities  of  energy  in 
peacetime  applications  of  powerful  energy  sources.  The  concept  of  confine¬ 
ment  or  containment  arises  in  explosion  phenomenology  itself,  apart  from 
any  protective  considerations,  when  a  sudden  release  of  energy  takes  place 
at  a  point  in  a  confining  medium  such  as  soil  or  water.  The  propagation  of 
a  strong  disturbance  in  a  condensed  medium  is  strongly  controlled  by  the 
ability  of  the  material  to  abstract  mechanical  energy  from  the  disturbance 
as  it  progresses  through  the  medium,  and  to  degrade  that  energy  to  a  non¬ 
destructive  form  through  dissipative  processes.  These  characteristics  of 
certain  solid  materials  and  liquids  can  result  in  the  rapid  attenuation  of  the 
shock  pressures  due  to  3.11  explosion,  in  the  mediumo 

There  is  limited  information  in  existence  on  the  type  and  magnitude 
of  disturbance  which  can  result  from  the  serious  malfunctioning  of  a  nuclear 
reactor^.  It  is  partly  because  of  this  element  of  uncertainty  that  stringent 
safety  and  precautionary  measures  are  generally  adopted,  particularly  in 
the  cases  of  high-output  thermal  reactor  power  generators.  Such  power 
facilities  will,  at  least  eventually,  be  necessarily  located  in  heavily  popu¬ 
lated  areas,  and  this  condition  denies  isolation  as  the  sole  safety  provision. 
^Superscript  numerals  refer  to  bibliography  section. 
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Thus  the  concept  of  the  complete  containment  of  an  incident  within  the 
confines  of  the  facility  arises  as  a  prime  safety  provision. 

The  concept  of  explosion  confinement  has  been  applied  to  nuclear 
reactor  safety  through  the  placement  of  blast  absorbing  components, 
referred  to  as  blast  shields,  surrounding  the  reactor  core,  which  is 
considered  as  the  potential  origin  of  a  catastrophic  disturbance.  Such 
shields  possess  the  ability  to  reduce  explosion  pressures  by  several  orders 
of  magnitude  and  thus  protect  an  enclosing  gas-tight  containment  envelope. 
One  suggested  design  for  a  blast  shielding  device  is  a  Seuidwich  type  of 
construction  of  alternating  rigid  and  compressible  layers.  The  essential 
features  of  this  design  have  been  discussed  elsewhere  ;  it  suffices  to  note 
that  its  shielding  characteristics  are  regarded  as  derived  principally  from 
the  blast  attenuating  properties  of  the  compressible  solid  material. 

In  order  to  place  the  concepts  of  blast  absorption  and  blast  shielding 
devices  on  a  firm  basis,  a  knowledge  of  the  growth  of  explosions  in  condensed 
media  is  necessary.  The  purpose  of  this  report  is  to  present  the  results 
of  a  theory  of  explosion  growth  eind  shock  propagation  in  a  solid  or  liquid 
medium,  in  a  form  useful  for  application  to  reactor  safety  analyses  and 
blast  shield  designs. 

The  present  analysis  attempts  to  predict  at  least  the  gross  features 
of  the  mechanical  output  from  a  point- source  explosion,  i.  e.  ,  the  behavior 
of  the  flow  variables  within  and  on  the  explosion-engulfed  sphere  as  the 
disturbance  advances  through  the  medium.  Certain  simplifying  assumptions 
restrict  the  analysis  to  the  region  of  strong  shocks,  where  shock  pressures 
are  high  enough  to  consider  the  material  as  fluid-like  in  behavior.  The 
analysis  rests  to  a  certain  extent  on  over-all  energy  considerations;  no 
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a.ttGmpt  is  made  to  predict  in  det3.il  the  history  of  conditions  n.e3r  the  point 
of  explosiono  Fortunately,  the  more  severe  assumptions  are  applied  in 
regard  to  behavior  near  the  explosion  source  and  involve  only  a  small 
fraction  of  the  explosion^engulfed  spherical  volume  during  most  of  the  time 
history  of  the  phenomena  of  interesto  Thus  the  resulting  peak  shock  pressure 
decays  and  arrival  times  are  not  strongly  affected  by  these  simplifications. 

The  techniques  utilized  in  the  present  report  have  been  developed 

3  4 

in  connection  with  research  on  the  effects  of  underwater  and  underground 
explosions.  Predictions  which  have  been  made  regarding  shock  pressure 
decays  and  times  of  arrival  using  these  techniques  have  met  with  substantial 
experimental  verification  in  large  =  scale  tests  of  nuclear  weapons.  The 
results  presented  here  are  only  those  directly  of  interest  to  the  evaluation 
of  blast  absorbing  materials. 

One  effect  of  the  shock  process  on  a  particle  of  material  is  to  make 
a  certain  fraction  of  the  energy  imparted  to  the  engulfed  particle  imre  cover - 
able  as  mechanical  work.  The  analytical  method  involved  here  consists  of 
evaluating  this  *'waste  heat"  and  abstracting  it  from  the  total  available 
energy  as  the  spherical  shock  front  advances  through  the  material.  This 
permits  the  writing  of  a  differential  relation  connecting  the  material  compres 
sion  to  distance  from  the  explosion  source  and  to  the  explosion  yield.  The 
evaluation  of  the  waste  heat  requires  assumptions  on  the  relationship  of 
shock  pressure  to  compression  of  the  particular  material  considered  and 
on  the  subsequent  expansion  experienced  by  the  material.  This  set  of 
relationships  then  admits  solutions  in  non-dimensional  form  for  shock 
pressure  and  shock  arrival  time  as  functions  of  distance  from  the  explosion 
source  for  arbitrary  values  of  parameters  related  to  the  compressibility 
of  the  condensed  medium. 
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The  solutions  apply  to  regions  in  which  the  material,  although 
crushed,  does  not  undergo  a  physical  chamge  of  phase.  An  assumption  on 
the  effects  of  the  early  history  of  the  explosion  permits  the  statement  of 
initial  conditions  on  the  problem  in  the  region  where  the  formulation  is 
applicable.  No  attempt  is  made  to  study  in  detail  conditions  at  very  early 
times,  since  the  materials  considered  are  of  interest  as  blast  absorbers 
reasonably  distant  from  the  explosion  source.  Methods  are  developed  for 
determining  the  values  of  the  constants  describing  the  compressibility  of 
the  materials  which  are  involved  as  parameters  in  the  solutions  from 
compressibility  data.  Thus,  the  rapid  comparison  of  the  effectiveness  of 
a  wide  variety  of  specific  materials  as  shock  pressure  reducers  is  possible. 

Lastly,  this  report  discusses  the  application  of  the  results  of  the 
calculations  to  real  materials.  Although  these  results  are  for  highly 
idealized  physical  situations,  they  provide  information  on  the  relative 
effectiveness  of  the  materials  as  shock  absorbers. 

II.  DISCUSSION  OF  THE  PROBLEM 

The  present  study  concerns  itself  with  the  growth  of  a  point  source 
explosion  in  a  solid  or  liquid  medium  in  the  range  of  shock  pressure  below 
that  required  to  induce  a  change  of  phase  of  the  material,  and  above  pressures 
at  which  the  propagation  can  be  considered  acoustic.  As  an  initial  state  for 
the  solution  in  this  range,  conditions  are  deduced  from  a  grossly  simplified 
representation  of  the  behavior  of  the  explosion  durings  its  early  history. 

Of  particular  interest  here  are  the  readily  crushable  solids  in  which 
efficient  acoustic  shock  transmission  does  not  occur  until  the  shock  pressure 
decays  to  quite  low  values.  Such  materials  are  worth  consideration  as 
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blast  absorbers,  particularly  those  which  are  resistant  to  ignition  at  very 
high  shock  pressures.  These  requirements  suggest  that  porous  solids,  such 
as  foamed  ceramic  or  siliceous  materials,  may  be  quite  effective  as  blast 
shielding  materials.  Some  of  the  assumptions  of  compressible  fluid  dynamics 
are  therefore  applicable  in  the  range  of  interest  of  the  present  analysis, 
where  the  explosion-generated  shock  pressures  are  high  enough  relative 
to  the  strength  of  the  medium  to  permit  consideration  of  the  material  as 
behaving  like  a  fluid.  The  limitations  of  the  fluid  dynamic  model  are  dis¬ 
cussed  in  Appendix  A.  The  objective  of  design  for  explosion  confinement 
itself  would  be  to  reduce  shock  pressures  of  this  magnitude  in  solids  to 
seismic  levels,  for  which  a  large  body  of  analysis  is  available. 

In  the  present  study,  we  consider  the  problem  of  the  sudden  intro¬ 
duction  of  energy  at  a  point  in  an  infinite  medium,  neglecting  the  e-fect  of 
the  mass  of  any  reaction  products.  The  spherical  symmetry  of  the  resulting 
wave  propagation  permits  the  description  of  phenomena  in  terms  of  radial 
distance  from  the  explosion  source,  and  time  measured  from  the  instant 

of  explosion^ 

III,  PROCEDURES  AND  ASSUMPTIONS 

The  usual  approach  to  the  problem  of  a  point-source  release  of 
energy  into  a  fluid  medium  employs  the  differential  equations  of  motion, 
continuity,  and  some  form  of  the  statement  of  conservation  of  energy  for  the 
region  engulfed  by  the  advancing  spherical  shock.  These  equations  can,  in 
principle,  be  simultaneously  integrated  at  least  numerically  by  machine 
calculation,  provided  an  appropriate  equation  of  state  for  the  material  is 
given.  The  integration  starts  from  a  specified  initial  state,  and  is  subject 
to  the  conditions  at  the  moving  boundary  defined  by  the  shock  front.  The 
shock  front  conditions  are  in  turn  established  by  the  laws  of  conservation 
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of  mass,  momentum,  and  energy  (Hugoniot  conditions)  across  the  shock  front. 

5  6i7  8 

There  are  a  number  of  solutions  to  the  point  blast  problem  in  air  and  water  * 
all  presiiming  equation~of~state  data  of  some  type  for  the  medium, 

A,  Flow  Variable  Wave  Forms 

In  the  absence  of  adequate  data  for  the  formulation  of  a  sufficiently 
general  equation  of  state  for  the  class  of  materials  considered,  we  proceed 
directly  to  the  application  of  the  equations  of  motion  and  continuity  for  the 
engulfed  region  through  an  assumption  for  the  distribution  of  material  density 
at  an  arbitrary  instant  during  the  motion.  The  assvimed  density  profile 
contains  a  time -dependent  parameter,  the  value  of  which  is  determined  by 
the  shock  front  conditions  through  consideration  of  the  conservation  of  mass 
over  the  entire  shock -engulfed  region. 

The  particle  velocity  and  pressure  profiles  consistent  with  the  assumed 
density  distribution  are  obtained  by  direct  integration  of  the  continuity  and 
motion  equations.  For  the  strong  shock  case  in  which  the  peak  density 
varies  relatively  slowly  with  time,  this  process  leads  to  relatively  simple 
variations  of  particle  velocity  and  pressure  with  distance  from  the  explosion 
source, 

B.  Shock  Compres  sion  Locus 

The  wave  forms  of  density,  particle  velocity,  and  pressure  contain 
the  shock  front  values  of  the  respective  quantities.  The  relationship  of 
particle  velocity  to  pressure  and.  density  at  the  shock  is  established  by 
Hugoniot  conditions.  We  characterize  the  relationship  of  shock  pressure 
to  shock  compression  or  density  ratio  (the  Hugoniot  relation)  by  a  general 
form  containing  three  constants  related  to  the  compressibility  of  the 
material.  The  form  is  sufficiently  flexible  to  describe  the  properties  of 
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a  wide  variety  of  materials  vinder  shock.  Methods  are  subsequently  given 
for  deducing  the  values  of  the  constants  for  specific  materials  from  compres¬ 
sibility  data, 

C.  Energy  Partition 

An  equation  is  written  describing  the  division  of  the  total  energy 
yield  into  two  parts  at  any  fixed  instant  during  the  motion.  One  part  consists 
of  kinetic  and  hydrodynamic  energies  available  to  perform  destructive  mechani' 
cal  work  on  the  materied  during  later  motion,  and  distributed  over  the  entire 
explosion -engulfed  sphere  at  the  instant  of  time  under  consideration.  The 
second  part  consists  of  energy  made  unavailable  to  do  mechanical  work  at 
a  later  time  by  non-equilibrium  processes,  referred  to  as  waste  heat,  and 
presumably  evidenced  as  a  temperature  rise  in  the  material.  Now  the  shock 
process  is  considered  as  setting  up  conditions  under  which  energy  is  degraded 
to  waste  heat;  it  remains  to  evaluate  the  quantity  of  energy  committed  to  the 
unavailable  fraction  as  the  shock  front  advances  through  the  material.  This 
means  that  the  waste  heat  term  in  the  energy  balance  equation  can  be  evalu¬ 
ated  over  the  previous  history  of  the  shock  front. 

D.  Expansion  Process  and  Energy  Dissipation 

For  the  condensed  media  of  interest  in  the  present  analysis,  it  is 
assumed  for  purposes  of  the  waste  heat  calculation  that  the  expansion  under¬ 
gone  by  shocked  particles  of  material  is  equivalent  to  expansion  to  the 
pre- shock  specific  volume  along  the  Hugoniot.  The  validity  of  this  approxi¬ 
mation  is  discussed  in  Appendix  B.  It  is  shown  that,  for  the  condensed 
materials  of  present  interest,  all  chauiges  of  state  are  represented  by 
nearly  coincident  traces  in  the  pressure-specific  volume  when  no  change 
of  phase  occurs.  The  assumption  on  the  expansion  process  permits  (1) 


-7- 


Report  No, 
ARF  4132-6 


calculation  of  the  energy  recoverable  as  mechanical  "works  and  (2)  calcu¬ 
lation  of  the  waste  heat  as  the  difference  between  the  internal  energy  im¬ 
parted  to  the  material  during  shock  and  the  recoverable  mechanical  energy. 

E,  Calculation  of  Pressures  and  Arrival  Times 

The  energy  balance  equation  is  reduced  to  a  differential  equation 
involving  only  shock  compression  or  density  ratio,  distance  from  the 
explosion  source,  and  explosion  yield  which,  in  "view  of  the  Hugoniot  relation, 
defines  the  correspondence  between  shock  pressure  and  distance.  The 
constants  describing  the  compressibility  of  the  material  under  shock  appear 
as  parameters  in  the  equation.  Integration  was  performed  numerically  on 
an  IBM  650  digital  computer.  The  results  are  presented  graphically  and 
tabularly  in  the  form  of  a  two -parameter  family  of  pressure-distance  curves, 
the  parameters  being  two  of  those  appearing  in  the  assumed  Hugoniot  relation. 
Integration  was  performed  from  an  arbitrarily  selected  starting  condition, 
but  it  is  shown  that  a  change  in  the  starting  condition  corresponds  to  a  simple, 
one-directional  translation  of  the  pres  sure -distance  curves  on  a  logarithmic 
chart. 

The  non-dimensional  time  of  arrival  of  the  shock  front  as  a  function 
of  distance  from  the  explosion  source  is  easily  obtained  as  the  resvilt  of 
integration  of  the  reciprocal  of  the  shock  front  velocity  over  the  distance 
traversed  by  the  shock  front.  This  integration  was  also  performed  numeri¬ 
cally  on  an  IBM  650  digital  computer;  the  results  are  again  a  two-parameter 
family  of  curves  and  are  obtained  for  the  arbitrarily  selected  starting 
condition,  it  is  shown  that  a  change  in  the  starting  condition  corresponds  to 
equal  translation  in  two  directions  of  the  arrival  time-distance  curves  on  a 
logarithmic  chart. 
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IV. 


CONDITIONS  IN  THE  ENGULFED  REGION 

In  this  section  conditions  at  the  shock  front  are  summarized.  The 
velocity  distribution  and  pressure  profile  consistent  with  an  assumed  density 
distribution  in  space  at  a  fixed  instant  of  time  after  the  instant  of  explosion 
are  deduced. 

A.  Shock  Front  Conditions 

When  a  discontinuity  in  state  variables  exists  in  a  continuous  medium, 
the  disturbance  propagates  into  the  region  of  lower  pressure  and  density. 

In  particular,  if  the  discontinuity  or  shock  moves  into  material  at  rest, 
application  of  conservation  of  mass  across  the  shock  front  leads  to  a  relation 
among  the  shock  front  velocity  U,  the  particle  velocity  u^  following  the 

shock  front,  and  the  densities  and  before  and  immediately  following 

the  shock: 

Po  U  .  (U  -  U3). 

The  subscripts  o  and  s  refer  to  ambient  and  shock  front  conditions, 

respectively. 

Application  of  conservation  of  momentum  across  the  front  leads  to 

a  relation  between  the  absolute  pressures  and  before  and 

immediately  following  the  shock,  and  the  shock  and  particle  velocities: 

p  -  p  =  p  U  u^  . 
s  o  s 

Combination  of  Eqs.  (4.  1)  and  (4.  2)  leads  alternatively  to 


SO 

=  V  -  \r 

o  s 


(4.  3) 


or  to 


(4.  4) 
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where  the  specific  volume  Y  =  l/p.  The  subscript  s  denotes  conditions 
immediately  following  shock. 

Equations  (4.  1)  through  (4,  4)  hold  locally  whether  or  not  the  shock 
front  is  a  plane  surface  and  irrespective  of  the  properties  of  the  medium, 
presuming  only  the  existence  of  shock  conditions. 

The  total  energy  input  E™  to  a  unit  mass  of  material  during  shock 
is  easily  obtained  as 

(^-5) 

2 

and  since  the  unit  mass  acquires  kinetic  energy  ~  ^ internal 

energy  increment  per  unit  mass  is  obtained  through  Eqs,  (4,  4)  and  (4,  5) 
as  the  difference  between  this  total  energy  and  the  kinetic  energy  of  net 
translation  in  the  form 

=1.  =  I  -  '"s’ 

In  further  work,  we  neglect  the  ambient  pressure  in  comparison  to  the 

strong  shock  pressure  P  .  hispection  of  Eqs,  (4,4),  (4,5),  and  (4,  6) 

s 

shows  that  the  total  energy  is  partitioned  equally  between  internal  and 
kinetic  energy  in  the  strong  shock  case. 

Now  if  one  could  specify  an  equation  of  state  connecting  internal  energy 
with  pressure  and  specific  volume  or  with  other  variables  related  through 
the  laws  of  thermodynamics  to  pressure  and  specific  volume,  Eq,  (4,  6) 
would  lead  to  a  Hugoniot  relation  defining  a  locus  of  states  in  the  P-V 
plane  which  can  be  reached  by  shock  from  the  given  ambient  conditions. 
Instead,  however,  we  assume  a  Hugoniot  in  the  following  form: 
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n 


(4.  7) 

In  this  expression  the  constants  C,  and  n  are  regarded  as  characteristic 
of  the  medium.  The  form  exhibits  the  properties  of  (1)  decreasing  compres¬ 
sibility  with  increasing  density  characteristic  of  a  wide  variety  of  materials, 
and  (2)  a  limiting  shock  compression  or  density  ratio  (J-  at  infinite  shock 
pressures.  The  incorporation  of  three  parameters  related  to  the  material 
allows  a  good  deal  of  freedom  in  empirically  describing  the  behavior  of 
actual  materials.  In  the  discussion  following,  the  ambient  pressure 
in  Eq.  (4.  7)  is  neglected. 

Interesting  geometric  interpretations  of  the  energy  partition  and 
the  shock  and  material  velocities  are  afforded  by  the  plot  of  Eq.  (4.  7)  in 
Fig.  1; 

1.  The  total  shock  energy  input  to  a  unit  mass  of  material  is 

numerically  equal  to  the  area  of  the  rectangle  of  altitude  and  base 

V  -  . 

o  s 

2.  The  chord  of  the  Hug oniot  joining  the  points  (P^,  V^)  and 
(0,  V^)  divides  the  recteingle  into  two  equal  triangles,  each  numerically 
equal  in  area  to  the  kinetic  and  internal  energies. 

3.  The  square  of  the  particle  velocity  is  numerically  equal  to 

the  area  of  the  entire  rectangle  in  view  of  Eq.  (4.  4)  and  the  neglection  of 

P  . 
o 

4.  Finally,  £q.  (4.  3)  indicates  that  the  square  of  the  shock 
velocity  is  proportional  to  the  negative  slope  of  the  chord  of  the  Hugoniot.,, 
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If  is  not  neglected  in  comparison  with  ,  observations 
(1)  and  (4)  above  still  hold  when  the  Hugoniot  is  constructed  through  the 
point  V’^)*  In  particular,  for  small  shock  overpressures  (P^  -  P^) 
the  chord  of  the  Hugoniot  approaches  the  tangent  to  the  curve  at  (P^,  V^), 

1/  2 

and  the  shock  velocity  approaches  the  acoustic  velocity  c^  =  -  (dP/dV)y_y 

o 

because  the  Hugoniot  coincides  locally  with  the  adiabat  through  (P^* 

B.  Density  Profile 

Following  the  discussion  of  Section  III-A,  we  assume  that  the  density 
p  as  a  function  of  distance  r  at  any  time  t  after  the  instant  of  explo¬ 


sion  is  given  by: 


P  =  P, 


r  ^ 

(4)  . 


(4.8) 


where  R  is  the  radius  of  the  spherical  shock  front  and  q  is  permitted  to  be 
a  function  of  time.  Neglecting  the  mass  of  any  reaction  products  in  compari¬ 


son  to  the  mass  of  the  engulfed  medium,  we  equate  the  total  mass  contained 
in  the  sphere  of  radius  -R  before  the  explosion  to  that  at  the  instant  of  time 
under  consideration  as  follows: 


R 

(4.9) 
o 

Inserting  the  assumed  density  distribution  into  Eq,  (4.  9)  and  performing  the 
indicated  operation  leads  tor 


4^ 

-3  P< 


2  , 

r  dr. 


(4.  10) 


and  q  is  a  function  of  time  in  view  of  the  as  yet  unknown  variation  of 
with  time. 
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c. 


Particle  Velocity  Profi!^ 


Now  consider  the  mass  (dif/ 3)p  r  enclosed  initially  by  the  sphere 

o  o 


of  radius  r^.  After  passage  of  the  shock  front  and  as  a  result  of  motion  of 
the  fluid  particles,  the  same  mass  is  enclosed  in  a  sphere  of  radius  r, 
where  r(t)  is  the  trajectory  of  the  particle  initially  at  r^  (Fig«  2),  Thus 


4tt  3  .  f  2 

^r  p  =  4  rr  p  r  dr. 
3  o  *^o  I 


(4.  11) 


o 


Insertion  of  the  density  assumption,  Eq.  (4.  8),  and  integration  gives 


p  r  3 

_  p  r _ 


(4«  12) 


3  q  +  3 

where  use  is  made  of  Eq.  (4.  8)  again  after  integration.  Recognizing  that 
r^  remains  fixed  during  differentiation  following  the  motion,  we  form  the 
derivative  of  Eq.  (4.  12)  with  respect  to  time  along  a  particle  trajectory. 
The  symbol  D/Dt  denotes  the  operator  9/9t  +  u  9/9r: 


0  =  ^ _ ^  +  p 

^  q+3  Dt  P 


from  which 


2b  c  p  djn  (q  __t_3]. 

Dt  P  dt 


3ur 
q  +  3 


3pu 

r 


(q  +  3) 


2  ar  ^ 


{4,  13) 


where  the  particle  velocity  u  =  Dr/Dt. 

The  equation  of  continuity  for  the  fluid  region  engulfed  by  the 

explosion  is 

^  +  p  ^  +  iei .  0. 

Dt  ^  or  r 


(4.  14) 
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Inserting  the  result  of  Eq.  (4.  13)  into  Eqc  (4. 
terms 


d  In  (q  +  3) 
3t  ~ 


we  obtain  after  rearranging 


(4,  15) 


The  right-hand  side  of  Eq,  (4,  15)  is  a  function  of  time  only  and  integration 
can  be  performed  with  respect  to  distance  at  fixed  time  as  follows: 

R  R 

r  r 

and  since  u  »  u  at  r  =  R, 
s 


r 

R 


(1  -  In  ~)  , 


(4.  16) 


in  which 


a  (t)  = 


R  d  In  (q  +  3) 
u  3t 


din  (p^  - 
“dTInR 


where  dR/ dt  is  the  shock  velocity  U,  and  use  has  been  made  of  EqSo  (4®  1) 
and  (4®  10)®  Now  the  over  density  p^  ”  ^  slowly  varying  function 

of  t  (and  hence  R)  in  the  strong  shock  region;  we  neglect  a  In  r/R  in  compari- 
son  to  unity  in  Eq®  (4o  16)  to  obtain  approximately 


r 

R 


0  <  r  <  R, 


(4®  17) 


and  the  velocity  profile  in  the  engulfed  region  is  nearly  lineaxo 

Do  Pressure  Profile 

The  equation  for  radial  motion  of  an, ideal  fluid, in  the  absence  of  body 

forces  is 


(4.  18) 
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where  P(r,  t)  is  the  absolute  pressure  in  the  engulfed  region  0  ^ 


The  derivative  of  Eq, 


following  the  motion  is 


d  u  u  U  u  ^  1 

s  s  s 


s  s  _r_ 

R  R  R 


r  <  R. 


(4.  19) 


where  use  is  made  of  Eq.  (4.  17)  again  after  differentiation,  and  dR/ dt  is 
the  shock  velocity  U. 

Upon  insertion  of  the  density  assumption  and  Eq.  (4,  19)  into  Eq.  (4.  18), 
integration  can  be  performed  with  respect  to  distance  at  fixed  time  as  follows: 


—  dr  -  p 
9r  ® 


■  XT  J 

u  U  du  u 

s  s  s / 

“R~  ■  dt  "  R  J 

.  J  d 


(i)  dr. 


since  the  quantity  in  brackets  is  a  function  of  time  only.  The  result  of  inte¬ 


gration  is 


P(r,t: 


fu  U  du  u  9+2 

*  S  S  S  R  fT 

R - at - R“  iT2  ^R  ^ 


(4.  20) 


Evaluating  Eq,  (4,  20)  at  r  =  R  and  noting  that  P(R,  t)  =  Pg(t)>  we  obtain 


p^(t)  -  PC0,t)  =  p. 


u  U  du^  u 
s  s  s 


s  R  ”  dt  R  q+2 


(4.  21) 


Dividing  Eq.  (4. 


Eq,  (4.  21)  results  in 


P(r,t)  -  PC0,t)  _  ,r  , 
p^(t)  -  P(o,tT" 


(4.  22) 


From  symmetry  of  the  spherical  explosion,  the  fluid  at  the  origin  must  be 
stationary.  In  view  of  the  assumption  for  density,  Eq,  (4,8),  we  assume 
P(0,t)  =  0  and  obtain  from  Eq,  (4,  22)  the  approximate  pressure  profile 


P  =  P 


s  'R' 


0  <  r  <  R  . 


(4.  23) 
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Vo 


ENERGY  BALANCE 


We  note  a  further  geometric  interpretation  of  the  explosion  energy- 
partition  in  Figo  lo  The  chord  of  the  Hugoniot  divides  the  rectangle  repre¬ 
senting  the  total  shock  energy  input  to  a  unit  mass  of  material  into  two  equal 

triangles,  the  upper  one  of  which  can  be  considered  as  representing  the 

2 

kinetic  energy  per  unit  mass  at  the  shock  front,  u  /2,  in  accordance  with 
Eq,  (4o  4)o  The  lower  triangle  then  represents  the  internal  energy  increment. 
If,  folio-wing  the  discussion  of  Section  III-D,  we  consider  the  expansion  of  the 
shocked  unit  mass  as  being  equivalent  to  expansion  along  the  Hugoniot  to 
pre- shock  conditions,  the  recoverable  fraction  of  the  internal  energy  incre¬ 
ment  is  given  by  the  area  under  the  Hugoniot  in  Fig»  1.  Then  the  lens- 
shaped  area  between  the  Hugoniot  and  its  chord  may  be  regarded  as  waste 
heat,  or  energy  per  unit  mass  made  imrecoverable  as  mechanical  work  as 
a  result  of  the  shock  process, 

A,  Unavailable  Energy 

When  an  amount  of  energy  W  is  released  into  the  medium  at  a 
point,  the  division  of  energy  into  waste  heat  and  kinetic  and  available 
mechanical  (or  hydrodynamic)  energies  when  the  spherical  shock  front  has 
radius  R  is  described  by  the  equation 

R  3 

w  -  4  IT  J  Q  R^  dR  =  i-3-^  (E^y  +  ’Sj^y) .  (5.  1) 

o 

In  this  equation  «^d  represent  average  hydrodynamic  and 

kinetic  energies  per  unit  volume,  respectively,  the  averages  being  taken 
over  the  entire  shock- engulfed  volume.  The  waste  heat  per  unit  volume 
is  denoted  by  Qo  Following  the  discussion  of  Section  III-G,  the  total 


-16 


Report  No. 
ARF  4132-6 


unavailable  energy  is  calculated  as  an  integral  over  the  previous  history 
of  the  shock  front,  since  that  energy  is  considered  as  commited  to  the  un¬ 
recoverable  state  upon  passage  of  the  shock  through  the  naaterial.  We  note 
that  after  sufficient  time  all  of  the  explosive  energy  must  appear  in  the  form 
of  thermal  energy  as  in  a  calorimeter  experiment,  so  that 

GO 

W=  4  -»r  JoR^dlL  , 

o 

The  quantity  Q  is  evaluated  as  the  lens-shaped  area  in  Fig.  1  multiplied 
by  the  ambient  density  ,  because  a  spherical  shock  wave  of  radius  R 
must  advance  a  distance  numerically  equal  to  V^/4irR^  in  order  to  engulf 
an  additional  unit  mass  of  material.  Therefore 

PdV.  (5. 3) 


B.  Kinetic  Energy  Profile 

The  kinetic  energy  of  radial  motion  per  unit  volume  of  the  engulfed 
fluid,  denoted  by  ,  can  be  written  using  the  density  and  velocity  distri¬ 

butions,  Eqs.  (4.8)  and  (4.  17),  as 


"KV 


1  2 
=  2  = 


^  u 
s  s 


2  (L.) 


q+2 


=  E 


KVs 


q+2 


(5.4) 


where  E„.,  denotes  the  kinetic  energy  per  unit  volume  at  the  shock  front, 
KVs 

The  average  kinetic  energy  density  E^^y  appearing  in  Eq,  (5.  1)  is  calculated 
from 


R 


ir  R- 

-j- 


"KV 


=  4  IT  E 


KVs 


r  2 

(|)  dr. 


o 
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ajid  the  result  is 


^KV 


(5.5) 


C.  Hydrodynamic  Energy  Profile 

The  recoverable  mechanical  (or  hydrodynamic)  energy  contained  in 
a  unit  mass  of  material  immediately  following  shock  is,  as  shown  in 
Appendix  B, 

^Hs 

^s 

where  the  integral  is  understood  to  be  taken  along  the  Hugonioto  Then  the 
hydrodynamic  energy  per  unit  volume  at  the  shock  front  is  gi^yen  by 


p  £ 


(5.  7) 


Assuming  the  expansion  to  take  place  along  the  Hugoniot,  one  could  calculate 
the  total  recoverable  energy  in  the  engulfed  volume  from 


where  is  the  average  hydrodynamic  energy  density  appearing  in  the 

rl  V 

energy  balance,  Eq«  (5.  l)o 

For  computational  purposes,  it  is  desirable  to  express  as 

the  corresponding  energy  density  at  the  shock  front,  Epjyg  >  ^nnltiplied  by 

a  form  factor.  We  observe  that  is  already  given  by  Eqo  (5p  5)  in 

xv  V 

such  a  form,  but  that  this  will  not  necessarily  result  for  if  the 

operations  of  Eq»  (5.8)  are  carried  out  after  insertion  of  the  assumed 
Hugoniot,  Eq.  (4,  7)  and  the  density  profile,  Eq.  (4.8),  into  the  first  inte¬ 
grand. 
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If  the  area  under  the  Hugoniot  in  Fig.  1  is  considered  as  a  function 
of  Pg  and  plotted  on  logarithmic  coordinates  for  several  values  of  the 
parameters  p  and  n  ,  one  observes  that  the  relationship  is  quite  closely 

approximated  by 

E  -Jpj. 

®Hs  -  ^s  ’ 

where  j  is  constant  or  at  most  a  slowly  varying  function  of  P3  (and  hence 

of  V  ).  Then  the  coefficient  J  is  correspondingly  either  constant  or  a 
s' 

slowly  varying  function  of  Pg  • 

Combining  Eq.  (5.7)  with  Eq.  (5.9)  results  in 


"HVs 


=  J  Ps 


(5.  10) 


Taking  logarithms,  one  has 


=HV6  “  + 


Differentiation  yields 


d  In  Ejjyg  _  4  J.  -p  d  j 

d  In  P.  s  d  In 


,  ,  _  d  In 
.  d  InJ  .  _ _ 

>“  ^  d  In  P^'  d  In  Pg 

s  s  s 


(5.  11) 


Because  j  and  J  are  observed  to  be  at  most  slowly  varying  functions  of 

P  and  because  p.  is  also  slowly  varying  with  respect  to  in  the 

s  s 

present  problem,  we  have  approximately 


d  In  E, 


(5.  12) 


Now  it  is  assumed  that  throughout  the  engulfed  region  the  hydrodynamic 
energy  per  unit  volume  is  related  to  the  pressure  in  the  same  manner  as 


Eq.  (5.  10)  for  the  shock  front  values.  That  is, 


%v  =  J  p- 


(5.  13) 
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and  j  is  calculated  from  Eq,  (5,  12).  Thus,  using  the  density  and  pressure 


profiles,  Eqs.  (4.  8)  and 


'  ps  i 


j(q+2)  +  q 


”  ^HVs 


j(q+2)+q 


(5.  14) 


The  average  hydrodynamic  energy  density  appearing  in  the  energy 


balance,  Eq, 


Q  A.  f  ^ 


is  calculated  from 


4w  E. 


j(q  +  2)  +  q  2 
{^)  rdr, 


and  the  result  is 


'hv-^  ir^iT(T+  2r+‘l, 


(5-15) 


where  is  to  be  obtained  from  the  definitions  of  Eqs<,  (5«  6)  and  (5,  7). 

Do  Peak  Pres  sure  Distance  Relation 

We  momentarily  regard  all  variable  quantities  appearing  in  Eq«  (5.  1) 

as  functions  of  the  single  variable  R  (as  they  must  be  for  any  preassigned 

energy  release  W),  and  differentiate  to  obtain,  after  removing  the  common 
2 

factor  4Ti'R  , 


-Q  =.  R  gd_  ^  _ 

Now  adding  and  subtracting  (R/3)  dQ/dR  from  the  right-hand  side  and 
arranging  yields,-  the  differential  form 


IfEH-y-  -f  +  Q) 

(Ehy  +  +  Q) 


^HV  ^KV  *  Q 


.  (5.16) 


Report  No. 
ARE  4132-6 


A  partial  integration  can  now  be  performed  from  a  prespecified  starting 
condition  R  -  Rj  to  an  arbitrary  upper  limit.  This  results  in 


In  (is-)  =  -  In  — - - 

1  (^HV  ^  ^KV  ^ 


^  ^KV  ^  ^  r  do  ,, 

+  Erv  ■*■  Q^i  J  ®hv'*‘®kv'‘’ ^ 


where  the  subscript  1  has  been  used  to  denote  conditions  when  tl.  =  Rj  • 

All  quantities  on  the  right-hand  side  of  Eq*  (5»  17)  can  be  calculated 
in  terms  of  through  Eqs,  (4«4),  (4«  7),  (4#  10),  (5^3),  (5#  5),  (d*  6),  (5e  12) 
and  (5.  15).  We  also  write  dQ  =  jdQ  (V)/dV3dV,  so  that  Eq.  (5.  iV)  can  be 


rewritten 


^  ^KV  ^^1 

^1  "  (Ejj^  +-Er^  +  Q) 


exp  r  -  .  ^5.  18) 

J  E„„  +  +  Q 


.V  ,  ■  HV  KV 

*sl 


Noting  from  Eq.  (5.  12)  that 


=HVs  d  P,  -  -  V,  d 


p  y. 

s  s 


V 

J 


list  the  functions  appearing  on  the  right-hand  side  of  Eq.  (5.  18)  as  follows 


s  V 


PdV 


=  P  d  V  PTV 

1  _  _ !_  1  +  ^  .1 

V  d  P  '  ^  rV 


PdV 


,  (5.  19) 


(3  -  1)  +  1 


3P  V 
s  o 


(5.  20) 
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(5.  21) 


V  dP 
P  3? 


(5.  22) 


where  P^  is  given  by  the  assumed  Hugoniot,  Eq.  (4.7),  and 


(5.  23) 


It  will  be  noted  that  the  integration  indicated  on  the  right-hand  side  of 
Eq.  (5.  24)  can  be  performed  analytically  for  integral  values  of  n.  Functions 
of  V,  which  is  now  regarded  as  a  variable  of  integration,  appear  in  the  ihte- 
grand  of  Eq.  (5«  18),  and.  are  obtained  by  replacing  by  V  in  the  formulas 
of  Eqs.  (5.  19)  to  (5.  21). 

E.  Initial  Conditions 

To  obtain  starting  conditions  (denoted  by  the  subscript  1)  for  the 
evaluation  of  Eq.  (5.  18),  we  consider  the  basic  relation,  Eq.  (5.  1),  with 
the  follovi^ing  simplified  representation  of  the  behavior  at  early  times:  Take 
R  j  such  that  the  waste  heat  P®3r  unit  mass  at  this  point  (i.  e.  ,  the 

energy  per  \mit  mass  converted  to  thermal  energy)  is  just  sufficient  to  induce 
a  change  of  phase  when  the  material  is  decompressed  to  the  pre-shock  pressure. 
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For  R  >  Rj  ,  therefore,  the  material  undergoes  no  change  of  phase  and 
Eq.  (5.  18)  is  applicable.  Now  in  0  <  R  <  R^  assume  conservatively  that 
the  waste  heat  Q  per  unit  mass  is  no  higher  than  its  value  at  R  =  Rj. 
With  Q  constant  at  the  value  in  0  <  R  ^  Rj  ,  the  basic  relation, 

Eq.  (5.  1),  can  be  integrated  to  yield 

4-n-R.  ^ 

“3  ^  ^KV 


which  gives  a  relation  between  the  energy  yield  W  and  quantities  with  sub¬ 
script  1  in  Eqp  (5p  18)o  This  is  equivalent  to  replacing  the  point  release  of 
energy  by  a  sphere  of  uniform  temperature  and  radius  R  containing  the 
same  energy®  For  condensed  matter  (so-called'Water-like"  substances)  it 
is  shown  in  Appendix  B  that  Q  =  c  (T  -  *^0^’  where  c  is  the  specific 
heat  of  the  material  and  T  is  the  absolute  temperature*  To  obtain  starting 
conditions  for  an  explosion  in  water,  therefore,  one  would  find  from 

Qf  =  c  (Tj  -  T^),  where  T^  is  the  boiling  temperature  of  water  at 
normal  pressure;  from  the  Hugoniot  curve,  one  would  then  find  the  pressure 
and  specific  volume  at  which  the  lens -shaped  area  is  equal  to  the  value  of 

Q,  V  so  obtained® 

1  o 

Inserting  Eq®  (5o  25)  into  Eqo  (5o  18)  results  in 


r  V 


47rR'' 


3W  ^  ^ 


exp 


^HV  ^  ^KV  ^  ^ 


(dQ/dV)  dV 

J  ^HV  ^KV 


Q 


(5.  26) 


s  1 


which  can  now  be  integrated  numerically  from  the  starting  value  of  the  specific 
volume  Vg  j  to  a  series  of  values  of  to  yield  a  relationship  between 

shock  radius  and  specific  volume  for  a  particular  material.  In  view  of  the 
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assumed  Hugordotj,  Ecqo  |4o  7}^  this  implies  a  relationship  between  radius 
and  shock  pressure  which  represents  the  pressure^dietance  decay*, 

Fo  Shock  Arrival  Time 

When  numerical  results  are  obtained  .for  the  pressure'^  and  volume- 

distance  relationships ^  the  time  of  arrival  t  of  the  shock  as  a  function 

s 

of  shock  radius  ;R  can  b®  obtained  by  numerical  integratioxi  of 

R 

*s=  /  •  (5-27) 

o 


where  .  IJ  ib  obtained  .from^  Eq©  (4*  3)o  The  shock  velocity  ,  U  will  be  quite 
large  at  R  -  Rj, ,  sc-  that  its  reciprocal  can  be  considered  as  essentially 
zero  at  R  —  0  and  varying  linearly  in  ,0  R  ^  R^  for  purposes  of  the 
calculationo  This  assumption  cannot  seriously  affect  the  resulting  arrival 
tillages,  and  is  implicit  in  a  numerical  integration  utilizing  trapezoidal 
appr  oximations  o 


VIo  SOLUTIOMS  FOR  PEAK  PRESSURE  AND  ARRI.VAL  TIME 

The  operations  indicated  in  the  equations  for  shock  pressure  and 
arrival  time  developed,  in  Section.^  represent  for  any  set  of  values  of  C, 


n^  and  corresponding  to  a  pairticular  material j,  and  for  a  particular 

energy  release  Wj,  a  good  deal  of  routine  numerical  or  semigraphical 
computatiosio  These  equations  how  eve  r^  are  readily  adapted  £ot  mauchine 


calculaticono 

Non-Dimensional  Forms 

We  introduce  the  non-dimensional  notation 


B 


\ 

(6-  1) 


j 
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s 


X  -  ^ 

P  dV“ 

Then  Eq*  (5.  26)  is  written 


so  that 
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w 


I 

V 


V  “  1 


)  (6.4) 


y  (v  -  1) 


The  Hugoniot  relations  Eq.  (4,,  7),  becomes 


P  = 


n  , 

V  -  I 

J*-  -  V 


(6o  5) 


from  which 


n 


n  V  V 

^  "  V*"  =  1  ^  “  V 


(6.6) 


After  Eqo  (5.  24}  is  put  in  non-dimensional  form.,  the  indicated  integration 
can  be  performed  by  standard  methods  to  yields  for  integral  values  of  n, 


z  = 


X 


( 1 )  In 


(p°l)v 


V 


(F^v)v 


2,  2 

.P  (v  -1) 


In 


1 

jx-v 


In  v=p 


V 


(6.7) 


z  = 


pV-1) 


ITX  ^ 

'•^•l)lni^^  =  In 'V“|x  (1  -  Jfx^  — i:-  (v^-1) 


Jx-v 


V 


(n£3)o 


These  results  make  possible  the  numerical  evaluation  of  Eqo  (6.  2) 
by  machine  calculation  for  selected  values  of  p  and  the  integer  n,  once  the 
upper  limit  integral  in  Eqo  (6.  2}  is  specified.  The 

results,  in  conjunction  with  the  numerical  evaluation  of  Eq,  (6.,  5)  for  the 
same  values  of  P  and  n,  constitute  a  two-parameter  family  of  dimensionless 
pres  sure -distance  data. 
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We  now  define  the  non-dimensional  time  of  shock  arrival 


T  =  Vc^  ‘s  • 

Then  the  arrival  time  definition,  Eq#  (5#  27),  can  be  written  in  the  following 
dimensionless  form  using  Eq.  (4*  3)  and  the  non-dimensional  notation  of 
Eqs,  (6.  1): 

7  = 

o 

The  early  time  history  (0  <  R  <  Rj)  pf  the  growth  of  the  explosion  is 
accounted  for  according  to  the  assumption  discussed  in  Section  V-F>  That 
is,  1/U  and  its  non-dimensional  equivalent are  taken  as  zero  at 
R  =  0  and  varying  linearly  in  O^R^Rj*  The  numerical  integration  of 
Eq*  (6.  9)  using  the  dimensionless  pressure  decay  data  obtained  by  evaluation 
of  Eq*  (6.  2)  is  now  straightforward,  and  can  easily  be  performed  by 

machine  calculation. 

B.  Numerical  Integration 

Programs  were  written  for  the  numerical  evaluation  of  the  specific 

volume -distance  relation,  Eq.  (6#  2),  and  the  non-dimensional  arrival  time, 

Eq.  (6.  9),  on  an  IBM  650  digital  computer.  As  an  initial  condition  for  the 

integration  of  Eq.  (6.  2)  the  upper  limit  v^  =  Tft/S  was  chosen.  The 

evaluation  was  performed  for  the  values  p  =  4,  5,  6,7,  and  8,  and  for 

integral  values  of  n  from  one  to  ten*  This  yielded  a  two-parameter  family 

of  50  data  curves  each  for  pres  sure -distance  decay  and  arrival  time* 

The  upper  limit  Vj  7|>x/8  corresponds  to  a  starting  specific  volume 

V  =  8V  /7|Jt,  and  was  selected  as  a  value  sufficiently  high  so  that  all 
s  1  o 

information  desired  in  any  application  of  the  results  is  likely  to  be  associated 


pv 


dX. 


(6.9) 
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with  values  of  v  less  than  Properties  of  the  results  discussed  in  a 

later  section  make  this  rather  arbitrary  selection  unimportant,  and  the 
results  turn  out  to  be  applicable  with  simple  modification  to  any  other 
starting  conditiono 

The  values  of  M  and  n  for  which  the  calculations  were  performed 
were  selected  from  experience  with  rea.1  materials  as  a  set  covering  a  wide 
range  of  possible  valueSo  The  selection  of  and  n  for  any  real  material 
of  interest  is  discussed  in  Section  VII  of  this  reports  In  applications, 
pressure-distance  and  arrival  time  data  for  any  specified  values  of  jx  and 
n  lying  within  the  range  covered  .can  be  obtained  by  straightforward  inter¬ 
polation  between  the  calculated  curveSo  The  param.eter  C  is  contained,  in 
the  non-dimensional  variables  and  does  not  appear  explicitly  in  the 
calculationSo 

lo  •  Pr e s sure ° Distance  Solutions 

The  program  for  evaluation  by  machine  calculation  of  the  specific 
volume  “distance  relationship,  Eqo  (6o  2),  and  the  associated  intermediate 
functions  of  v,  was  written  in  conventional  computer  languagCo  The  density 
ratio  V  was  regarded  as  the  independent  variable  for  purposes  of  the  step¬ 
wise  calculation,  a  uniform  increment  of  Oo  1  in  v  being  used  for  all  cases. 

Standard  sub-routines  were  used  for  the  evaluation  of  logarithms 
appearing  in  the  expressions  for  z,  EqSe  (6o  7|o  The  numerical  integration 
of  the  exponent  in  Eq,  (6o  2)  was  performed  using  Simpson’s  rule  to  approxi¬ 
mate  the  integrand  over  the  increments  of  v.  The  integration  proceeded 
downward  from  v  ~  7p-/8,  the  terminal  point  being  in  the  range  1,  025  <  v:^l. 

For  each  value  of  v,  the  quantities  p,  ph,  pk,  and  pw  were 
printed  out.  The  latter  three  are  seen  from  their  definitions  in  Eqs,  (6,  1) 
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to  be  the  non-dimensional  hydrodynamic,  kinetic,  and nnavailable  energies 
per  unit  voltime.  In  addition,  the  values  of  the  integrand,  the  integral,  and 
the  cube  root  of  the  coefficient  of  the  exponential  in  Eq.  (6»  2)  were  also 
printed  out  on  computer  data  sheets. 

2.  Arrival  Time  Solutions 

The  machine  program  for  calculation  of  the  non-dimensional  tim.e 
of  shock  arrival  by  Eq.  (6,  9)  was  also  written  in  conventional  language. 

The  numerical  integration  was  performed  using  a  trapezoidal  rule  to  approxl 
mate  the  integrand  over  the  increments  of  X.  These  increments  were  nou 
uniform,  but  were  determined  by  the  imiform  increment  in  v  of  0.  1  and 
taken  from  the  results  of  the  preceding  calculations. 

In  accordance  with  the  discussion  of  Section  V-F,  the  integrand  was 
assumed  to  be  zero  for  X  =  0,  and  the  integration  proceeded  upward  from 
X  =  0.  Implicit  in  the  integration  by  a  trapezoidal  rule  is  the  linear  approxi¬ 
mation  for  the  integrand  in  0  ^  X  ^  Xj,  where  X  -  Xj  corresponds  to 
V  =  Vj. 

For  each  value  of  X,  the  values  of  T  and  the  integrand  were 

printed  out  on  computer  data  sheets. 

C.  Graphical  and  Numerical  Results 

The  shock  pressure  and  arrival  time  data  derived  numerically  by 

the  procedures  described  earlier  are  summarized  in  Tables  1  through  20. 

We  define  X  as  the  cube  root  of  the  coefficient  of  the  integral  in 
e 

Eq,  (6.  2): 

,  -1/3 

Xe  =  ^(h  +  k  +  w)j  (6.  10) 

This  quantity  is  referred  to  as  the  "envelope  radius".  It  is  evident  from 
inspection  of  Eq,  (6.  2)  that  if  an  initial  condition  v  other  than  v^  is 
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assumed,  the  applicable  pressure-distaince  data  for  some  particular  values 
of  p  and  n  would  begin  at  X  =  X  for  those  values  of  p  and  n«  Hence 
X  regarded  as  a  function  of  v  for  given  values  of  P  and  n  represents  a 
locus  of  initial  conditions  on  the  pressure-distance  data*  The  values  of  X^ 
are  included  in  Tables  1  through  20*  The  application  of  this  information  to 
account  for  the  effect  of  an  initial  condition  differing  from  that  used  in  the 
machine  calculations  is  discussed  in  the  following  section* 

The  non-dimensional  shock  pressure  p  is  plotted  in  logarithmic 
coordinates  as  a  function  of  the  non-dimensional  shock  front  radius  X  for 
selected  values  of  p  and  n  in  Figs.  3  through  7*  The  values  of  X  are 
cdso  plotted  as  the  dash  curves  in  these  figures,  and  these  curves  of  X^ 
are  referred  to  as  ''envelope''  curves  because  each  represents  a  locus  of 
initial  pressure-distance  conditions. 

The  non-dimensional  arrival  time  T  is  plotted  in  logarithmic 
coordinates  as  a  function  of  X  for  selected  values  of  p  and  n  in  Figs.  8 
through  12. 

D.  Properties  of  the  Solutions 

The  solutions  for  shock  pressure  and  arrival  time  as  fiinctions  of 
distance  from  the  explosion  source  have  interesting  properties  with  regard 
to  the  effect  of  the  selection  of  initial  conditions.  Close  in  to  the  explosion 
source,  the  pressure  decays  are  quite  rapid  with  distance,  these  rates  of 
decay  decreasing  rapidly  beyond  the  region  of  effective  shock  absorption, 

1,  Variation  of  Initial  Conditions 

An  arbitrary  selection  of  the  initial  condition  Vj  =7 p/8  has  been 

made  to  permit  solution  in  non-dimensional  form  of  the  equations  governing 

« 

shock  pressure  decay  and  arrival  time.  This  corresponds  to  a  starting 


-30- 


Report  No. 
ARF  4132-6 


value  of  the  specific  volume  which  is  in  general  not  realized  in  an  appli¬ 
cation  to  an  actual  situation.  The  effect  of  this  difference  is,  however, 
easily  accounted  for. 

Consider  a  material  for  which  the  starting  value  of  the  specific 
volume  corresponds  to  a  value  of  v  other  than  Vj.  Denote  this  value 
by  V,.  The  actual  non-dimensional  distance  would  be  given  in  terms 

of  V  by  Eq.  (6.  2)  with  the  upper  limit  of  integration  w 2; 


We  may  compare  the  values  X  and  X^  given  by  Eqs.  (6.  2)  and  (6.  11), 
respectively,  at  the  same  dimensionless  pressure  p.  Since  p  and  v 
are  uniquely  related  by  Eq.  (6.  5)  for  fixed  values  of  jJ-  and  n,  this 
amounts  to  comparison  at  the  same  value  of  shock  compression  or  density 
ratio  V.  The  envelope  radius  X^  =  j^p(h  +  k  +  w)j  ^  is  thus  the  same 
in  both  cases.  Dividing  one  equation  by  the  other  leads  to 


dv  .  (60  12) 

The  value  of  the  exponential  factor  in  this  equation  is  a  function  of  the  limits 
of  integration  v^  and  v^  only,  and  is  thus  fixed  as  far  as  the  present 
comparison  is  concerned. 

Taking  the  logarithm  of  Eq.  {6.  12),  we  obtain 
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In  logarithmic  coordinates  then,  the  effect  of  an  initial  condition  differing 
from  the  value  of  v^^  used  in  the  present  calculations  is  represented  by- 
translation  of  the  pressure-distance  curve  in  the  direction  of  the  X-axis* 
Each  envelope  curve  represents  a  locus  of  initial  conditions  for  given  values 
of  JA  and  n,  and  so  the  dista;nce  of  translation  is  equal  to  the  horizontal 
distance  between  the  pressure-distance  curve  and,.its  envelope  at  p  =  p^ 
corresponding  to  v  ^"^2®  translation  is  to  the  right  in  Figs*  3 

through  7  if  ^2  ^  and  the  translated  curve  is  meaningful  only  for 
p  <  p^o  The  value  v^  -  7|Ht/8  has  been  selected  as  a  value  high  enough 
so  that  v^  ^  most  cases  of  practical  interest*  In  the  event  ^  Vj 

for  some  application,  both  the  curve  and.its  envelope  may  be  extrapolated 
in  the  direction  .of  increasing  pressure,  and  a  translation  to  the  left  must 
then  be  performed* 

The  operation  of  translation  is  shown,  schematically  in- Fig*  ISa* 

In  general,  the  values  of  p  and  n  .for  any  application  will  differ  from 
the  particular  values. for  which  the  present  calculations  were  performed* 
Strictly  spe adding,  the  necessary  interpolations  for  both  the  curve  and  its 
envelope  should  be  accomplished  before  the  translation,  although  the  order 
of  operations  makes  little  or  no  difference* 

The  effect  of  an  initial  condition  V2  differing  from  the  value  Vj 
used  in  the  present  calculations  is  represented  also  by  a  translation  of  the 
arrival  time  curve*  The  actual  non-dimensional  time  of  shock  arrival  T 

a 

■would  be  given  b'y;  Eq«  (b*  9)  with  the  integration  following  the  shock  front  in 
the  actual  case ; 


(6.  13) 
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A  simplifying  assumption  on  the  early  history  of  the  explosion  similar  to 
that  discussed  in  Section  V-F  is  implied.  Now  define  K  such  that 


3K  = 


(6o  14) 


v^(h+k+w) 


Then  the  pressure-distance  curve  on  logarithmic  coordinates  is  translated 
a  distance  K  in  the  direction  of  the  X-axis.  Thus  X^  =  X  e  •  Recognizing 
that  the  values  of  p  and  v,  and  hence  of  the  integrand  of  Eq«  (6.  13), 
remain  unchanged  during  such  a  trajxslation  of  the  pressure-distance  curve. 


we  may  write 


f  d(Xe^) 

-  J 


<6.  15) 


where  the  non-dimensional  shock  velocity  4*  is  given  by 


Changing  the  variable  of  integration  to  X,  we  obtain  from  Eq»  (6.  15) 


(6,  16) 


T  =  e 
a 


vw 


^  K 
Te  , 


<6,  17) 


Taking  the  logarithm  of  Eq.  (6,  17)  results  in 


In  T  =  In  T  +  K  , 

a. 

which  indicates  that,  at  the  same  dimensionless  pressure  p,  the  logarithms 
of  the  actual  cind  calculated  arrival  times  differ  by  the  constant  K.  On  a 
logarithmic  arrival  time  chart  (Figs.  8  through  12),  a  point  representing 
a  given  non-dimensional  pressure  p  is  also  translated  a  distance  K  along 
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the  non-dimensional  distance  scale.  Ib  logarithmic  coordinates  then,  a 
time-of-arrival  curve  is  translated  a  distance  K  in  each  direction,  the 
translation  being  to  the  right  and  upward-ih  Figs,  8  through  12  if 
V2  Vjj  if  V2  ^  Vj,  the  trainslation  is  to  the  left  and  downward.  The 
translated  curve  is  meaningful  only  for  p  <  P20 

A  family  of  envelopes  to  the  time -of -arrival  curves  could  be 
constructed  with  the  non-dimensional  envelope  arrival  time  being 
given  by  T  -X  y24*(X  )  consistent  with  the  assumption  that  1/U  varies 
linearly  in  0  ^  R  ^  Rj.  The  information  would  however  be  redundant, 
since  for  any  case  of  interest  the  translation  distance  K  is  already 
determined  from  translation  of  the  pressure-distance  curve. 

The  operation  of  translation  of  a  time-of-arrival  curve  is  shown 
schematically  in  Fig,  13b.  Once  again,  any  necessary  interpolations  on 
p  and  n  should  be  accomplished  before  the  translation. 

2,  General  Properties 

From  inspection  of  the  pressure-distance  curves,  Figs,  3  through  7 

it  is  evident  that  the  pressure  decays  during  the  early  portion  of  the 

explosion  history  are  extremely  rapid.  The  pressure  in  .this  range  varies 
6>  T 

as  1/R^  or  1/R  «  The  corresponding  arrival  times  in  the  same  range 

2 

(Figs*  8. through  12)  appear  to  vary,  approximately  as  R  , 

Near  P/C  “  0*  1,  the  rate  of  pressure  decay  slows  quite  rapidly, 
as  is  evidenced  by  the  relatively  rapid  change  in  slope  of  the  pressure- 
distance  curves  near  this  .dimensionless  pressure*  At  lower  pressures, 
the  shock  energy  transmission  becomes  progressively  more  efficient  and 
the  pressures  appear  to  vary  at  rates  approaching  l/R  •  This  is  in  agree 
ment  with  acoustic  theories  of  wave  transmission  from  a  point  source* 
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An  acoustic  wave  transmission  solution  should  in  fact  be  an 
asymptotic  approximation  to  the  present  strong  shock  solution  fpr  large 
Ra  The  present  non- dimen sionaJ  formulation,  however,  becomes  highly 
indeterminate  at  v  ^  io  0,  and  numerical  solution  even  by  machine  calcula¬ 
tion  becomes  quite  difficult  for  values  of  v  close  to  1*  0. 

Vila  PROPERTIES  OF  MATERIA!^ 

The  present  approximate  theory  of  explosion  growth  in  condensed 
media  has  been  developed  in  such  a  way  as  to  avoid  the  formulation  of  a 
general  equation  of  state  for  these  materials,  and  replace  it  by  the 
assumption  of  a  generaliz;ed  Hugonict,  Eqo  (4o  7)o  This  assumed  form 
contains  three  arbitrary  parameters,  C,  H**  und  n,  the  determination 
of  which,  for  any  particular  material  of  interest,  permits  application  of 
the  present  generalized  pressure  and  arrival  time  data* 

The  assumed  form  for  the  Hugoniot  is  such  that  dP /dV  0  and 
d^P /dV^  >  0  for  all  values  of  P  0,  The  corresponding  model  of  the 
material  is  thus  a  solid  or  a  liquid  in  which  an  increase  in  hydrostatic 
pressure  results  in  a  decrease  in  specific  volume,  and  whose  compres¬ 
sibility  under  shock  decreases  with  increasing  shock  pressure*  Nearly 
all  real  materials  exhibit  this  kind  of  behavior  at  sufficiently  high  pressures. 
The  conditions  on  the  slope  and  curvature  of  the  pres  sure -specific  volume 
function  are  necessary  for  the  formation  of  shocks  from  finite -amplitude 
compression  wavesc 

The  propagation  of  a  large -amplitude  hydrostatic  pressure  disturbance 
in  a  solid  or  a  liquid  has  been  studied  in  the  present  report*  The  material 
has  been  assumed  to  be  fluid-like  in  behavior,  or  equivalently,  its  strength 
in  shear  has  been  assumed  negligible  in  comparison  to  the  pressures  of 
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interest,  A  longitudinal  wave  of  compression  in  a  real  solid  will  in 
general  contain  both  dilational  and  distortional  (or  bulk  and  shear) 

components.  When  the  amplitude  of  a  longitudinal  wave  is  sufficiently 
large  so  that  its  shear  component  cannot  be  supported  by  the  material, 
at  least  a  portion  of  the  energy  must  be  transmitted  as  a  wave  of  hydro¬ 
static  pressure  and  a  shock  wave  can  form.  In  the  material  model 
which  is  used  in  the  present  theory,  the  passage  of  the  wave  is  assumed 
to  destroy  the  strength  of  the  material  in  shear,  so  that  the  material 
is  indeed  fluid-like  in  behavior.  Thus  in  the  present  model,  no  energy 
is  transmitted  by  signals  which  could,  under  certain  conditions  in  a  real 
solid,  run  ahead  of  a  trailing  shock  wave.  The  assumptions  implicit  in 
the  fluid  dynamic  model  are  discussed  further  in  Appendix  A, 

A.  Non-Porous  Materials 

The  best  values  of  C,  and  n  for  the  present  description  of 
3J^y  given  medium  can  be  obtained  by  a  process  of  fitting  to  an  experi¬ 
mentally  derived  Hugoniot  for  the  particular  material  of  interest. 

Standard  fitting  criteria,  such  as  that  of  the  minimum  sum  of  square 
deviations,  can  be  applied  to  numerical  laboratory  data^^  For 

example,  for  lead  C  “  2,  86  kilobars,  jx  =  4,  and  n  =  4®  64  provide  an 
excellent  fit  of  the  data  obtained  by  Walsh  et  al,  For  most  metals 
the  Walsh  data  indicate  that  the  exponent  n  should  be  about  5;  the 
exponent  for  water  is  close  to  ?• 

Alternative  to  the  use  of  conventional  fitting  criteria,  one  may 
choose  to  satisfy  certain  properties  of  compressibility,  such  as  the  limiting 
compression  ratio  and  the  initial  rate  of  change  of  specific  volume  with 
pressure®  The  presence  of  the  three  parameters  C,  1^*,  and  n  permits. 


-36» 


Report  No® 
ARF  4132-6 


in  general,  the  satisfaction  of  three  such  conditions.  In  what  follows, 

we  present  methods  of  deducing  values  of  these  parameters  based  on  the 

satisfaction  of  general  properties  of  compressibility  data. 

The  quantity  M-  represents  a  limiting  value  of  the  density  ratio 

V  /V  when  the  shock  pressure  P  tends  to  infinity.  As  such,  it  plays 
o 

the  same  role  as  the  factor  (■Y+ 1)/(‘V- 1)  in  the  Hugoniot  for  a  perfect 
gas,  where  y  is  the  ratio  of  specific  heats.  An  appropriate  value  for 
this  quantity  in  any  given  problem  can  readily  be  inferred  from  shock 
compression  data.  For  most  solids  a  value  of  p  about  4  appears  appli¬ 
cable  the  influence  of  the  exact  value  of  P  on  the  final  results  is  of 
course  not  particularly  strong,  since  C  and  n  remain  to  be  selected 
so  as  to  provide  a  best  fit  of  compression  data. 

With  an  assumption  for  the  value  of  P,  two  further  conditions 
on  compressibility  can  be  satisfied  by  proper  choice  of  C  and  n.  One 
scheme  of  presentation  of  experimental  hydrostatic  compressibility 

data  at  pressures  up  to  10^  bars,  at  least  for  static  results  on  relatively 

12  13 

incompressible  solids  ’  ,  is  in  the  form 


AV 

--y-T 


AP  -  BP^  , 


(7.  1) 


o 

where  AV  is  defined  as  V  “  empirical  constants  A  and  B 

are  tabulated  for  a  number  of  solids. 

The  constants  A  and  B  can  be  related  to  C,  P,  andn;  C  and 
n  can  then  be  written  as  functions  of  A,  B,  and  the  assumed  value  of 
p.  Dropping,  the  subscript  s  on  P  and  V  for  brevity  in  the  assumed 
Hugoniot,  Eq.  (4.  7),  we  define  €  =  AV/  and  write  Eq,  (4,  7)  in  the 
form 
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(7.  Z) 


The  right-hctnd  side  of  Eq.  (7.  Z)  can  be  expressed  as  a  Taylor's  series 
in  €  about  €  =  0: 

€  +  o  •  •  (7»  3) 

From  Eq«  (7c  3)  we  observe  that  the  rate  of  change  of  P  with  €  at 
€  =  0  is  given  by 


pi-1  P 


(n-1) 


where  primes  on  f  denote  differentiation  with  respect  to  its  argument  €. 
These  identities  provide  relations  between  the  given  Taylor  coefficients  of 
Eq.  (7.  3)  and  those  of  Eq.  (7.  5);  the  required  inverse  series  is  therefore 
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(7.6) 


AV  fi-l  P 

V  ■  n  C 
o 


fi(n+l)  -  (n-1). 


P, 

C 


•  •  • 


Now,  by  comparing  the  coefficients  of  like  powers  of  P  in  Eqs,  (7.  1)  and 
(7.6),  we  obtain  C  and  n  in  terms  of  A,  B,  and  p  as  follows: 


C  = 


(P-1)^  A 


2(|x-l)  B  -  (p+1)  A 


2(p-l)  B  -  (p+l)A 

n  =  - 2 - • 

(p-1)  A 


)  (7.  7) 

J 


Typically,  A  is  of  the  order  of  10  ^  bars  (1  bar  =  10  dynes /cm  ) 
and  B  is  of  the  order  of  10  bars  ^  for  the  less  compressible  homogeneous 
solids  Assuming  p  =  4.  0,  C  is  of  the  order  of  10  bars  and  n  is  of 
the  order  of  unity. 

When  limited  compressibility  data  is  available  in  numerical  form,  a 

satisfactory  procedure  for  the  determination  of  C  and  n  for  a  non-porous 

solid  essentially  consists  of  a  numerical  fitting  process.  The  product  Cn 

is  related  to  the  bulk  modulus  of  the  material  through  Eq.  (7.  4);  hence 

2 

Cn  =  (p-1)^.  (7.8) 


where  c^  is  the  velocity  of  sound  in  the  medium.  Equation  (7.  8)  thus 
provides  one  condition  on  C  and  n  if  the  value  of  p  is  assumed,  A 
second  condition  on  C  and  n  for  the  assumed  value  of  p  can  then  be 
obtained  by  inserting  a  pair  of  experimental  values  of  P  and  V  into  the 
analytical  expression  for  the  Hugoniot,  Eq.  (4.7), 

B,  Porous  Materials 

As  in  the  case  of  non-porous  materials,  the  best  values  of  C, 
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Ji,  and  n  for  a  given  porous  material  would  be  obtained  by  a  process  of 
fitting  to  an  experimentally  derived  Hugoniot. 

The  Hugoniot 

In  what  follows,  we  present  an  approximate  method  for  deducing 
the  behavior  of  a  porous  solid  under  shock  from  considerations  of  the 
behavior  of  its  gas  and  solid  fractions.  The  method  consists  of 

1.  Describing  the  behavior  of  the  solid  fraction  under  shock  by 
an  equation  of  the  same  form  as  the  assumed  Hugoniot,  Eq.  (4,  7),  in 
which  the  constants  C,  and  n  apply  to  the  solid  fraction  alone. 

2.  Describing  the  behavior  of  the  gas  fraction  under  shock  by 
the  usued  Hugoniot  derived  on  the  assumption  of  the  perfect  gas  law. 

The  ambient  pressure  is  considered  as  being  very  small  in  comparison  to 
the  pressures  of  interest. 

3.  Considering  the  gas  and  solid  fractions  as  being  shocked 
independently,  and  permitting  heat  transfer  from  the  gas  to  the  solid  at 
constant  pressure.  The  temperature  rise  in  the  pure  solid  due  to  shock 
compression  alone  is  assumed  to  be  insignificant,  and  the  gas  fraction 
is  again  assumed  to  behave  as  a  perfect  gas.  The  change  in  volume  of 
the  solid  due  to  a  rise  in  temperature  is  neglected. 

4.  Characterizing  the  resulting  composite  shock  pressure- 
specific  volume  function  by  the  form  of  the  assumed  Hugoniot,  Eq.  (4.  7), 
in  which  the  constants  C,  |a,,  and  n  now  refer  to  the  mixture. 

Let  the  subscripts  a  and  b  denote  quantities  relating  to  the  gas 
and  solid  fractions,  respectively,  and  denote  ambient  conditions  by  the 
subscript  zero.  For  brevity  we  omit  the  subscript  s  on  P  and  V  in 
the  assumed  Hugoniot,  Eq,  (4.  7).  We  introduce  the  following  notation; 
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aO 


'bO 


am 


V 


PbO 

T 


am 


T 

^a 

‘^b 


=  initial  specific  volume  of  porous  material 
=  initial  volume  fraction  of  gas 
=  a  ~  initial  volume  of  gas 

=  (1  -  a)  V  =  initial  volume  of  solid 

=  volume  of  gas  after  shock  and  preceding  heat  exchange 

=  volume  of  gas  after  shock  and  heat  exchange 

s  volume  of  solid  after  shock 

-  initial  density  of  gas  fraction 

=  initial  density  of  solid  fraction 

absolute  initial  temperature  of  porous  material 

=  absolute  temperature  of  gas  after  shock  and  preceding  heat 
exchange 

=  final  absolute  temperature  of  the  shocked  porojus  material 
specific  heat  of  gas  fraction  at  constant  pressure 
=  specific  heat  of  solid  fraction  at  constant  pressure 
^  gas  constant 


We  study  a  volume  of  the  porous  material  equal  to  its  specific  volume, 
so  that  BJid  V^Q  do  not  represent  specific  volumeSo  We  restate 

Eq.  (4*  7)  as  it  applies  to  the  solid  fraction: 


V. 


'bO 


(7.9) 


V, 


b 


where  the  constants 

If  the  gas  fraction  is  shocked  separately  to  the  pressure  P,  the 
volume  is  given  by  the  usual  Hugoniot: 


and 


refer  to'the  solid  fraction  alone. 
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(7.  10) 


V 

am 


y  +  1 

m. 


p  +  p 


p  + 


•y  +  1 

-^-TT 


9 


•where  y  is  the  ratio  of  specific  heats  of  the  gas.  We  neglect  the  ambient 
pressure  P^  in  Eq.  (7«  10),  and  obtain  approximately 


V 


am 


•y  +  1 
y  ~  1 


(7.  11) 


Asstiming  that  the  ideal  gas  equation  of  state  holds  for  the  gas  fraction, 
we  may  write 

^  V  _  =  p  o  a  V  R*  (T  -  T  ), 
am  '^aO  o  *  '  am  o' 


after  neglecting  the  product  P  V  „  in  comparison  to  the  product  PV 

o  au  ^  am 

The  factor  P_^q  a  represents  the  total  mass  of  gas  involved  and  R^j.  is 
the  difference  of  the  specific  heats  of  the  gas  at  constant  pressure  and  constant 
volume,  respectively.  T*he  temperature  of  the  gas  fraction  after  shock  and 
before  mixing  with  the  solid,  T  ,  is  then  approximately 


~  O  R^^  ‘t'+  1 

where  Eq.  {7.  11)  and  the  definition  V^q  -  ®-'^o  been  used.  The 


(7.  12) 


temperature  of  the  solid  fraction  after  shock  and  before  mixing  with. the 


yas  is  assumed  to  be  T  . 

o 


We  now  consider  the  solid  and  gas  as  being  mixed  at  the  shock  pressure 
P,  and  assume  that  the  mixture  reaches  thermal  equilibrium: 


aV  c  (T  -  T) 
o  a  am 


V  c.  (T  -  T  .). 
o  b  '  o 


(7.  13) 


The  factor  p^Q  (1  -  represents  the  total  mass  of  solid  involved. 
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Substituting  Eq.  (7.  12)  and  rearranging*  we  obtain 


o  a 


y  -  V 

Y+  1 


P  = 


PbO<'  -  PaO  ^ 


c^(T  -  T^).(7.  13) 


Now  c  and  c,  can  be  expected  to  be  of  roughly  the  same  order  of  magnitude, 
a  b 

and  the  ratio  Pg^Q^PbO  values  like  5  x  10 for  typical  mineral  solids  and 

air  at  atmospheric  pressure.  Consequently,  the  factor  in  brackets  on  the 
right-hand  side  of  Eq.  (7.  13)  is  very  nearly  equal  to  the  total  mass  of  gas 
and  solid  enclosed  in  the  volume  V^,  or  unity.  The  final  temperature  of  the 
mixture  is  then  given  approximately  by  the  expression 


aV  c  , 

-r  -  T  +  °  ^  ^  ^ 

“  A  T  -  --  ...  JL  H 


y  +  1 


(7.  14) 


If  the  volume  of  the  gas  fraction  is  taken  to  be  proportional  to  the 
absolute  temperature  at  constsuit  pressure,  we  may  write 


am 


T 


(7.  15) 


am 


Hence,  using  Eqs.  (7.  11),  (7.  12),  (7.  14),  and  the  definition  ^aO  -  aV^,  we 


have 


a  y  -  1 

\  ^  yTT 


T  +  ^ 
o  R, 


!a  y  -  1  p 
Ck  y  +  1 


T  + 
o 


R 


y  -  1 
y  +  1 


(7.  16) 


We  now  write  Eq.  (7.  9)  describing  the  solid  fraction  using  the 
definition  =  (1  -  a)  V^: 


a 


%  V  % 
(1  -  a)  (-^ 


Hi,  -  (1  -  a) 


V 


(7.  17) 


V, 


b 
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where  C^s  and  can  be  obtained  for  the  pure  solid  by  methods 
suggested  in.  the  preceding  sectiono 

In  accordance  with  the  assumption  that  the  volume  of  the  solid  is 
not  affected  by  temperature  changes ,  we  write  the  ratio  of  specific  volumes 


of  the  poroiis  material  before  and  after  shock  as 


Ob  .  UcT- 

o  o 


(7o  18) 


and  each  .term  on  the  right ‘-hand  .side  of  Eq©  (7©  18)  is  determined  as  a  function 
of  shock  pressure  through  Eqs©  (7©  16)  and. (7©  17)©  The  composite  Hugoniot 
may  now  be  obtained  numerically  or  by  graphical  addition  of  the  shock 
pres  sure -volume  curves  for  the  gas  and  solid  components©  Such  a  graphical 
a^ddition  is  shown  schematically  in  Fig©  14  for  a  0.  25  and  1©  4©  .  The 
neglection  of  in  the  Hugoniot  for  the  gas  fraction  is  evidenced  by  the  fact 

that  the  composite  curve  does  not  pass  through  the  point  (0,  The 

structure  of  the  solid  matrix  containing  the  gas  in  voids  has  of  course  been 
assumed  to  have  negligible  strength^  so  that  the  material  is  converted  to  a 
mixture  of  gas  and  solid  particles  at  pressures  which  are  negligible  in 
comparison  to  the  shock  pressure  of  interest©  The  resulting  composite 
Hugoniot  thus  has  a  flat  portion  in  .the  neighborhood  of  (Og 

The  condition  of  thermal  equilibrium,  in  the  mixture  immediately 
following  shockg  which  is  asf^umed  in  writing  Eq©  (7©  13)^  depends  on  the 
mixture  being  sufficiently  finely  divided  so  that  the  required  transfer  of  heat 
can. take  place©  The  large  heat  capacity  of  the  solid  thus  effectively  causes 
•  the  gas  fraction  to  be  adsorbed  into  the  solid  mass© 

It  remains  now  to  determine  the  values  of  G,  and  n  in  the  general 
analytical  form  ©f  the  Hugoniot g  Eq©  (4©  7)g.  for  a  best  fit  of  the  composite 
curve©  It  may  be  iivOted  that  a  r.equirement  of  matching  the  initial  slope  of 
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the  composite  curve  through  Eq©  (7®  4)  is  mapplicable  because  of  the  effect 
of  the  presen.ce  of  a  gas  fraction  as  evidenced. by  the  initial  flat  portion  of 
the  composite  curve* 

2*  Compressibility  Parameters 

Consider  the  limit  of  m  Eq*  (7o  16)  as  the  shock  pressure  P 

a  o 

becomes  infinite: 


V 

lim  =  a 

P^oo  ''o 

Similarly  from  Eq. 


X:-! 

y .  +  1 

(7.  17), 


(7.  19) 


lim  y 
P-^*"00 


(7.  20) 


The  limit  of  the  fraction  V/  V  is  then,  obtained  as  .the  sum  of  the  limits  for 

o 


the  gas  and  solid  fractions,  Eqs<.  (7.  19|  and  (7<,  20)  abovco  Now  the  factor 
PaO^^o  total  mass  of  gas  in  the  specific  volume  of  the  mixture. 

If  c  /c^  is  of  the  order  of  unity,  lim(V^./V^)  is  entirely  negligible  in 
comparison  to  limfV^/V^)  for  gas  fractions  up  to  about  0.,  85. 

4 

For  example,  for  a  s  0.85,  ®  !•  0»  “  1.4,  ^  ^  ^  * 

and  s  4.  0,  we  have  lim  011°  Thus  we  obtain  the  limiting 

density  ratio  for  the  porous  material 

The  largest  value  of  P-  for  which  numerical  results. have  been  obtained 
is  p.  s  8,  which  corresponds  .to  a  50  per  cent  gas  fraction  if  p^  -  4, 

Values  of  C  and  n  in  a  generalized  Hugoniot  for  the  porous 
material  can  now  be  obtained  by  a  fitting  process  consisting  of  inserting 
two  pairs  of  values  of  P  and  V  from  the  curve  obtained  .by  graphical 
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addition  of  the  gas  and  solid  pres  sure  “^vcdurne  relationships  (Figo  14)  into 
-the  generalized  Hugoniot,  Eqo  {4o7}s  with  jx  given  by  Eqo  (7c.  21)o 

Vnio  APPLICATIONS 

The  shock  pressure  decay,  and  arrival  time  results  presented  in  this 
report  are  intended  for  application  to  explosions  in  liquids  and  in  porous  and 
non-porous  solidso  The  results  are  given  in  dimensionless  form  in  order  to 
accommodate  a  variety  of  materials  having  generally  similar  compressibility 
characteristics  under  shocko  The  techniques  described  herein  have  been 
successfully  employed  in  predictions  on  the  growth  of  high-yield  explosions 
in  water  and  soilo  In  these  cases j,  values  of  the  constants  in  the  general 
form  of  the  Hugoniot  for  the  materials  were  obtained  by  fitting  to  experi¬ 
mentally  obtained  shock  compressibility  datao  The  table  below  gives  values 
which  have  been  used  for  water  and  a  sandy  soilo 


C(bars)  | 

1  n 

Water 

n 

1  9,000 

4 

] 

LJ _ 

Sandy  Soil 

1  12,900 

4 

i  3 

Ao  Behavior  of  Real_Mate rials 

Liquids  and  non-porous,  non  “Crystalline  solids  under  hydrostatic 
compression  can  be  expected  to  exhibit  continuously  decrea-sing  compres¬ 
sibility  with  increasing  pressureo  This  is  a  characteristic  of  the  assumed 
general  Hugoniot ^  Eqo  (4o  and  the  behavior  of  such  materials  is  readily 
represented  by  this  formo 

.  The  porous  solids  under  confined  compression  may  show  a  distinct 
‘’crushing  pressure”  at  which  the  structure  of  the  solid  matrix  containing 
the  gas  in  voids  fails  suddenly^  and  the  porous  material  is  converted  to  a 
mixture  of  solid  and  gas  pa^rticleso  This  results  in  a  rapid  decrease  of 
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specific  volume  with  increasing  pressure  at  crushing,  followed  by  a  slow 
decrease  of  specific  volume  with  increasing  pressure  when  crushing  is 
complete.  If  results  are  desired  in  a  shock  pressure  range  well  above  the 
crushing  strength  of  the  solid  matrix,  the  description  of  the  porous  solid  as 
an  idealized  mixture  is  still  adequate  for  the  analysis  of  strong  shock  wave 
propagation.  In  this  case  the  constants  ,  C,  jJ-,  and  n  must  be  obtained  from 
a  gross  fit  of  compressibility  data,  disregarding  details  of  the  data  at  low 
pressures.  In  particular,  one  can  no  longer  make  use  of  the  initial  rate  of 
change  of  pressure  with  specific  volume  as  a  condition  on  the  constants, 
since  this  rate  of  change  is  related  to  the  acoustic  velocity,  or  speed  of 
low-amplitude  pressure  waves  in  the  material.  The  sound  wave  may  be 
regarded  as  being  transmitted  as  a  low -amplitude  stress  wave  through  the 
solid  matrix  without  inducing  crushing  of  the  matrix  structure,  and  causing 
only  negligible  compression  of  the  gas  contained  in  the  voids.  One  would 
find  that  the  shock  velocity,  which  is  related  to  the  slope  of  the  chord  of  the 
actual  Hugoniot,  for  the  pressure  range  exceeding  the  crushing  strength  can 
in  fact  be  considerably  lower  than  the  sound  speed,  which  is  related  to  the 
initial  slope  of  the  actual  Hugoniot.  This  behavior  gives  rise  to  the  possibility 
of  an  acoustic  signal  running  ahead  of  a  slower  moving  strong  shock  from  an 
explosion  in  a  porous  material,  but  such  an  acoustic  signal  would  carry  a 
relatively  small  amount  of  energy  compared  to  the  trailing  main  shock. 

The  crushing  phenomenon  is  related  to  the  initial  flat  portion  of  the  composite 
Hugoniot  of  Fig.  14,  where  it  has  been  assumed  that  ambient  atmospheric 
pressure  and  the  crushing  strength  of  the  material  are  negligible  in  magnitude. 
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The  Ideal  Absorber 


Clearly,  the  ideal  absorber  would  be  one  which  maxi.mizes.  the 
waste  energy  per  unit  volume*,  An  upper  limit  on  the  waste  heat  QV 

o 

per  unit  mass  of  material  is,  from  Figo  1,  the  total  internal  energy  imparted 
to  the  unit  mass  by  the  shocko  This  would  be  realized  by  a  material  whose 
Hugoniot  exhibited  extreme  upward  concavity;  it  follows  that  in  the  ideal 
absorber,  half  the  shock  energy  imparted  to  a  unit  mass  of  material  is 
dissipated  as  waste  heato  It  is  likely  that  a  porous  solid  having  proper 
characteristics  of  grain. size  and  thermal  conductivity,  could  approach. this 
sort  of  behavioro  An  increase  in  porosity  is,  however,  accompanied  by  a 
decrease  in. average  density,  so  that  a  unit  mass  of  the  material  occupies 
an  increasing  volume®  A  large  bulk  .of  very  low  density  porous  material 
would  thus  be  necessary  to  attenuate  appreciably  the  shock  pressure  from 
an  .explosion®  .  In. design  for  explosion  containment,  therefore,  it  is  necessary 
to  effect  a  proper  balance  betwe.en  porosity  and  total  mass® 

-  In  general,  it  is  desirable  for  efficient  shock  attenuation  to  use 
materials  in  a  pressure  range  in  which  the  rate  of  decrease  of  pressure  with 
■distance  is  extremely  rapid,  i®  e®  ,  in  the  range  in  which, the  pressure“ 
distance  curves  of  Figs®  3  through  7  are  very  steep®  Although  the  present 
results  are  for  the  case  of  a  point*- source  explosion  in, an  infinite  medium, 
they  make  possible  the  evaluation  of  a  wide  variety  of  solids  for  their  relative 
blast  absorbing  characteristics® 
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Time 


Fig.  2  DISTANCE -TIME  RELATIONS 
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Fig.  4  PRESSURE -DISTANCE  CURVES,  =  5 
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Fig.  5  PRESSURE -DISTANCE  CURVES,  fx  =  6 
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Fig.  6  PRESSURE- DISTANCE  CURVES,  =  7 
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ARRIVAL  TIME  CURVES, 


ARRIVAL  TIME  CURVES 


ARRIVAL  TIME  CURVES, 


SHOCK  ARRIVAL  TIME 


Fig.  12  ARRIVAL  TIME  CURVES, 
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Pres  sure- Distance  Curve  b.  Arrival  Time  Curve 

Fig,  13  VARIATION  OF  INITIAL  CONDITIONS 


Shock  Pressure 


Fig.  14  HUGONIOT  FOR  A  POROUS  SOLID 
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TABLE  2 
Values  of  \ 
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TABLE  4 

Values  of  X 

e 


JJi  -  1 


\' 

] 

3 

9 

r, 

t 

f 

C) 

L;: 

i .  1 

8.41 

6.  53 

5.59 

4.97 

4.52 

4.16 

3.86 

3.60 

3.  34 

2.91 

1 

5.32 

4.05 

3.40 

2.  97 

2.64 

2.39 

2.17 

1.992 

1.831 

1.680 

1.  K  j 

4.07 

3.05 

2.51 

2.15 

1.882 

1.667 

1.488 

1.337 

1.206 

1.091 

1.  -iO.; 

3^36 

2.48 

2.  01 

1.692 

1.453 

1.263 

1.106 

0.974 

0.861 

0.763 

1 .  :oo 

2.  89 

2.  10 

1.680 

1.391 

1.174 

1.002 

0.861 

0.743 

0.644 

0.559 

1 . 6 00 

2,56 

1.835 

1.443 

1.176 

0.976 

0.818 

0.691 

0.585 

0.497 

0.423 

1  .  700 

2.30 

1,  630 

1.264 

1.015 

0.828 

0.683 

0.566 

0.471 

0.  392 

0.  327 

1  .  soo 

2.10 

1.467 

1.122 

0.888 

0.714 

0.578 

0.471 

0.385 

0.  315 

0.258 

i .  900 

1.928 

1.  333 

1.007 

0.785 

0.622 

0.496 

0.397 

0.319 

0.257 

0.207 

.  0  00 

1.785 

1.221 

0.  911 

0.701 

0.547 

0.429 

0.  339 

0.268 

0.212 

0.1677 

0.  i  30 

1.663 

1.126 

0.  830 

0.  630 

0.484 

0.375 

0.291 

0.226 

0. 1764 

0. 1374 

.... 

1.555 

j 

1.042 

1 

0.759 

0.569 

0.432 

0.329 

i 

0.252 

0. 1931 

0. 1481 

0.1137 

0.  300. 

1.459 

0.968 

0.698 

1  0.516 

0.387 

0.291 

0.219 

0.1657 

0. 1253 

0.0948 

•  L  9  00 

1.372 

0.  902 

0.643 

'  0.470 

0.348 

0.258 

0. 1920 

0.1431 

0.1067 

0.0796 

0.  500 

1.293 

0.  843 

0.594 

0.430 

0.  314 

0.230 

0.1688 

0. 1241 

0;0913 

0.0672 

7 .  OOO 

1.219 

0.788 

0.550 

0.  393 

0.284 

0.205 

0.1489 

!  0.1080 

0.0785 

1 

0.0570 

i:.  rOi.i 

1.151  1 

i 

0.738 

0.510 

0.361 

0.257 

0.1838 

0. 1317 

0.0944 

0.0677 

0.0486 

Hoo 

! 

1.087 

0.691 

0.473 

0.331 

0.233 

0. 1648 

0. 1167 

0.0827 

0.0586 

i 

0.0415 

cT  07  C 

1.025  1 

0.647 

0.439 

0.304 

0.212 

0. 1480 

0.1036 

0.0725 

0.0580 

0.0360 

3.  OQO 

0.966 

0.605 

0.407 

0.279 

0. 1922 

0. 1329 

0.0920 

0.0637 

0.0441 

0.0306 

100 

0.909 

0.565 

0.  376 

0.255 

0.1744 

0.1193 

0.0817 

0.0560 

0.0383 

0.0263 

L  .:oo  , 

0.853 

0.526 

0.  348 

0.234 

0. 1579 

0.1069 

0.0725 

0.0491 

0.0333 

0.0226 

3.  300 

0.797 

0.488 

0.320 

0.213 

0. 1425 

0.0955 

0.0641 

0.0430 

0.0289 

0.01939 

3.  iOo 

0.740 

0.450 

0.293 

0.1930 

0.1280 

0.0850 

0.0565 

0.0375 

0.0249 

0.01658 

' 3. StO 

6.682 

' 

0.412 

0.265 

0, 1736 

0.1141 

0.0750 

0.0494 

0.0325 

0.0214 

0.01409 
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Values  of  p 


1 

> 

■1 

J  • 

1  ' 

•  ,  ■  '■ 

0.0191 

0.0396 

0.0617 

0.0854 

0. 1109 

0. 1384 

0. 1679 

0. 1996 

0.234 

0.270 

1  .  i  . 

0,0457 

0.0995 

0. 1626 

0.237 

0.324 

0.427 

0.547 

0.688 

0,855 

1.050 

i 

0.0738 

0. 1679 

0.288 

0.441 

0.636 

0.885 

1.202 

1.606 

2.  12 

2.78 

]  .  3  r  • 

0. 1034 

0.245 

0.441  . 

' 

0.710 

1,080 

1.588 

2.29 

3.25 

4.57 

6.  38 

3  .  i ,  5 

0, 1347 

0.333 

0.626  ■ 

1.059 

1.697 

2.64 

4.03 

6.07 

9.09 

13.54 

1  .  •  >  ■  S 

0. 1678 

0.432 

0.848 

1,505 

2.54 

4.17 

6,73 

10.76 

17,  12 

27.1 

1  .  L  5 

0.203 

0.543 

1.113 

2.07 

3.67 

6.34 

10.82 

18.  33 

30.9 

52.0 

1  .  "  3 

0,240 

0.667 

1.424 

2.77 

5.  15 

9.39 

16.90 

30.2 

53.9 

96.0 

1  .  -  -  r3 

0.280 

0.805 

1.789 

3.64 

7. 10 

• 

13.58 

25.8 

48.6 

91.3 

171.5 

1 , 9  /  3 

0.  322 

0,959 

2.22 

4.70 

9.60 

19.29 

38.4 

76.2 

150.8 

298 

V  5 

0.367 

1. 130 

2,71 

6.00 

12.  81 

27.0 

56.3  , 

117.2 

243 

506 

9.  175 

0.416 

1.321 

3.29 

7.57 

16.88 

37.1 

81.2 

176.9 

385 

838 

.:  7  5 

0.468 

1.532 

3.95 

9.46 

22.0 

50.5 

115.4 

263 

599 

1362 

37  5 

0.524 

1.768 

4.72 

11.74 

28.4 

68.0 

162.0 

385 

916 

2170 

-17  5 

0.5  84 

2.03 

5.61 

14.46 

36.4 

90.6 

225 

557 

1380 

3420 

7. .  /  S 

0.649  : 

2.  32 

6.63 

17.72 

46.  3 

119.8 

309 

797 

2050 

5280 

7: .  b  i’  5 

0.720 

2.  65 

7.80 

21.6 

58.5 

157.2 

421 

1127  1 

3120 

8068 

...  .  -  5 

0.798  1 

3.01 

9.16 

26.2 

73.5 

205 

569 

1580 

1  i 

4390  • 

I  ; 

12170 

i 

0.  882  1 

i 

3.42 

10.71 

31,7 

i 

92.0 

265 

764 

2200 

6310 

18160 

>-'5 

0.975 

3.88 

12.51 

38.2 

114.6 

342 

1 

1018 

3030  1 

1  9010 

1 

26800 

5 .  L  7  5 

1 

1,078 

4.  39 

14.58 

1  45.9 

142.3 

439 

1350 

4150 

j  12770 

39300 

3,  1  . 

1.192  i 

4.98 

16.99 

55.  1 

176.2 

561 

1782 

5-660 

1  17960 

1 

57000 

3 . 7  5 

1.319  I 

5.64 

19.78 

66.1 

218 

715 

2340 

7670 

25100 

82300 

3.  375 

1.462 

6.39 

23.0 

79.2 

269 

909 

i  3070 

10360 

35000 

118000 

3.475 

1.623 

7.26 

26.9 

95.0 

332 

1154 

i 

4010  ! 

13940 

48500 

168400 

3  .  5  /  .5 

1.807 

8.27 

31.3 

113.9 

409 

1464  j 

5240 

18720 

67000 

240000 

3 . 5  /  5 

2.02 

9.44 

36.7 

136.9 

1 

505 

1858 

6830 

25100 

92300 

339000 

3.  77  5 

2.27 

10.82 

43,1 

164.9 

625 

2360 

8920 

23670 

127100 

480000 

3.  S7.5 

2.56 

12.46 

50.8 

199.5 

776 

3010 

11660 

45200 

175100 

679000 

3.  97  3 

2.90 

14.44 

60.3 

243 

967 

3850 

15300 

60800 

242000 

961000 

•I  .  0  r'  5 

3.  32 

16.  87 

72.  1 

297 

1214 

4950 

20200 

82000 

335000 

4  .  1  /  3 

3.  85 

19.92 

87.0 

367 

1536 

6420 

26800 

11190C 

1  467000 

1 

■>  r., 

4.52 

23.8 

106,4 

459 

1968 

8420 

36000 

15390C 

1 

i  658000 

4  .  3  I  .3 

5.40 

29.0 

132.4 

585 

2560 

11220 

49100 

1 

1 _ 

21480C 

i  940000 

1 
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TABL.E  6 


Values  of  X 


/X-  5 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1  .  ovs 

1,17  5 

5.44 

3.23 

3.49 

2.06 

2.42 

1.421 

1.721 

1.004 

1.238 

0.718 

0.  897 
0.517 

0.650 

0.  374 

0.464 

0.271 

0.19^ 

1.  Z75 

1.375 

2.48 

2.09 

1,577 

1. 321 

1,080 

0,900 

0.759 

0.628 

0.539 

0,444 

0.  386 

0.  315 

0.278 

0.225 

0.200 

0.1614 

0. 1442 

0.  1158 

0. 1025 
0.0829 

1.  175 

1 .  57  5 

1.844 

1.671 

1.159 

1.045 

0.785 

0.704 

0.545 

0.485 

0.382 

0.339 

0.270 

0.238 

0.1915 
0. 1673 

0.1362 

0,1181 

0,0970 

0.0835 

0,0692 

0.0591 

1 . 67  5 

1 . 7  5 

1.541 

1,440 

0.960 

0.  892 

0.643 

0.594 

0,441 

0.405 

0.  305 
0.279 

0.213 

0. 1934 

0.1489 
0. 1343 

0. 1043 
0.0934 

0.0732 

0.0651 

0.0515 

0,0454 

1,  6  75 

1, 975 

1.358 

1.290 

0.873 

0.792 

0.556 

0.522 

0.377 

0.  352 

0.258 

0.240 

0. 1774 
0. 1639 

0. 1223 
0.1123 

0,0845 

0.0771 

0,0585 

0.0529 

0.0405 

0.0364 

Z.  075 

7..  175 

1.231 

1.181 

0.753 

0.719 

0.494 

0.490 

0.332 

0.  314 

0.224 

0.211 

0. 1524 
0.1425 

0. 1038 
0.0964 

0.0707 

0.0627 

0.0483 

0.0442 

0.0329 

0.0300 

.7.  77  5 

Z.  375 

1.137 

1.097 

0.690 

0.  663 

0.449 

0.429 

0.298 

0.284 

0.1992 
0. 1887 

0. 1337 
0.1259 

0.0899 

0.0842 

0.0605 

0.0563 

0.0407 

0.0377 

0.0275 

0.0252 

Z.  475 

Z.  575 

1.061 

1.029 

0.639 

0.618 

0.412 

0.396 

0.271 

0,259 

0. 1792 
0. 1706 

0.1189 
0. 1126 

0.0703 

0.0744 

0.0526 

0.0492 

0.0350 

0,0326 

0.0233 

0.0216 

Z.  67  5 

Z .  7  7  5 

0.999 

0.972 

0.598 

0.579 

0.  382 
0.369 

0.249 

0.239 

0.1628 

0.1556 

0. 1069 
0. 1016 

0.0703 

0.0665 

0.0462 

0.0435 

0.0304 

0.0285 

0,0200 

0.01863 

Z.  S7  5 

Z.975 

0.946 

0.922 

i 

0.562 

0.546 

0.356 

0.345 

0.230 

0.221 

0,1490 
0. 1428 

0.0968 

0.0924 

0.0630 

0.0598 

0,0410 

0,0387 

0,0267 

0.0251 

0,01739 

0.01626 

3.  075 

3,  1  75 

0.899 

0.877 

0.531 

0.516 

0.334 

0.323 

0.213 

0.206 

0. 1371 
0. 1317 

0.0882 
0. 0844 

0.0568 

0.0541 

0.0366 

0.0347 

o!0236 

0. 0223 

0.01523 

0.01429 

3.  Z75 

3.  375 

0.  856 

0.  836 

0.502 

0.489 

0.  314 
0.  304 

0. 1988 
0, 1921 

0. 1266 
0. 1218 

0.0807 

0.0773 

0.0515 

0.0491 

0,0329 

0.0312 

0.0210 

0.01986 

0.01342 

0 .01263 

3.4/5 

3.  5  /5 

0.  816 
0.797 

0.476 

0.464 

0.  295 
0.286 

0. 1856 
0. 1794 

0. 1172 
0. 1128 

0.0741 

0.0710 

0.0469 

0,0448 

0, 0297 
0.0282 

0.01878 

0.01189 

0.01121 

3.67  5 

3.  775 

0.778 

0.760 

0.452 

0.440 

0.278 

0.270 

0. 1734 
0. 1677 

1  0.1086 
0. 1046 

0. 0681 
0,0653 

0.0427 

0.0408 

0.0268 

0.0255 

0.01684 

0.01595 

0.01570 

0.01000 

3.  875 

3.  97  5 

0.741 

0.722 

0.428 

0,416 

0.262 

0.535 

0. 1620 
0. 1565 

0.1007 

0.0969 

0.0624 

0.0600 

0.0390 

0.0372 

0.0243 

0.0231 

0.05112 

0,01432 

0.00941 

0.00888 

4.075 

4 .  1/5 

0.703 

0.684 

0.404 

0.392 

0.245 

0.237 

0. 1510 
0. 1455 

0.0931 

0.0894 

0.0575 

0,0550 

0.0355 

0.0339 

0.0219 

0.0208 

0.01356 

0.01283 

0.00838 

0,00790 

4.  Z7  5 

4.  37  5 

0.664 

0.627 

0.379 

0.  366 

0,229 

0.220 

0. 1399 
0. 1343 

0.0857 

0.0820 

0.0526 

0.0501 

0,0322 

0.0306 

0.01977 

0.01870 

0. 01212 
0.01143 

0.00744 

0.00699 
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TABLE  7 


Values  of  T 


5 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 . 0  <■  5 

7.92 

3.  38 

1.  805 

1.051 

0.642 

0.403 

0.257 

0. 1616 

0. 0563 

1  .  i7f. 

3.82 

1.557 

0.794 

0.443 

0.259 

0. 1561 

0.0959 

0.0595 

0.0361  . 

0.0521 

1  .  27b 

2.  51 

0.  985 

0.484 

0.260 

0. 1464 

0.0850 

0.0503 

0.0303 

0.01826  ; 

0.01075 

1.37  5 

1. 852 

0.705 

0.  335 

0. 1738 

0.0945 

0.0529 

0.0302 

0.01752 

0.01024 

0,00600 

1.475 

1.460 

0.540 

0.249 

0, 1249 

0.0657 

0.0355 

0,01961 

0,01097 

0.00619  ! 

0.00353 

1  .  575 

1.198 

0.431 

0. 1930 

0.0940 

0.0479 

0.0251 

0.01337 

0.00722 

■0.00394  i 

0.00217 

1.675 

1.011 

0.  355 

0. 1546 

0.0731 

0.0361 

0.01832 

0,00946 

0.00494 

0.00260  ; 

0. 001385 

1  .  77  5 

0.872 

0.299 

0. 1268 

0.0583 

0.0280 

0.01375 

0,00687 

0.00347 

0.001767 

0.000913 

1. 87  5 

0.763 

0.256 

0. 1060 

0.0474 

0.0221 

0.01054 

0,00511 

0.00250 

0.001231 

0.000618 

1. 975 

0.677 

0.223 

0,0899 

0.0392 

0.01774 

0.  00822 

0.00386 

0.00183: 

0.0;00875 

0.000429 

>;07  5 

0.606 

0. 1954 

0.0771 

0.0328 

0.01445 

0.00650 

0.00297 

0.00136'i 

0.000633 

0.000303 

2.  175 

0.548 

0. 1732 

0.0669 

0.0277 

0.01190 

0.00521 

0.00231 

0.00103! 

,0.000465^ 

0.000219 

7,  275 

0.498 

0. 1548 

0.0585 

0. 0237 

0.00990 

0.00422 

0.001822 

0.000793 

>0.000347' 

0.000161 

■  2.375 

0.456 

0. 1392 

0.0515 

0.0204 

0,00831 

0.00345 

0.001451 

0.00061^ 

0.000261; 

0, 000120 

2.47.5 

0.419 

0, 1259 

0.0457 

0.01765 

0.00703 

0.00284 

0.001166 

0.00048] 

0.000199i 

0.000092 

2.  575 

0.  387 

0. 1443 

0.0407 

0.01540 

0.00599 

0.00237 

0,00094£ 

0.00038C 

0,000153| 

0.000071 

7.675 

0.  358 

0.1045 

0,0365 

0.01351 

0.00514 

O'.  001981 

0.000772 

0,000303 

0.000119 

0.000056 

2.  775 

0.333 

0.0957 

0.0328' 

0.01191 

0.00443 

0.00167C 

0.00063f 

0.000243 

i  0.000093 

0.000045 

2.675 

0.311 

0.0880 

0.0297 

0.01055 

0.00383 

0.001415 

0.  000526 

1  0.000197 

0.000074 

0.000036 

2.  975 

0.290 

0.0811 

0.0269 

0.00939 

0.00335 

0.001207 

0.000439 

0..00016C 

0.000059 

0.000030 

3.  075 

0.272 

0.0750 

0.  0245 

0.00838 

0.00293 

0.001034 

0..  000366 

0.000131 

O.0OOO47 

0.000026 

3,175 

0.255 

0,0695 

0.0223 

0.00751 

0,00257 

0.000891 

0.00031C 

0.000106 

0.000038 

Q,  000022 

3.  275 

0„239 

0.0645 

0.0204 

0,00676 

0.00227 

0.000771 

0.000263 

0.000090 

0.000031 

6.000019 

3.  37  5 

0.225 

0.0600 

0.01872 

0.00609 

0.00201 

0.000670 

0.00022^ 

0. 000075 

0.000025 

0.000017 

3. 47  5 

0.212 

0.0558 

0.01720 

0.00551 

0.001790 

0.000586 

0. 000192 

0.000063 

0.000021 

0.000015 

3 .  57  5 

0. 1992 

0.0520 

0.01583 

0.00500 

0.001597 

0.000514 

0.000166 

0.000053 

0.000014 

3.  67  5 

0. 1876 

0.0485 

0,01459 

0.00454 

0,001430 

0.000452 

0. 000143 

0.000045 

0.000014 

0.000013 

3.  775 

0. 1767 

0.0453 

0.01346 

0.00414 

0.001284 

0.000400 

0.00012E 

0.  00003S 

0.000012. 

0. 000012 

3.  875 

0. 1665 

0.0423 

0.01244 

0.00378  i 

0 .001156 

0.000355 

0.000109 

0.000033 

0.000010 

0.000012 

3. 975 

0.  1567 

0.0395 

0.01150 

0.00345  1 

0.001044 

0.000317 

0.000096^ 

0.00002S 

0.000008 

0.000011 

•1 .  0  /  5 

0. 1474 

0.0369 

0.01064 

0.00316 

0.000945 

0.000283 

0.000083. 

0. 000025 

0.000007 

0.000011 

4.1/5 

0.  1385 

0.0344 

0.00984 

0.00290 

0.000858 

0.000254 

0. 000075 

0, 000022 

0.000006 

0,000010 

4.  27  5 

0. 1298  ' 

0.0321 

0.00911 

0.00266 

0,000780 

0.000229 

0.00006V 

0.000019 

0.000005 

0.000010 

4  .  3  5 

Q.  1214  i 

0.0298 

0.00842 

0.00243 

0.000710 

0.000207 

1 

_ L 

0.000066 

0.000017 

0.000005 

0.000009 

71- 
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TAF5I-K  8 

Values  of  X 

e 


\ 

1 

3 

4 

in 

1  .  N  7  5 

11.22 

8.77 

7,54 

6.74 

6. 15 

5.69 

5.28 

4.81 

i  .  1  . 

6.43 

4.93 

4. 15 

3.64 

3.26 

2.95 

2.70 

2.47 

17225” 

_ 

1 .  :-r: 

4.79 

3.61 

2.98 

2.57 

2.25 

2,01 

1.799 

1.622 

1.464 

1.299 

1.3.  f. 

3.92 

2.90 

2.36 

1.994 

1.720 

1.501 

1,321 

1.168 

1.036 

0.919 

1 ,  i  /  5 

3.36 

2.45 

1.962 

1.631 

1,382 

1.184 

1.022 

0.887 

0.771 

0.673 

1  .  5  /^. 

2.97 

2.13 

1.683 

1.376 

1.146 

0.965 

0.818 

0.696 

0,594 

0.508 

1  .  '  '5 

2.67 

1.896 

1,474 

1. 187 

0.973 

0.805 

0.670 

0.559 

0.468 

0.  393 

i  .  7  /  S 

2.44 

1.710 

1.311 

1.040 

0.839 

0.682 

0.558 

0.458 

0.  376 

0.310 

1  .  o-  F  3 

2,25 

1.559 

1.180 

0.922 

0,732 

0.586 

0.472 

0,  380 

0.  307 

0.248 

1  .  9  /  3 

2.10 

1.434 

1.071 

0.826 

0,  646 

0,509 

0.403 

0.320 

0.254 

0.202 

9/3 

1.961 

1 . 328 

0.980 

0,745 

0.575 

0.446 

. 

0.348 

0.271 

0.212 

0. 1660 

7 .  17  3 

1.844 

1.236 

0.901 

0.677 

0.515 

0,  394 

0.302 

0.232 

0. 1790 

0.1379 

7.  ...  7  3 

1.742 

1.156 

0.  833 

0.618 

0.464 

0.350 

0.265 

0,201 

0. 1522 

0.1155 

7.  37  S 

1.651 

1.085 

0,773 

0.566 

0.420 

0.  313 

0.233 

0. 1742 

0.1303 

0,0975 

7. .  4  .  '3 

1.568 

1.021 

0.720 

0.521 

0.381 

0.280 

0.206 

0. 1522 

0.1123 

0.0829 

7 ,  3  y  3 

1.494 

i 

0.964 

0.673 

0.482 

0.348 

0.252 

0.1836 

0,1336 

0.0973 

0.0709 

OF':. 

1 

1.425 

0.912 

0.630 

0.446 

0.318 

0.228  1 

1 

0.1639 

0,1178 

0.0848 

0.0610 

...  ,  '  F 

1.362  j 

|0. 864 

0.591 

0.414 

0.292 

0.207 

0.1469 

0.1044 

0.0742' 

0,0527 

...  '  /  D 

1.  303 

0,820 

0.556 

0.385 

i 

0.269 

0.1883 

Oo 1321 

0.0927 

0.0651 

0.0458 

...  h  3 

1.248 

0.779 

0.523 

0,359 

0.248 

0.1717 

! 

0,1191 

0,0827 

0.0574 

0.0399 

:> .  /  S 

1.196 

1 

0.741 

0.493 

0.  335 

0.229 

i 

0,1569  j 

i 

0.1076 

0,0739 

0:0508 

0.0349 

3,  1  ' 

1.146 

0.705 

0.465 

1 

0.313 

0.212 

0.1436 

j 

'  0.0975 

0.0662 

0,0450 

0.0306 

.3.7/5 

1.099 

|0.671 

0.439  ' 

0.294 

0.1959 

0,1316 

1  0,0885 

0.0595 

0.0400 

0.0269 

3  .  3  .*  5 

1.054 

0.639 

0.415 

0.274 

0.1816 

0. 1208 

0.0804 

0.0535 

0.0357 

0.0237 

5,47  5 

1.010 

0.  609 

0.392 

0,256 

0.1684 

0,1109 

0,0731 

'0.0482 

0.0318 

0.0210 

3.5/  ■> 

0.  968 

0.579 

0.370 

’ 

0.240 

0. 1562 

0.1019 

0.0666 

0,0435 

i  ..  .  .  .7 

0.01858 

3  .  b  /  5 

0.927 

0.551 

0.  349 

0.224  ! 0.1449 

0.0937 

'  0.0606 

0.0393 

0.0254 

0.01646 

3 .  /  7  3 

0.  886 

0.524 

0.329 

0.210 

0.1343 

0.0861 

0.0552 

0.0354 

0.0228 

0.01460 

5.  K73 

0.846 

0.497 

0.310 

1.960 

0. 1244 

. 

0.0791 

0.0503 

0.0320 

0.0204 

0,01296 

^.9.5  j 

1 

0.  807 

0.471 

0.291 

1.828 

0.1151 

0.0726 

0.0458 

0,0289 

0,01822 

0,01150 

•i .  7  /  5 

0.767 

0.445 

0.273 

1,702 

0. 1063 

0.0665 

0,0416 

0.0260 

0.01628 

0.01019 

4.1:  5 

0.726 

0.419 

0.256 

j 

1.580 

0.0979 

0.0607 

0.0377 

0.0234 

0.01452 

0.00902 

•1  .7/3' 

0.685 

0.  393 

0.238 

1 

1.460 

0.0898 

0.0553 

0,0340 

0.0210 

0,01291 

0.00796 

4.33  1 

. 0.643 

0.  366 

0.220 

1, 343 

0.0820 

0.0501 

_ 1 

0.0306 

0.01870 

0.01143 

0.00699 

.72- 
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rAllLE  9 
Values  of  p 


I  ^ ^ ^ - 


X 

1 

•> 

3 

^  -1 

f - 

') 

r.i 

I  .  0  >0 

0,0101 

0.0207 

.0,0318 

0.0435 

0.0558 

0.0687 

0.0822 

0.0964 

0.  1113 

1  .  130 

0.0309 

0.0664 

0.1073 

0. 1544 

0.209 

0.271 

0.  342: 

0.424  . 

0.519; 

1.  330 

0.0526 

0, 1184 

0.201 

0,303 

0.432 

0.  593 

0.793 

1.044 

1.358 

1  ,  ^50 

0.0752 

0. 1768 

0.314 

0.499 

0.749 

1.087 

1.542 

2.16 

2.99 

1.130 

0.0989 

0.242 

0.450 

0.752 

1.189 

1,  823 

2.74 

4.08 

6.01 

8.  81 

1  .  330 

0. 1235 

0.  315 

0.612 

1.072 

1.786 

2.  89 

4.61 

7.26 

11.38 

17.76 

1 . 6  50 

0. 1494 

0.396 

0.803 

1,474 

2.58 

4.41 

7.42 

12.40 

20.6 

34.2 

1.  7  50 

0.1764 

0.485 

1.026 

1.972 

3.63 

6,52 

11.59 

20.5 

36,0 

63.2 

J  .  S50 

0.205 

0.584 

1.285 

2.58 

4.98 

9.42 

17.63 

32.  8 

60.9 

112.9 

1. 950 

0.235 

.0.692 

1.584 

3.  32 

6.72 

13.3 

26.2 

51,4 

100.4  . 

196.0 

3.050 

0.266 

0.811 

1.928 

4.22 

8.91 

18.5 

38.3 

78,7 

161.6 

332 

3.  150 

0.299 

0.941 

2.32 

5.29 

11.67 

25.4 

54.9 

118.3 

255 

548 

2  50 

0.333 

1.083 

2.77 

6.57 

15.11 

34.3 

77,6 

174.9 

394 

.  887 

2.  350 

0.370 

1.239 

3.28 

8,08 

19.  36 

45.9 

108.2 

255 

599 

1407 

2 . 4  50 

0.408 

1.409 

3.86 

9.87 

24.6 

60.6 

149.0 

365 

896 

2190 

2,55  0 

0.449 

1.595 

4.52 

11.97 

30.  1 

79.4 

203 

518 

1321 

3370 

2.  650 

0.493 

1.798 

5.  26 

14.42 

38.7 

103.1 

274 

726 

1924 

5010 

a.  7  50 

0.538 

2.02 

6.09 

17.29 

48, 1 

132.8 

336 

1006 

2770 

7610 

2.  850 

0.587 

2.26 

7.03 

20.6 

59.4 

169.8 

485 

1381 

3940 

11220 

2.  9  50 

0.639 

2.53 

8.09 

24.5 

72.9 

216 

637 

1880 

5550 

16360 

3. 0  50 

0.695 

2.  81 

9.28 

29.0 

89.  1 

273 

i  832 

2540 

7740 

23600 

j 

3.  150 

0.754 

3.  13 

10.62 

34.2 

108.5 

342  . 

1079 

3400 

10710 

j  33700 

3,  250 

0.818 

3.48 

12. 12 

40,2 

131.5 

428 

1392 

4530 

14700 

i  47800 

3.  350 

0.  887 

3.  86 

13.81 

47.2 

158.8 

533 

1786 

4990 

20100 

I  67200 

3.  ISu 

0.961 

4,28 

15.71 

55.2 

191.3 

661 

2280 

7870 

27200 

j  93700 

3.  550 

1.041 

4.74 

17.  85 

64.  4 

230  ’ 

817 

2900 

10300 

36500 

1  129800 

3.  6  50 

1  1.128 

1 

5.24 

20,  3 

75.1 

275’ 

1006 

3670 

13400 

48900 

178600 

3.'V50 

1.222 

5.82 

23.0 

87.5 

1  329 

1236 

4630 

17380 

65200 

244000 

3.  8  50 

1.326 

6,43 

26.  1 

101.7 

1  393  ’ 

1514 

5830 

22500 

86400 

333000 

3. 950 

1.439 

7.  12 

1  29. 6 

118.3 

'  469 

1852 

7320 

28900 

114200 

451000 

9 . 0  50 

1.564 

7.90 

33.  6 

137.5 

:  558 

2260 

9160 

37100 

150300 

609000 

■1 .  1  50 

1.703 

8.77 

38.  1 

'  159.8 

665 

2760 

11460 

47600 ! 

197400 

8190^0 

4.  25U 

1.857 

9.75 

43.  3 

185.9 

792 

3370 

14310 

60800 

259000 

4.  3  50 

2.03 

10.86 

49.3 

216 

943 

4110 

17860 

77700 

33800 

4 . 4  50 

2.23 

12.  L3 

.56.2 

252 

1125 

5010 

22300 

99200 

441000 

4.  550 

2.25 

13.59 

64.  3 

295 

:  1344 

6120 

27800 

12670C 

1  576000 

4  .  b  5  0 

2.70 

15.28 

73.7 

346 

1610 

7490 

34800 

16190( 

1  753000 

4 .  V  50 

3.00 

17.25 

84.9 

407 

1934 

9190 

43600 

20700( 

)  985000 

4.0  59 

3.  35 

19.59 

98.3 

480 

2330 

11320 

54900 

26600( 

i  ■ 

4.  <‘50 

3.76 

22.4 

114,6 

571 

2830 

14010 

69300 

34300( 

5.  CSO' 

4.26 

25.8 

134.5 

684 

3460 

17460 

88200 

44500( 

5.  1  50 

4.88 

30.0 

159.5 

826 

4260 

21900 

113000 

58200( 

73- 
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TABLK  10 

Values  of  \ 


I 


\  , 

\ 

\ 

I 

} 

5 

/ 

10 

1 .  (' 

7,31 

4.57 

3.08 

2.  12 

1.482 

1.045 

0.749 

0.571 

0.434 

1 .  i  ^>0 

3.68 

2.29 

1.530 

1.049 

0.727 

0.508 

0.355  : 

0.244 

1 . ::  ;-.o 

•2.72 

1.685 

1*  121 

0.763 

0.526 

0.  365 

0.254 

0. 1766 

0.1133 

2.26 

1.  389 

0.918 

0.  622 

0.426 

0.294 

0.203 

0.  1411 

0.0974 

1 . 

1.974 

1.208 

0.794 

0.535 

0.  364 

0.249 

0. 1714 

0. 1182 

0.0814 

0.0553 

i .  S...0 

1.781 

1.084 

0.709 

0.474 

0.  321 

0.  218 

0. 1491 

0.  1020 

0.0700 

0,0480 

1  .  h  1 0 

1.638 

0.992 

0.645 

0.429 

0.289 

0. 1950 

0.1322 

0.0899 

0.0612 

0.0417 

\  0 

1.528 

0.921 

0.596 

0.394 

0,263 

0. 1767 

0.1190 

0.0803 

0.0543 

0.0367 

1  .  -S^O 

1.439 

0.  864 

0;556 

0.  366 

0.243 

0. 1619 

0.1083 

0.0726 

0.0487 

0.0327 

i  . 

1.367 

0.817 

0.523 

0.342 

0.226 

0. 1496 

0.0994 

0.0661 

0,0440 

0.0294 

oso 

1. 305 

0.777 

0.495 

0.322 

0,211 

0.1390 

0.0918 

0.0607 

0.0401 

0.0266 

Z.  150 

1.252 

0.742 

0.471 

0.  304 

0. 1985 

0. 1300 

0.0852 

0.0560 

0.0368 

0.0242 

;:50 

1.206 

0.712 

0.449 

0.289 

0.1875 

0. 1220 

0.0795 

0.0519 

0.0339 

0.0221 

3  50 

1.165 

0.685 

0.430 

0.276 

0, 1777 

0.  1150 

0.0745 

0,0483 

0.0314 

0.0204 

z.  -ibo 

1. 129 

0.661 

0.414 

0.263 

0. 1690 

0. 1087 

0.0700 

0.0452 

0.0291 

0.01880 

Z.  550 

1.096 

0.640 

0.  398 

0.252 

0. 1611 

0.  1030 

0.0660 

0.0423 

0. 0272 

0.01742 

Z .  6  50 

1.066 

0.620 

0.384 

0.242 

0. 1539 

0.0979 

0.0624 

0.0398 

0.0254 

0.01619 

Z.  V50 

1.038 

0.  602 

0.  372 

0.233 

0. 1473 

0.0933 

0.0591 

0.0375 

0.0238 

0.01510 

Z.  850 

1.013 

0.  585 

0.  360 

0.225 

0. 1413 

0.0890 

0.0561 

0.0354 

0.0223 

0.01411 

a.  950 

0.989 

0.569 

0.  349 

0.217 

0. 1357 

0.0851 

0.0534 

0.0335 

0.0210 

0.01321 

3.  050 

0.967 

0.555 

0.339 

0,210 

0, 1305 

0.0814 

0.0508 

0. 0318 

0.01985 

0.01240 

3.  150 

0.946 

0.541 

0.  329 

0.203 

0.1257 

:  0.0781 

0.0485 

0.0302 

0.01875 

0.01166 

3.  3  50 

0.926 

0.528 

0.320 

0.1965 

^  0.1212 

0.0749 

'  0.0463 

0.0287 

0.  01.775 

j  0.01099 

3,  3  50 

0.908  ■ 

0.516 

0.311 

0. 1904 

0. 1170 

0.0720 

0.0443 

0.0273 

0.01683 

1  0.01037 

3.-:50 

0.890 

0.504 

0.  303 

0. 1847 

0.1130 

0.0692 

0.0424 

0,0260 

0.01597 

0,00980 

■  3.550 

0.  873 

0.493 

0.295 

0. 1793 

0. 1092 

0, 0666 

0.0407 

0.0248 

0,01518 

0.00927 

3.  6  50 

0,856 

0.483 

0.288 

0. 1741 

0. 1057 

0.0642 

0.0390 

0,0237 

0.01444 

0.00878 

3,  750 

0.  840 

0.472 

0.281 

0. 1692 

0. 1023 

0.0619 

0.0375 

0.0227 

0. 01375 

0.00833  , 

3.  :::50 

0.825 

0.462 

0,274 

0, 1645 

0.0990 

0.0597 

0.0360 

0,0217 

0.01310  ; 

0.00791 

3. 950 

0.  810 

0.453 

0.268 

0. 1600 

0, 0959 

0.0576 

0.0346 

0.0208 

0.01250 

0.00751 

4.0  50 

0.795 

0.443 

0.261 

0, 1556 

0.0930 

0.0556 

0.0333 

0.01993 

0. 01193 

0.00714 

-  .  1  50 

0.781 

0.434 

0,255 

0. 1514 

0.0901 

0.0537 

0.0320 

0.01910 

0.01139 

0.00680 

•1 .  .150 

0.766 

0.425 

0.249 

0. 1473 

0.0874 

0.0519 

0.0308 

0.01832 

0.01089 

0.00647 

4.  550 

0.752 

0.416 

0.243 

0, 1433 

0.0847 

0.0501 

0.0297 

0.01757 

0.01041 

0.00616 

4 . 4  50 

0.738 

0.408 

0.237 

0. 1394 

0.0821 

0.0484 

0.0286 

0,01686 

0. 01000 

0.00587 

.  550 

0.724 

0.  399 

0.231 

0. 1356 

0.0796 

0.0468 

0.0275 

0.01618 

0.00952 

0.00560 

4 . 6  50 

0.710 

0.390 

0,226 

0. 1318 

0.0772 

0.0452 

0.0265 

0.01553 

0.00910 

0.00533 

4.  7  50 

0,696 

0.  382 

0,  220 

0. 1281 

0.0748 

0.0437 

0.0255 

0.01490 

0. 00870 

0.00508 

4.5  5Q 

0.682 

0.  373 

0.214 

0. 1244 

0.0724 

0.0421 

0.0245 

0.01429 

0,00832 

0.00484 

.  OZT; 

0.667 

0.  364 

0.209 

0, 1208 

0.0700 

0.0406 

0,0236 

0.01370 

0. 00795 

0.00461 

5.  0  50 

0.652 

0.  355 

0.203 

0. 1171 

0.0677 

0.0392 

0.0227 

0.01311 

0.00759 

0.00439 

5.  1  50 

0.636 

0.345 

0. 1970 

0.1133 

0.0653 

0.0377 

0.0217 

0.01254 

0. 00724 

0,00417 

-74- 
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i  A  ;  M  .  ;  ; 

Values  of  T 


X 

1 

! 

■ 

1 

j 

1  ■ 

1  . 

12.  81 

5.  39 

2.  829 

1.616 

0.968 

0.599 

0.384 

0.271 

0. 1769 

!  .  .1  A'A 

5.14 

2.06 

1.029 

0.561 

0.321 

0. 1886 

0.1128 

0.0661 

0.0944 

1 . vA 

3.24 

1.247 

0.599 

0.  314 

0. 1726 

0.0978 

0.0564 

0.0327 

0.01594 

1  . 

2.  35 

0.  857 

0.406 

0.205 

0. 1090 

0.0596 

0.0332 

0.01874 

0.01054 

]  .  1 

1.840 

0.664 

0.298 

0. 1460 

0.0749 

0.0396 

0.0213 

0.01161 

0.00637 

0.00337 

1 .  ^:  ao 

1.504 

0.528 

0.230 

0. 1092 

0.0543 

0.0277 

0.01443 

0.00760 

0.00406 

0.00216 

1  .  r  30 

1.267 

0.433 

0. 1836 

0.0846 

0.0408 

0.0202 

0.01016 

0.00517 

0.00267 

.0.001378 

i  .  ''  A  0 

1.092 

0.  364 

0. 1504 

0.0673 

0.0315 

0.01511 

0.00737 

0.00363 

0.001819 

0.000907 

1 .  >’o0 

0.957 

0.312 

0. 1256 

0.0547 

0.  0249 

0.01157 

0.00547 

0.00261 

0.001271 

JD.  000613 

1 .  ')  "0 

0.850 

0,271 

0,  1065 

0.0452 

0.0200 

0.00902 

0.00414 

0.001916 

0.000908 

0.  000424 

•i .  0  30 

0,763 

0.238 

0.0914 

0.0378 

0.01626 

0.00715 

0.00318 

0.001430 

0.000662 

0.000299 

130 

0.691 

0.212 

0.0794 

0.0320 

0.01340 

0.00573 

0.00248 

0.001083 

0.000491 

0.  000216: 

0.630 

0. 1895 

0,0695 

0.0274 

0.01117 

0.00465 

0, 001958 

0.000832 

0.000371 

0.000158 

Z.  330 

0.578 

0. 1708 

0.0613 

0.0236 

0,00939 

0.00381 

0.001562 

0.000646 

0.000284 

0,000118 

0.  ■•130 

0.534 

0.1549 

0.0545 

0.0205 

0.00795 

0.00315 

0.001258 

0.000507 

0.000221 

0.000089 

.  3  3  0 

0.495 

0. 1412 

0.0487 

0.01790 

0.00679 

0.00262 

0.001022 

0.000401 

0.000174 

0.  000069 

3.  0  3f- 

0.460 

0. 1293 

0.0437 

0.01573 

0.00583 

0.00220 

Q. 000836 

0.000320 

0.000139 

0.000053 

Z .  V  3  J 

0.430 

0. 1189 

0.0945 

0.01390 

0.00504 

0.001855 

0.000690 

0,000258 

0.000112 

0.000042 

2  .  .  3  0 

0.403 

0. 1097 

0.0357 

0.01234 

0.00438 

0.001575 

0.000572 

0.000209 

0.000092 

o! 000034 

0.378 

0. 1016 

0.0325 

0.01100 

0.00382 

0.001345 

0.000478 

0.000171 

0.000077. 

0.000028 

^ ,  0  *.o 

0.356  i 

0.0943 

0.0297 

lo. 00985 

0.00335 

0.001154 

!  0,000401 

t  , 

0.000140 

0.000065 

0.000023 

-j .  130 

'  0.336 

0.0877 

0.0272 

0.00885 

0.00295 

0.000995 

1  ' 

0.000339 

0.000116 

0.000055 

1  0,000019’ 

3 .  30 

0.317 

0.0819 

0.0249 

0.00797 

:  0.00261 

0.000862 

I  0.000287: 

0.000096 

0.000048 

j.  0.000017 

t  y..  -  ' 

3.  3  30 

0.300 

0.0766 

0.0229 

0.00721 

1 

0.00231 

0.000749!  0.000245 

j 

0.000080 

0.000042 

i  0,000015 

3.‘150. 

0.285 

0.0717 

0,0212 

1 

0.00654  ! 

0.00206 

1 

0.00065410.000210 

0.000067 

i 

10.000037 

0,000013 

3.  :-''30 

0.270 

0.0673 

0.0196 

0.00594  j 

0.001836 

0.00057310.000180 

0.000057 

:  0.000033 

0,000011 

3 . 6  30 

0.257 

0.0632 

0.0181 

0.00542  j 

0.001645 

0.000504 

0.000155 

0.000048 

:  0.000030 

0.000010 

3  .  7  3  0 

0.244 

0.0595 

0.0168 

0.00495 

0,001476 

0.000446 

0.000135 

0,000041 

0.000028 

0.000010 

3  .  o  5  0 

0.235 

0. 0560 

0,0156 

0,00453 

0.001331 

0.000394! 

0.000117 

0.000035 

0.000026 

0.000009 

0.221 

0.0528 

0.0146 

0.00416 

0.001202 

0.000350 

0.000102 

0.000030 

0.000024 

0.000008 

4  .  (  3  0 

0.211 

0.0499 

0.0136 

0.00382 

0.001086 

0.000312 

0.000090 

0i000026 

0.000023 

0, 000008 

J .  1  -0 

0.201 

0.0471 

0.0127 

0.00352 

0.000986 

0.000279;  0.000079 

0.000022 

0.000022 

0, 000008 

4.  OrJ 

0.1919 

0.0445 

0,0118 

0.00325 

0.00089^ 

0.000250'  0,000070 

1 

0.000019 

0.000021 

0. 000007 

^30 

0.1830 

0.0421 

0..0111 

0.00300 

0.00082C 

0,  0002251  0.  0000621 

0.000017 

0.000020 

0.  000006' 

•; .  4  30 

0.1745 

0.0398 

0.0104 

0.00278  ' 

0.00074S 

0.000203 

0.000055 

0.000015 

0,000019 

0,000006 

..3  30 

0.1665 

0.0377 

0.00972 

0.00257 

0.00068^ 

0.000183 

0.000049 

0.000013 

0.000019 

0.000006 

•'i  .  3  0 

0.1587 

0.0357 

0.00911 

0.00239 

0.000626 

0.000166 

0.000044 

0.000011 

0.000018 

0.000005 

4 .  3(-' 

0.1512 

0.0337 

0.00854 

0.00222 

0.000576 

0.000151 

0.000039 

0.000010 

0.000018 

0.000005 

■1  .  :  ■ 0 

0, 1440 

0.0319 

0.00802 

0.00206 

0.000532 

0.000137 

0.000035 

0.000009 

0.000017 

0. 000005 

4.  0  3'./ 

0.1370 

0.0302 

0.00752 

0.001917 

0.00049] 

0.00012( 

>0.000032 

0.000008 

0.000017 

0.000004 

5 .  ' '  3  • 

0.1302 

0.0285 

0.00706 

0.001785 

0.00045: 

O.OOOllf 

0. 000029 

0.000007 

0.000016 

0.000004 

1  M 

0.  1235 

0.0269 

0.00662 

0.001662 

0. 000416 

0.000105 

0.000026 

0.000006 

0.000016 

0.000004 
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Values  of  X 
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\ 

\ 

A 

1 

> 

- - 

i - 

]  S> 

1 . 

i  V) 

15.  82 
7.69 

12.42 

5.92 

10.74 

5.01 

9. 65 
4.41 

8.86 

3.97 

8,27 

3.61 

7.84 

3.31 

7.74 

2.96 

■ 

5.52 

4.17 

3,47 

3.00 

2.64 

2.  36 

2.  13 

1.909 

1.563 

1 

i . 

V'.O 

4.45 

3.  31 

2.70 

2.29 

1.985 

1.740 

1.538 

1.364 

1.204 

I . 

Ail 

3.79 

2.77 

2,23 

1.858 

1.581 

1.361 

1, 180 

1.028 

0,896 

0.769 

1  . 

3'>r; 

3.33 

2.40 

1.902 

1.561 

1.306 

1.104 

0.940 

0.803 

0.  688  ■! 

1 

0.591 

1 . 

..i: 

3.00 

2.13 

1.663 

1.343 

1. 105 

0.918 

0.767 

0.643 

0,  541 

0.456 

] . 

7  ::,  ^ 

2,74 

1.922 

1.477 

1.176 

0.951 

0.777 

0.638 

0.526 

0.434 

0.359 

)  . 

K-  3 '  ■ 

2.53 

1.752 

1.329 

1.042 

0.831 

0.667 

0.539 

0,436 

0.354 

0.288 

1 . 

'AO 

2.35 

1.613 

1.208 

0.934 

0.733 

0.580 

0.461 

0.367 

0.293 

0.234 

0  50 

2.21 

1.496 

1. 106 

0.844 

0.653 

0.508 

0.398 

0.  312 

0.245 

0. 1920 

130 

2.08 

1.395 

1.019 

0.767 

0.585 

0.449 

0.346 

0.267 

0.206 

0. 1596 

> 

1.969 

1.  307 

. 

0.944 

0.702 

0.528 

0.400 

0.303 

0.231 

0. 1757 

0. 1338 

AO 

1.872 

1.203 

0.879 

0.645 

0.479 

0.358 

0.268 

0.201 

0. 1507 

0.1132 

A  V 

1.784 

1.  162 

0.  821 

0.595 

0.436 

0,  322 

0.238 

.  '  ■ 

0. 1758 

0. 1301 

0.0964 

0  7'} 

1.705 

1.100 

0.769 

0.551 

0,399 

0.291 

0.212 

0. 1547 

0. 1130 

0.0826 

■j  SO 

1.634 

1,045 

0.723 

0.512 

0,367 

0.264 

0.1899 

0.1369 

0.0987 

0.0713 

<  'jO 

1.568 

0.994 

0.681 

0.477 

0.338 

0.240 

0.1708 

0. 1216 

0.0867 

0.0618 

1.507 

0.948 

0.643 

0. 446 

0.312 

0.219 

0. 1541 

0.1085 

0, 0764 

0.0538 

. 

‘AO 

1.451 

0.905 

0.608 

0.418 

0.289 

0.201 

0. 1396 

0.0972 

0.0677 

0.0471 

oso 

1.398 

0,  866 

0.576 

1  0.392 

0.268 

. 

0.1843 

0.1268 

0.0873 

0.0601 

0.0414 

'} . 

ISO 

1.349 

0.829 

0.547 

0.368 

\ 

0,250 

0.1697 

0.1155 

0.0787 

0.0536 

0,0365 

^S. 

oso 

1.302 

0.794 

:  0.520 

0.346 

0.233 

0. 1565 

0. 1054 

0.0711 

0. 0479 

0.0323 

[)■.} 

1.258 

0,762 

0.494 

0.327 

0.217 

0,1446 

0.0965 

0.0644 

0.0430 

0,0287 

3 . 

‘i  S'v 

1.216 

0.732 

0.471 

0.  308 

0.203 

0. 1339 

0.0885 

0,0585 

0.0387 

0.0256 

i . 

5  so 

1.177 

0.703 

0.448 

0.291 

0.1899 

0. 1241 

i  0.0812 

0.0532 

0.0348 

0.0228 

S. 

'.SO 

1.138 

0.675 

0.427 

0.275 

0. 1779 

0. 1153 

i  0.0747 

0.0485 

0,0315 

0.0204 

j . 

7  S  ‘J- 

1.101 

0.649 

0.408 

0.260 

0. 1668 

0. 1071 

0.0689 

0.0443 

0.0285 

0.01832 

‘j . 

o.so 

1. 066 

0.624 

0.389 

0.246 

0.1566 

0.0997 

0.0635 

0.0405 

0.  0258 

0.01646 

■/SO 

1.031 

0.600 

0.  371 

0.233 

0.1470 

0.0928 

0.0587 

i  0.0371 

0.0234 

0.01482 

S'  S  0 

0.998 

0.577 

0,355 

0.221 

0.1382 

0.0865 

0.0542 

0.0340 

0.0213 

0.01336 

ISO 

0.965 

0.555 

0.339 

0.209 

0. 1299 

1  0.0807 

0.0501 

0.0312 

I  0.01939 

0.01206 

,1  . 

0”.0 

0.933 

0.534 

0.  323 

0. 1983 

0.1221 

0.0753 

0.0464 

0.0286 

1  0.01767 

0.01090 

3  50 

0.902 

0.513 

0.  308 

0.1878 

0.1148 

0.0702 

0,0430 

0.0263 

0.01611 

0.00986 

. 

0  •') 

0.871 

0.492 

0.294 

0.1779 

0. 1079 

0.0655 

0.0398 

0.0242 

0.01469 

0.00893 

S;-  0 

0.  840 

0.472 

0.280 

0. 1863 

0.1014 

0.0611 

0.0368 

0.0222 

0.01340 

0.00809 

O'  SO 

0.809 

0.452 

0.267 

0. 1592 

0.0952 

0.0570 

0.0341 

0.0204 

0.01223 

0,00732 

•1  . 

.  SO 

0.779 

0.433 

0.254 

0. 1504 

0.0893 

0.0531 

0,0315 

0.01876 

0.01115 

0.00663 

'  J 

0.748 

0.414 

0.241 

0.1418 

0.0837 

0.0494 

0.0292 

0.01721 

0.01017 

0.00600 

';S0 

0.717 

0.394 

0.228 

0. 1335 

0.0782 

0.0459 

0.0269 

0.01578 

0. 00925 

0.00543 

r 

IS-;. 

.0.685 

0.375 

0.216 

0.1254 

0.0730 

0.0425 

0.0248 

0.01443 

0.00841 

0.00490 

- 

1  so 

0.653 

0.355 

0.203 

0.1174 

0.0679 

0.0393 

0.0228 

1  0.01317 

0.00762 

0.00441 
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Values  of  p 

iJi  -  ^  TABLE  13 

- 1 — ^ ^ - 1 - 1 - T - 1 - T - ^ - 1 - - - 


. . . . . 

1 

'  3 

4 

5 

6 

1 

A. 

9 

10 

1.  IZS 

0.0212 

0.0452 

0.0721 

0.1024 

0.1365 

0. 1748 

0.218 

0.267 

mm 

0.383 

1.  ZZ5 

0.0389 

0.0866 

0.1451 

0.217 

0.305 

0.412 

0.544 

msm 

0.902 

1.145 

1.  3Z5 

0;0572 

0.1331 

0.234 

0.367 

0,543 

0.777 

1.087 

1.498 

2.04 

2.76 

1,4Z5 

0.0762 

0.1848 

0.340 

0.560 

0.875 

1.323 

1.961 

2.87 

4.17 

6.01 

1 . 

0.0958 

0.242 

0.465 

1.324 

2,112 

3.  32 

5.16 

7.97 

12.24 

1.625 

0.1162 

0.305 

1.111 

1.922 

3.24 

5.38 

8.86 

14.51 

23.7 

1.  :25 

0. 1374 

0.375 

0.783 

1.489 

2.71 

4,81 

8.43 

14.67 

25.5 

44.0 

I  .  r.zs 

0. 1594 

0.450 

0.981 

1.950 

3.72 

6.95 

12.84 

23.6 

43.2 

79.0 

1 . 92^' 

0.  1822 

0.533 

1.209 

2.51 

5.01 

9.83 

19.10 

37.0 

71.3 

137.5 

Z.  0  23 

0.206 

0.623 

1.468 

3.18 

6.64 

13.66 

27.9 

56,6 

114.9  • 

233 

\  Z  \. 

0.231 

0,721 

1.763 

3.98 

8.68 

18.68 

39.9 

85. 1 

181.0 

385 

■  ?  3 

0.257 

0.827 

2.10 

4.92 

11.21 

25.2 

56.3 

125.6 

280 

623 

.Z.  3  2!. 

0.283 

0.942 

2.47 

6.04 

14.32 

33.6 

78.3 

182.4 

424 

987 

2.2Z.9 

0,  311 

1.067 

7.  34 

18.11 

44,2 

107.6 

261 

634 

1537' 

r:.  625 

0.341 

1.201 

3.37 

8,86. 

22.7 

57.7 

146.0 

369 

932 

2350 

2.  625 

0.371 

1.346 

3.91 

10.62 

28.3 

74.6 

196. 1 

515 

1352 

3550 

2.  i'Z5 

0.404 

1.503 

4.50 

12.66 

34.9 

95,5 

261 

711 

1938 

5280 

Z.  825 

0.437 

1.672 

5.16 

15.02 

42.9 

121.5 

344 

971 

2740 

7750 

2.  925 

0.472 

1.854 

5.90 

17.72 

52.3 

153.4 

449 

1315 

3850 

11250 

3. 025 

0.509 

2.05 

6.71 

20.8 

63.5 

192.5 

583 

1764 

5340 

16140 

3.  125 

0.548 

2.26 

7.62 

24.4 

76.7 

240 

751 

2350 

7330 

22900 

3. 225 

0.589 

2.49 

8.62 

28.4 

92.2 

298 

961 

3010 

10000 

32200 

3.  325 

0.633 

2.74 

9.73 

110.3 

367 

1222 

4060 

13520 

44900 

3. 8  25 

0.678 

3.00 

10.96 

38.2 

1M.6 

451 

,  1546 

5300 

16140 

62100 

3.  525 

0.727 

3.29 

12.32 

44.1 

156.3 

552 

1  1946 

6860 

24200 

85200 

3.  625 

0,778 

3.60 

13.82 

50.9 

185.2 

672 

2440 

8830 

32000 

116100 

3.  V25 

0.832 

3.93 

15.48 

58.5 

21^ 

815 

3040  1 

11320 

42200 

157100 

3.  825 

0.  890 

4.29 

17.  31 

67,1 

258 

986 

3770 

14430 

55200 

211100 

3 .  9  2  5 

0.951 

4.69 

19.34 

76.9 

303 

1189 

4670 

18320 

71900 

282200 

4 , 0  2-5 

1.017  i 

5. 11 

21.6 

87.9 

355 

1429 

5750 

23200 

93200 

375000 

4 .  125 

1.087 

5.57 

24. 1 

100.4 

415 

1713 

7070 

29200 

120300 

496000 

4  .  225 

1. 162 

6.07 

26.  8 

114.5 

485 

2050 

8660 

36600 

155000 

653000 

4.  325 

1,243 

6.  62 

29.9 

130.4 

565 

2450 

10580 

45800 

198000 

856000 

4.4  25 

1.330 

7.22 

33.  3 

148.5 

659 

2920 

12900 

57100 

253000 

4.  525 

1.424 

7.87 

37.0 

169.0 

766 

3470 

15690 

71000 

321000 

4.  o25 

1.526 

8.59 

41.2 

192.2 

891 

4120 

19060 

88200 

408000 

4  .  :  2  5 

1. 637 

9.  37 

45.9 

219 

4890 

23100 

109200 

516000 

4  .  -3  2  5 

1.759 

10.24 

51.2 

249 

wmm 

5800 

28000 

135000 

652000 

4.  925 

1,892 

11.21 

57.1 

283 

1396 

6880 

33900 

167000 

822000 

5.^;>2  5 

2.04 

12.28 

63.7 

322 

1622 

8150 

.  41000 

206000 

5.  125 

2.20 

13.48 

71.3 

367 

1885 

9660 

49500 

254000 

:  .  225 

2.  38 

14.  82 

79.8 

419 

2190 

11460 

59900 

31300C 

5.  525 

2.58 

16.33 

89.6 

479 

2560 

13610 

72500 

38600C 

5.425 

2.81 

18.05 

549 

2980 

16180 

87800 

47600C 

) 

5.  525 

3.07 

20.0 

113.7 

631 

3490 

19280 

106500 

58900C 

* 

5.  625 

3.36 

22.3 

128.7 

727 

■■■ 

23000 

129600 

72900C 

»  ■ 

5.  7  25 

3,71 

24.9 

146.4 

842 

4820 

27600 

158100 

90500C 

) 

-1.  .s25 

4.11 

28.0 

167.4 

979 

5710 

33200 

193700 

i 

1 

5. 925 

4.58 

31.7 

192.6 

1145 

6790 

40200 

238500 

0.  U25 

5. 15 

36.2 

223 

1350 

8140 

49100 

295600 
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Values  of  X 


\ 

1  ’ 

■  3 

i  i 

i 

■ 

u 

1 

i  .  IZb 

4.23 

2,58 

1.682 

1.124 

0.759 

0.516 

0.  346 

0.0911 

I . 

2.99 

1.806 

1.  172 

0.778 

0.523 

0.353 

0.239 

0. 1571 

0,0692 

1 . 

2-.  43  .  . 

1.460 

0.942 

b.  622 

0.415 

0.279 

0. 1881 

0.1262 

0.0574 

1 .  -1 ... 

2.10 

1,257 

0.806 

0.529 

0.351 

0.234 

0. 1570 

0.1052 

0.0686 

'0.0488 

L  ".,!n 

1.884 

1.121 

0.715 

0.467 

0.308 

0.204 

0.1358 

0.0906 

0.0633 

0.0374 

1  .  tJ.b 

1.727 

1.023 

0.649 

0.421 

0.276 

0. 1816 

0,1200 

0.0795 

0,0527 

0.0347 

1.  7.?. 3 

1.607 

0.947 

0.598 

0.385 

0.251 

0. 1642 

0. 1078 

0.0709 

0.0467 

0,0307 

1 .  e.7:'- 

1.512 

0,887 

0.557 

0,  357 

0.231 

0.1501 

0.0979 

0.0639 

0.0418 

0.0274 

!  . 

1.434 

0.  837 

0.523 

0,  333 

0.214 

0.1385 

0.0897 

0,0582 

0.0378 

0.0246 

7 .  v;  7 :  ‘ 

1.369 

0.796 

0.495 

0.  313 

0.200 

0. 1286 

0.0828 

0.0533 

0.0344 

0,0222 

7.  177 

1. 313 

0.760 

0.470 

0.296 

0. 1884 

0.1202 

0.0768 

0.0492 

0.0315 

0.0202 

’  •'  '  c; 

1.265 

0.729 

0..449 

0.281 

0.1779 

0,1128 

0.0717 

0.0456 

0.0290 

0.01849 

7.  .^77 

1.222 

0.702 

0.430 

0.268 

0.1686 

0.1063 

0.0671 

0.0424 

0.0269 

0.01700 

7 , 4  7  b 

1. 184 

0.678 

0.413 

0.256 

0. 1603 

0. 1005 

0.0631 

0.0397 

0.0249 

0.01569 

7.  52 S 

1.150 

0.656 

0.398 

0.246 

0.1529 

* 

0.0953 

0.0595 

0.0372 

0.0232 

0,01454 

1.120 

0.636 

0.384 

0.236 

0. 1461 

0.0906 

0.0563 

0.0350 

0,0217 

0.01352 

2.  /25 

1.092 

0.618 

0.  372 

0.227 

0. 1400 

0.0863 

0.0533 

0.0330 

0. 0204 

0.01260 

7.  >:25 

1.066 

0.601 

0.  360 

0.219 

0. 1343 

0,0824 

0.0507 

0.0311 

0.01915 

0.01178 

2  .  2  3 

1.042 

0.586 

0.3.50 

0.212 

0.1291 

0.0789 

0.0482 

0.0295 

0.01804 

0.01104 

A  02v> 

1,020 

0.571 

0.340 

0.205 

0.1243 

0,0756 

0.0460 

0.0280 

0.01703 

0.01037 

'i.  175 

1.000 

0.558 

0.  330 

0. 1985 

0.1199 

0.0725 

0.0439 

0.0266 

0.01611 

0.00976 

3, 2.25 

0.980 

0.545 

0,  322 

0. 1925 

0. 1157 

0.0697 

0.0420 

0.0253 

0.01527 

0.00921 

-3.  .:-;5 

0.961 

0.533 

0.  313 

0. 1868 

0,1118 

0.0670 

0.0402 

0.0241 

0.01449 

0.00870 

7  .  *1 2  5 

0.944 

i  0,522 

0.306 

0.1815 

0, 1081 

0.0645 

0.0385 

0.0230 

0.01377. 

0.00823 

j».  52  5 

0.927 

0.511 

0.298 

0.1765 

0.1047 

0.0622 

0.0370 

0.0220 

0.01310 

0.00780 

.625 

0.911 

I  0.501 

0.291 

0.1716 

0.  1014 

0.0600 

1  0.0356 

0.0210 

0.01248 

0.00740 

3.7Z5 

0.896 

0.492 

0.285 

0.1671 

0.0984 

0.0580 

!  0.0342 

0.0202 

0.01191 

0.00703 

3.  625 

!  0.882 

0.482 

0.279 

1  0.1628 

0.0955 

0.0560 

1  0.0329 

0.01935 

0.01137  : 

0.00668 

3. 923 

1  0.868 

0.473 

0.272 

0.1586 

0.0927 

0.0542 

1  0.0317 

10.01857 

0.01087 

0,00636 

'1 . 0  2  5 

0.854 

0.465 

0.267 

0.1547 

! 

i  0.0900 

0.0525 

1  0.0306 

0.01783 

i 

0.01040  ! 

0.00606 

4.  125 

0.841 

0.456 

0.261 

i  0. 1509 

0.0875 

0,0508 

I  0.0295 

0,01714 

0.00996  : 

0.00579 

4.  22  5 

0.828 

0.448 

0.256 

1  0.1473 

0.0851 

0.0492 

0.0285 

0.01649 

0.00954  1 

0,00553 

4.  323 

0.816 

0.440 

0.250 

0,1438 

0.0828 

0,0477 

;  0.0275 

0.01587 

0.-00915 

0.00528 

4  .  2  5 

0.804 

0.433 

0.245 

'  0.1404 

0.0806 

0.0463 

1  0.0266 

0.01528 

0.00878 

0.00505 

i .  5  2  5 

0,792 

0.425 

0.240 

0.1371 

0.0784 

0,0449 

0.0257 

1 

0.01472 

0.00843 

0,00483 

-1  .  2  5 

0.780 

0.418 

0.236 

0, 1340 

0.0764 

0,0436 

j 

0.0249 

0,01419 

0.00810 

0.00463 

4  .  /  i-1 5 

0.769 

0.411 

0.231 

0. 1309 

0.0744 

0.0423 

0.0241 

0,01369 

0.00779 

0.00443 

■K  .-.••2  5 

0.757 

0.404 

0.226 

0. 1279 

0.0724 

0.0411 

0.0233 

0.01320 

0.00740 

0.00425 

4.4  25 

0.746 

0.397 

0.222 

0. 1250 

0.0706 

0.0399 

0.0225 

0.01274 

0.00720  i 

0.00407 

5.  .:?  3 

0.735 

0.  390 

0.217 

0. 1221 

0.0687 

0.0387 

0.0218 

0,01229 

0.00693 

0.00391 

123 

0.724 

0,383 

0.213 

0.1193 

0.0670 

0.0376 

0.0211 

0.01187 

0.00667 

0,00375 

.  7.25 

0.712 

0.  377 

0.209 

0. 1166 

0.0652 

0.0365 

0.0205 

0,01146 

0.00642 

0. 00360 

.  i.'3 

0.701 

0,  370 

0.204 

0.1139 

0.0635 

0.0355 

0.01981 

0,01106 

0.00618 

0.0034S 

■- .  -125 

0.690 

0.  363 

0.200 

0.1112 

0.0619 

0.0344 

0.01918 

iO. 01068 

0.00595 

0,00331 

525 

0.678 

0.356 

0.1959 

0. 1085 

0,0602 

0.0334 

0.01856 

10.01031 

0.00572 

0.00318 

5. 625 

0.667 

0.  350 

0. 1917 

0.1059 

0.0586 

0,0324 

0.01796 

0.00995 

0.00551 

0.00305 

'  ■■  4 

0.655 

0.343 

0. 1874 

0. 1033 

0.0570 

0.0315 

0.01737 

0.00959 

0.00530 

0.00293 

3 .  ^•;  3 

0.643 

0.336 

0.1831 

0.1006 

0.0554 

0,0305  i 

0.01679 

0.00925 

0.00500  ' 

0.00281 

5. 

0.630 

0.  328 

0.1787 

0.0980 

0.0538 

0.0295 

0.01622 

0.00891 

0,00489 

0.00269 

6. 02  3 

0.617 

0.321  ' 

0. 1742 

0.0953 

0.0522 

0.0286 

0.01565 

0.00857 

0.00470 

0.00257 
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Values  of  T 

n 

i.  ''  Nv 

\  \ 

1 

3 

4 

‘5 

6 

i:.' 

i— 

1.125 

6*.86 

2.710 

1.332 

0.714 

0.401 

0.231 

0.1329 

0.1349 

0.0332 

0.0206 

1. 2^5 

4.08 

1.546 

0.729 

0.  375 

0.203 

0. 1126 

0,0635 

0,0340 

0.01315  • 

0. 01034 

1.  325 

2.89 

1.059 

0.482 

0.239 

0. 1246 

0.0668 

0,0365 

0.01989 

0.0278 

0.00620 

1.425 

2.23 

0.792 

0,  349 

0. 1676 

0.0843 

0,0437 

0,0230 

0. 01225 

0.00609 

0; 00398 

1  .  525 

1.817 

0.  625 

0.267 

0. 1243 

0.0605 

0.0303 

0.01548 

0. 00800 

0.00414 

0. 001747 

1.625 

1.525 

0,511 

0.212 

0.0958 

0.0452 

0,0220 

0.01084 

0. 00542 

0.00272' 

0.001350 

1  .  725 

1.311 

0.428 

0.  1731 

0.0760 

0.0348 

0.01638 

0,00783 

0. 00379 

0.001848 

0.000904 

1 .  H2  5 

1.148 

0.  366 

0. 1442 

0,0616 

0.0274 

0.01251 

0.00580 

0.00272 

0.001285 

0.000610 

1 . 925 

1.020 

0.  318 

0. 1221 

0.0508 

0.0220 

0.00974 

0.00438 

0. 001994 

0.000914 

0.000423 

2 .  0  2  5 

0.915 

0.280 

0. 1048 

0.0425 

0.01789 

0,00771 

0.00337 

0.001489 

0.000663 

0,000298 

2 .12  5 

0.  830 

0.248 

0.0909 

0.0359 

0.01475 

0.00618 

0.00262 

0. 001129 

0.000489 

0. 000214 

2.  2  25 

0.758 

0.222 

0.0797 

0.0307 

0.01228 

0.00501 

0.00207 

0.000868 

0.000366 

0.000157 

2 . 3  2  5 

0.  696 

0.200 

0.0703 

0.0265 

0.01033 

0.00411 

0,001652 

0.000675 

0.000278 

0.000116 

2.  425 

0.644 

0.  1820 

0,0625 

0.0230 

0.00876 

0,00340 

0,001332 

0. 000532 

0.000214 

0.000088 

2.  525 

0.598 

0. 1662 

0.0560 

0.0201 

0.00748 

0,00283 

0.001083 

0.000423 

0.000167 

0.000067 

2.  625 

0.  557 

0. 1524 

0.0503 

0.01170 

0.00643 

0.00238 

0.000887 

0.000339 

0.000131 

0.000053 

2.  '25 

0.522 

0. 1403 

0.0454 

0.01566 

0.00556 

0.00201 

0. 000733 

0. 000275 

0.000104 

0.000042 

2.  825 

0.490 

0.  1297 

0.0412 

0.01392 

0.00484 

0.001708 

0.000609 

0.000224 

0.000084 

0.000033 

2. 925 

0.461 

0.  1203 

0.0375 

0.01242 

0.00423 

0.001460 

0.000509 

0.000184 

0.000068 

0.000027 

3. 025 

0.435 

0. 1120 

0.0343 

0.01113 

0.00371 

0.001255 

0.000428 

0.000153 

0.000056 

0.000022 

3.  125 

0.412 

0. 1045 

0.0315 

0.01001 

0.00327 

0.001083 

0.000362 

0.000127 

0.000046 

0.000019 

3,  225 

0.391 

0.0977 

0.0289 

3.00904  ■ 

0.00289 

0.000939 

0.000307 

0.000107 

0.000039 

0.000016 

3.  325 

0.  371 

o\0916 

0.0267 

D. 00818 

0.00257 

0,000818 

0.000262 

0.000090 

0.000033 

0.000014 

3.425 

0,  353 

0.0860 

0,0247 

D.  00743 

0.00229 

0.000715 

0.000225 

0, 000077 

0.000028, 

0.000012 

3.  525 

0.  336 

0.0809 

0,0228 

0,00677 

0.00205 

0.000627 

0.000193 

0.000066 

,0.000024, 

0.000011 

3. 625 

0.321 

0.0763 

0.0212 

0.  00618 

0.001836 

0.000552 

0.000167 

0.000057 

0.000021 

0.000010 

3.  725 

0.306 

0.0720 

0.01973 

0.00565  , 

0.001651 

0,000488 

0.000145 

0.000049 

0.000018 

0.000009 

3. 825 

0.293 

0.0681 

0.01839  ! 

0.00518  * 

p. 001489 

0.000432 

0.000126 

0.000043 

0.000016 

0.000008 

3.925 

0,280 

0.0644 

0.01717.  ! 

0.00476 

p. 001345 

0.000384 

0.000110 

0.000037 

0.000015 

0,000008 

4,025  . 

0.269 

0.0610 

0.00439 

b. 001219 

0,000342 

0.000096 

0,000033 

0.000013 

0.000007 

.  4,125 

0.257 

0,0579 

0.00405 

b. 001107 

0.000306 

0.000085 

0.000029 

0.  000012  '1 

0.000007 

4.  225 

0.247 

;  0.0550 

0.  00374 

p. 001008 

0.000274 

0.000075 

0. 000026 

0.000011  i 

0.000006 

4.  325 

0.237 

0. 0523 

0.01323 

0.  00346 

0.000919 

0.000246 

0.000066 

0.000023 

O.OOOOlO 

0.000006 . 

4.  425 

0.228 

0.0497 

0.01244 

0.  00321 

), 000840 

0.000221 

0.000058 

0.000021 

0.000009  1 

0.000006 

4.  525 

0.219 

0.0473 

0,00298 

}. 000769 

0. 000200 

0.000052 

0.000019 

o.obbbbb 

0.000005 

4 . 625 

0.210 

0.0451 

0. 01102 

0.00277 

0.000706 

0.000180 

0.000046 

0. 000018 

6.000008 

0.000005 

4 . 3  .2  5 

0,202 

0.0430 

0.01039 

0.00258 

0,000649 

0.000164 

0,000041 

0.000016 

0.000008 

0.000005 

4.823 

0. 1944 

0.0410 

0.00241 

0.000597 

0. 000149 

0.000037 

0.000015 

0.000008 

0.000005 

4.925 

0. 1869 

0.0391 

0.00926 

0. 00225 

0.000551 

0.000135 

0.000033 

0.000014 

0. 000007 

0.000004 

5.02  5 

0.1798 

0.0373 

0.00875 

0.00210 

0.000509 

0.000123 

0,000030 

0.000013 

0,000007 

0.000003 

5.125 

0.1729 

0.0356 

0.00827 

0.001967 

0.000471 

0,000113: 

0.000027 

0.000012 

0.000006  I 

0.000003 

5.  225 

0. 1662 

0.0340 

0,00783 

0.001842 

0.000436 

3.000103 

0.000024 

0.000011 

0.000005 

0.000003 

5.  323 

0, 1598 

0.0325  : 

0.00741 

0.001727 

0.000405 

3. 000095 

0.000022 

0. 000011 

0.000005 

0.000003 

3.425 

0. 1536 

0.0310 

0.00702 

0.001621 

0.000376 

3.  000087 

0.000020 

0.000010 

■0.000004 

0.000002 

5.  525 

0. 1476 

0.0296 

0.00665 

0.001522 

0.000350 

0.000080 

0,000018 

0.000010 

(0.000004 

0.000002 

5 . 625 

0. 1417 

0.0282 

0.00630 

3.001430 

0,000326 

0.000074 

0.000017 

0.000009 

0.000003'^ 

0.000002 

1 

5.  725 

0.1359  ' 

0.0269 

0.00597 

).  001 345 

0,000304 

0.000069 

0.000015 

0.000009 

0.000003 

0.000001 

5. 825 

0,0257 

0.00565 

).  001265 

0.000284 

0.000063 

0.000014 

0.000008 

0.000002 

0.000001 

5.  725 

0. 1248 

0.0245 

0.00535 

0.  001191 

0,000265 

0.000059 

0.000013 

0.000008 

0.000002 

0.000001 

6.025 

0.1194 

0.0233 

0.00507 

0.001121 

0.000248 

0.000055 

0.000012 

0,000007 

0.000002 

0.000001 
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1.  IZS 
1  .  225 
1  .  325 
1.425 
1.525 


9*22 

6.30 

4.98 

4.20 

3.68 


7.13 

4,78 

3,72 

3.09 

2.66 


^3  4 


6.07 

3.99 

3.05 

2.49 

2.11 


5.  37 
3.47 
2.  60 
2.08 
1.741 


5 


4.85 

3.07 

2,26 

1.781 

1.461 


1  !  / 


4.43 

2.76 

1.990 

1,540 

1.241 


4.00 

2,48 

1.764 

1,340 

1.061 


2.  15 
1.556 
1.168 
0.910 


0.'9’86 

0.780 


2.33 
1.643 
1.222 
.0.  916 
0.614 


1.625 
1.  7  25 

1 .  825 
1.925 

2.  (>25 


3.  30 
3.01 
2,77 
2.58 
2.42 


2.35 

2.12 

1.928 

1.774 

1.645 


1.840 

1.632 

1,467 

1.332 

1.219 


1.492 
1.303 
1. 154 
1.033 
0.933 


1.232 

1.059 

0.923 

0.814 

0.724 


1.028 
0.868 
0.744 
0,  646 
0.566 


0.863 

0.716 

0.603 

0.515 

0.444 


0.727 

0.592 

0.491 

0,412 

0.350 


0.613 
0.491 
0.400 
0.  330 
0.275 


0.514 

0.408 

0.326 

0.265 

0.217 


2.125 
2.  2?.5 
2.  525 
2.425 
2.  525 


2.29 

2,17 

2.06 

1.968 

1,884 


1.535 

1,440 

1.356 

1.282 

1.216 


1. 124 
1..042 
0.970 
0.  907 
0.851 


0,848 

0.776 

0.714 

0.660 

0.612 


0.649 

0.586 

0.532 

0.485 

0.444 


0.500 
0.445 
0.  398 
0.359 
0.  324 


0.  387 
0.339 
0.299 
0.266 
0.237 


0.  300 
0.259 
0.225 
0. 1973 
0.1738 


0.232 
0. 1977 
0. 1696 
0. 1465 
0* 1274 


0. 1804 
0. 1512 
0.1278 
0. 1089 
0.0934 


2. 625 
2.  ,  25 
825 
2.  925 
3. 02d 


1.808 

1.739 

1.675 

1.617 

1.562 


1. 157 
1.103 
1.054 
1.009 
0.967 


0.802 

0.756 

0.716 

0.6.79 

0.644 


0.569 

0.532 

0.498 

0.467 

0,439 


0,409 
0.  377 
0.  349 
0.  324 
0.301 


0.295 

0.269 

0,246 

0.225 

0.207 


0.213 
0.1916 
0. 1733 
0. 1572 
0.1431 


0. 1539 
0. 1369 
0. 1223 
0. 1097 
0.0988 


0, 1113 
0. 0978 
0.0864 
0.0766 
0.0682 


0.0806 

0.0700 

0.0610 

0.0535 

0.0471 


3.125 
3. 225 
3.  325 
3. 425 
3,  525 


1.511 

1.464 

1.419 

1.377 

1.336 


0.928 
0.892 
0.  859 
0.  827 
0.798 


0.613 

0,584 

0.557 

0.532 

0.509 


0.413 
0.  390 
0.369 
0.349 
0.331 


0.281 

0.262 

0.246 

0.230 

0.216 


0.1914 
0. 1770 
0. 1641 
0. 1523 
0. 1417 


0.1306 
0. 1195 
0.1097 
0. 1009 
0.0930 


0.0892 

0.0808 

0.0734 

0.0668 

0,0610 


0.0609  0.0416 

0.0546  0.0369 


0.0491 

0.0443 

0.0400 


0.0329 

0,0294 

0.0263 


3. 625 
3.  r  25 
3.  825 
3.  925 
■4.025 


1.298 

1.262 

1.227 

1.194 

1.161 


0.770 

0.743 

0.718 

0.694 

0.671 


0,487 

0.467 

0.448 

0.429 

0,412 


0.314 

0.298 

0.284 

0.270 

0.257 


0.203 
0.1916 
0.1806 
0. 1705 
0. 1610 


0.1321 
0. 1233 
0. 1152 
0. 1078 
0. 1010 


0.0858 

0.0794 

0.0736 

0.0683 

0.0634 


0.0558 

0.0512 

0.0470 

0.0432 

0.0398 


0.0363 
0.0330 
0. 0300 
0.0274 
0.0250 


0. 0236 

0,0213 

0.01918 

0.01734 

0.01572 


3.  125 

4.  225 
•i.  325 
4.  125 


1.130 
1.  100 
1.071 
1,043 
1.015 


0.649 

0.628 

0.608 

0.588 

0,570 


0,396 
0.380 
0 ,  366 
0.352 
0.  338 


0.245 

0.234 

0.223 

0.213 

0.203 


0.1523 
0.1441 
0.1364 
0. 1292 
0.1225 


0.0947 

0.0890 

0.0836 

0,0786 

0.0740 


0.0590 

0.0550 

0.0512 

0.0478 

0.0447 


0.0368 

0.0340 

0.0314 

0.0291 

0,0270 


0.0229 

0.0210 

0.01927 

0.01772 

0.01630 


0.01427 

0.01298 

0.01182 

0.01079 

0.00985 


4.625 
4.7  25 
4.  525 
4,925 
025 


0.988 

0.961 

0,935 

0.910 

0.884 


0.551 

0,534 

0.516 

0.500 

0.483 


0.  325 
0.  313 
0.301 
0.289 
0.278 


0. 1941 
0. 1854 
0. 1771 
0. 1692 
0. 1616 


0.1162 
0. 1102 
0.1046 
0.0992 
0.0942 


0.0696 

0.0656 

0.0618 

0.0583 


0.0418 
0.0391 
0.0365 
0.0  342 


0.0549 


0.0320 


jo. 0250 
0.0233 
0.0216 
0, 0201 
0.01869 


0.01502 

0.01385 

0.01278 

0.01181 

0.01091 


0.00901 

0.00825 

0.00756 

0.00694 

0.00636 


5.  123 
215 
3.  325 
■'.4  25 
5.  525 


0.  859 
0.834 
0.809 
0.785 
0,760 


0.467 

0.451 

0.436 

0.420 

0.405 


0.267 

0.257 

0.246 

0.236 

0.227 


0. 1543 
0.1474 
0. 1406 
0. 1341 
0. 1278 


0.0894 

0.0848 

0.0804 

0.0762 

0.0722 


0.0518 

0.0488 

0.0460 

0.0433 

0.0408 


0.0300 

0.0281 

0.0263 

0.0246 

0.0231 


0.01740 

0,01619 

0.01507 

0.01402 

0.01304 


0. 01008 
0.  00933 
0,00862 
0.00798 
0.00737 


0.00585 
0.00537 
0. 00494 
0.00454 
0.00417 


^-:.625 
5.  r25 
3  ,  2  5 
5,  925 
('  .  0  25 


0.735 

0.710 

0.684 

0.658 

0.631 


0.  390 
0.  375 
0.  360 
0.  344 
0,  329 


0.217 
0.274 
0.1979 
0.1885 
0. 1791 


0. 1216 
0. 1156 
0.1097 
0.1039 
0, 0982 


0.0683 

0.0645 

0.0609 

0.0574 

0.0539 


0.0384 

0.0361 

0.0338 

0.0317 

0.0296 


0.0216 

0.0201 

0.01880 

0.01751 

0.01627 


0.01212 

0.01126 

0.01044 

0.00967 

0,00894 


0.00681 
0.00629 
0. 00580 
0, 00534 
0.00491 


0.00383 

0.00352 

0.00322 

0.00295 

0,00270 
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lABLI:,  17  Values  of  'p 


\ 

1 

1 

o 

-1 

- 

9 

10 

1*  ]0k  ' 

0.0144 

0,0304 

0.0479 

0.0672 

0.0884 

0.  1118 

o 

if- 

0. 166 

0. 1968 

0.231 

1 . 200 

0,0294 

0.0647 

0.  1079 

0. 1578 

0.219 

0.292 

0,380 

0,485 

0.  612 

0.763 

1  .  50C 

0.0447 

0. 1029 

0.1786 

0.277 

0.405 

0.571 

0.787 

1.068 

1.434 

1.908 

1. 

0,0606 

0. 1454 

0.264 

0.431 

0.663 

0.989 

1.446 

2.09 

2.98 

4.23 

1 .  ooo 

0,0769 

0.  1923 

0.  365 

0.625 

1.014 

1,599 

2.48 

3,79 

5.76 

8.72 

J  .  600 

0.0937 

0.244 

0.484 

0.  868 

1.482 

2.47 

4.04 

6.56 

10.58 

17.02 

1 .  -OU 

0.1111 

0.  300 

0.621 

1.167 

2.10 

3.67 

6.  36 

10.91 

18.66 

31,8 

1 , 800 

0, 1290 

0.361 

0,779 

1.532 

2.89 

5,33 

9.71 

17.61 

31.8 

57.4 

1  .  900 

0.1475 

0.428 

0.960 

1.972 

3.90 

7.55 

14.49 

27.7 

52.7 

100.4 

0.  000 

0, 1666 

0.500 

1. 167 

2.50 

5. 17 

10.50 

21.17 

42.5 

85.2 

'  170.5 

2.  100 

0. 1864 

0.578 

1,400 

3. 13 

6.75 

14.  37 

30.4 

63.9 

134.5 

283 

2. 20 0 

0.207 

0.  662 

1.663 

3.87 

8.71 

19,38 

42.  8 

94.4 

208 

458 

M)o 

0.228 

0.753 

1.959 

4.73 

11.  17 

25.8 

59.6 

137.2 

316 

727 

2.900 

0.250 

0.850 

2.29 

5.75 

14.04 

34.0 

81.7 

196.4 

472 

1132 

.  500 

0.273 

0.955 

2.66 

6.92 

17.57 

44.2 

110.8 

277 

693 

1734 

i.  600 

0,296 

1.067 

3.07 

8.28 

21.8 

57.0 

148.6 

387 

1005 

2610 

2.  700 

0.321 

1.  187 

3.53 

9.  84 

26.9 

72.9 

197.2 

533 

1439 

3880 

2.800 

0.  346 

1.315 

4.03 

11.63 

32,9 

92.5 

259 

726 

2030 

5700 

2.  7H:0 

0.  373 

1.453 

4.59 

13.67 

40.0 

116.4 

338 

981 

2840 

8250 

.5. 000 

0.400 

1.600 

5.20 

16.00 

48.4 

145.6 

437 

1312 

3940 

11810 

1.  ll'JO 

0.429 

1.757 

5.  88 

18.64 

58.2  ' 

180.9 

561 

1740 

5400 

16730 

3.  200 

0,458 

1.925 

6.  62 

21.6 

69.7 

224 

716 

2290 

7330 

23500 

5.  300 

0,489 

2.10 

7.43 

25.0 

83.1 

275 

907 

2990 

9870 

32600 

i.  iOO 

0.522 

2.  30 

8.33 

28.8 

98.6 

336 

1141 

3880 

13200 

44900 

3.  3G0 

0,556 

2,50 

9.31 

33.1 

116.5 

408 

1430 

5000 

17510 

61300 

j.  600 

0.591 

2,72 

10.38 

38.0 

137.2 

495 

1781 

6410 

23100 

83100 

i.  ;oo 

0.628 

2.95 

11.55 

43.3 

161,0 

596 

2210 

8170 

30200 

11180  , 

3.  800 

0.667 

3.20 

12.83 

:  49.4 

188,4 

717 

2720 

10350 

39300 

149500 

3.900 

0.707 

3.47 

14.22 

i  56.2 

220 

858 

3350 

13050 

50900 

198500 

4.  GUO 

0.750 

3.75 

15.75 

1  63.8 

256 

1024 

4100 

16380 

65500 

262000 

4.  100 

0.795 

4.05 

17.42 

'  72.2 

297 

1218 

4990 

20500 

83900  i 

344000 

9  .  200 

0,  842 

4.38 

19.23 

81.6 

364 

1444 

6070 

25500 

107000  1 

449000 

4.  300 

0.  892 

4.73 

21.2 

i  92.1 

397 

1708 

7350 

31600 

135800* 

584000 

4  .  400 

0.944 

5.10 

23.4 

;  103.8 

458 

2020 

8870 

39000 

171700 

755000 

4. 500 

1.000 

5.50 

25.8 

116.9 

527 

2370 

10680  ! 

48000 

216000 

973000 

4.600 

1.059 

5.93 

28.3 

131.4 

606 

2790 

12820  , 

59000 

271000  1 

4,  700 

1.121 

6.39 

31.2 

147.  6 

695 

3270 

15350  1 

72200  1 

339000 

4. 800 

1.188 

6.89 

34.3 

165.6 

1  796 

3820 

18350 

88000 

423000 

4. 900 

1.258 

7.42 

37.6 

185.6 

911 

4460 

21900 

107200 

525000 

i 

5.  (>00 

1.333 

8.00 

41.3 

208 

1041 

5210 

26000 

130200 

651000 

5.  100 

1.414 

8.62 

45.4 

233 

1189 

1  6070 

30900 

157800 

805000 

\  [ 

5.  200 

1.500 

9.  30 

49.9 

261 

1358 

7060 

36700 

190900 

i93000_^ 

; 

5.  300 

1.593 

10.03 

54.8 

292 

1549 

8210 

43500 

231000 

'  1 

S.  4  00 

1.692 

10.83 

60.2 

'327 

1766 

9540 

51500 

278000 

i 

5.  500 

1.  800 

11.70 

66, 2 

366 

2010 

11070 

60900 

335000 

5.600 

1.917 

12.65 

72.8 

409 

2290 

12850 

1  72000 

403000 

I 

5.  7  00 

2.04 

13.69 

80,1 

459 

2620 

14910 

85000 

484000 

5.  :'00 

2.18 

14.84 

88,2 

514 

2980 

1  17300 

100400 

582000 

5.  900 

2.33 

16.10 

97.3 

577 

3400 

20100 

118500 

699000 

6.  000 

2.50 

17.50 

107.5 

648 

3890 

23300 

140000 

840000 

, 

6.  100 

2.68 

19.06 

118.9 

728 

4440 

27100 

165400 

. 

6.  2C0 

2.89 

20.  8 

131.8 

820 

5090 

31600 

195600 

; 

6.  300 

3.  12 

22.8 

146,5 

926 

5840 

36800 

232000 

1 

6 ,  1-00 

3.  38 

25.0 

163.2 

1048 

6710 

42900 

275000 

6. 500 

3.67 

27.5 

182.4 

1189 

7730 

50300 

327000 

6.600 

4.00 

30.4 

205 

1355 

8940 

59000 

390000 

6 . 7  0  0 

4.  38 

33.8 

231 

1549 

10380 

69600 

466000 
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Value  8  of  A. 


\ 

\' 

1 

- 

'1 

G 

/ 

0 

JG 

I  .  100 

4.98 

2.98 

1.902 

1.244 

0.822 

0.546 

0.356 

0. 1696 

0.0808 

1  .  ^00 

3.28 

1.947 

1.237 

0.804 

0.528 

0.349 

0.229 

0.1035 

0.0488 

1.  300. 

2.61 

1.537 

0.971 

0.627 

0.410 

0.269 

0.1771 

O.llOS 

0.0834 

0.0420 

1.0  00 

2.30 

1.308 

0,  821 

0.527 

0.342 

0.224 

0.1464 

0.0949 

0,0709 

0.0365 

1,  300 

1.985 

1.159 

0,724 

0.462 

0.298 

0. 1933 

0. 1258 

0.0817 

0.  0444 

0.0322 

1  .  ^:00 

1.812 

1.052 

0.653 

0.415 

0.266 

0.1713 

0.1108 

0.0718 

0.0462 

0.0295, 

1  .  voo 

1.681 

0.600 

0.379 

0.241 

0. 1544 

0.0992 

0.0638 

0.0411 

0.0259 

1 . 800 

1.579 

0.908 

0.558 

0.350 

0.221 

0.1409 

0,0899 

0.0575 

0.0367 

0.0234 

1. 9G0 

1.495 

0.856 

0.523 

0.  326 

0.205 

0,1298 

0.0823 

0.0522 

0.0332 

0.0211 

0.  OOG 

1.426 

0.813 

0.494 

0,307 

0. 1917 

0. 1204 

0.0758 

0.0478 

0,0302 

0.01906 

0.  100 

1.  367 

0.776 

0.470 

0.290 

0. 1801 

0.1124 

0.0703 

0.0440 

0.0276 

0.01734 

Z,  ZOO 

1.  316 

0.744 

0.448 

0.275 

0.1699 

0.1054 

0.0655 

0.0408 

0.0254 

0.01585 

Z,  300 

1.272 

0.716 

0.429 

0.262 

0.1610 

0.0993 

0.0614 

0.0380 

0.0235 

0.01456 

400 

1.233 

0.691 

0.412 

0.250 

0.1531 

0.0939 

0.  0577 

0.0355 

0.0218 

0.01343 

0.  500 

1. 197 

0.669 

0.397 

0,240 

0, 1459 

0.0890 

0.0544 

0.0332 

0.0203 

0.01244 

600 

1.166 

0.649 

0.384 

0.231 

0,1395 

0.0846 

0.0514 

0.0312 

0.01901 

0.01156 

700 

1.137 

0.630 

0,  371 

0.222 

0. 1337 

0.0806 

0.0487 

0,0295 

0.01782 

0,01078 

2. 800 

1.  no 

0.613 

0.360 

0.214 

0. 1283 

0.0770 

0.0463 

0.0278 

0.01675 

0.01008 

Z.  900 

1.086 

0.598 

0.349 

0.207 

0. 1234 

0.0737 

0.0441 

0.0264 

0,01578 

0.00945 

■  5. 000 

1.063 

0.584 

0.  339 

0.200 

0. 1188 

0.0706 

0.0420 

0.0250 

0.01490 

0.00888 

1.  lf>0 

1.042 

0.570 

0,  330 

0. 1941 

0. 1146 

0.0678 

0.  0402 

0.0238 

0.01410 

0.00836 

3.  ZOO 

1.023 

0.558 

0.  322 

0. 1883 

0. 1107 

0.0652 

0.0384 

0.0227 

0.01336 

0.00788 

3.  300 

1.004 

0.546 

0.  314 

0. 1829 

0.1070 

0.0627 

0.0368 

0.0216 

0.01269 

0.00745 

3. 400 

0.987 

0.535 

0.  306 

0.1778 

0. 1036 

0.0605 

0.0353 

0.0206 

0. 01206 

0.00705 

o 

o 

0.971 

0.524 

0.299 

0. 1730 

0, 1004 

0.0583 

0.0339 

0.01974 

0. 01149 

0.00669 

3.  600 

0.955 

0.514  . 

0.293 

0. 1684 

0.0973 

0.0563 

0.0326 

0.01890 

0.01095  . 

0.00635 

3.  700 

0.940 

0.505 

i  0.286 

0.1641 

0.0945 

0.0544 

0.0314 

0.01812 

0.01046 

0.00604 

3.  800 

0.926 

0.496 

0.280 

0. 1600 

0.0917 

0,0527 

0.0303 

0.01739 

0.01000 

0.00575 

3.  900 

0.912 

0.487 

0.274 

0.1561 

0.0892 

0.0510 

0.0292 

0.01670 

0.00956 

0.00548 

4.000 

0.899  i 

0.479 

0.269 

0. 1524 

0.0867 

0.0494 

0.0282 

0.01606 

0.00916 

0,00522 

4.  100 

0.  887  ; 

0.471 

0.263 

0.1489 

0.0844  ' 

0,0479  1 

0. 0272 

0.01545 

0.00878 

0,00499 

4 ,  ZOO 

0,874 

0.463 

0.258- 

6. 1455 

0.0822  ; 

0.0465  : 

0,  0263 

0.01488 

0. 00842 

0.00477 

4.  300 

0.863 

0.456 

0.253 

0. 1422 

0.0801 

0.0451  ‘ 

0.0254 

0,01434 

0.00809 

0.00456 

4.4  00 

0.851 

0,449 

0.249 

0, 1391 

0.0780 

0.0438 

0,0246 

0.01383 

0.00777 

0,00437 

4.  SOO 

0.  840 

0.442 

0.244 

0. 1361 

0.0761 

0.0426 

0.0238 

0.01335 

0.00748 

0.00419 

•i,  600 

0.  829 

0.435 

0.240 

0. 1332 

0.0742 

0,0414 

0.0231 

0,01289 

0.00719 

0.00402 

4.  700 

0.819 

0.428 

0.235 

0.1304 

0.0724 

0.0403 

0,0224 

0.01245 

0.^0693 

0.00385 

4.  MOO 

0.  808 

0.422 

0.231 

0. 1277 

0.0707 

0,0392 

0,0217 

0.01204 

0.00668 

0.00370 

4.  900 

0.798 

0,416 

0.227 

0. 1250 

0.0690 

0.0381 

0.0211 

0.01164 

0.00644 

0.00356 

.5. 000 

0.788 

0.410 

0.  223 

0.1225 

0.0674 

0.0371 

0.0204 

0.01126 

0,00621 

0.00342 

5.  100 

mmm 

0,404 

0.219 

0. 1200 

0.0658 

0.0361 

0.01985 

0.01090 

0.00599 

0.00329 

5.  200 

nSM 

0.  398 

0.215 

0. 1176 

0.0643 

0.0352 

0.01928 

0.01056 

0.00578 

0,00317 

8.  300 

0.759 

0.  392 

0.212 

0. 1153 

0.0629 

0.0343 

0.01873 

0.01023 

0.00558 

0.00305 

6. 4U0 

0.750 

0.386 

0.208 

0.1130 

0.0614 

0,0334 

0.01820 

0.00991 

0.00539 

0.00294 

5.  500 

0.741 

0.381 

0.205 

0.1108 

0,0601 

0.0326 

0.01769 

0,00960 

0. 00521 

0.00283 

o 

o 

0.731 

0.375 

0.201 

0. 1086 

0.0587 

0.0318 

0.01719 

0.00931 

0.00504 

0.00273 

5.  700 

0.722 

0.370 

0. 1976 

0.1064 

0.0574 

0.0310 

0.01672 

0.00902 

0.00487 

0.00263 

5 .  800 

0.713 

0.364 

0»1942 

0. 1043 

0.0561 

0.0302 

0.01625 

0,00875 

0,00471 

0.00254 

5.  900 

0.703 

0.359 

0.  1908 

0. 1022 

0.0548 

0.0294 

0.01580 

0.00848 

0.00456 

0.00245 

6.000 

0.694 

0.353 

O. 1875 

0. 1002 

0.0536 

0.0287 

0.01537 

0.00823 

0.00441 

0.00236 

6.  100 

0.685 

0.348 

0. 1842 

0.0982 

0.0524 

0.0280 

0.01494 

0.00798 

0.00426 

0,00228 

6.  ZOO 

0.675 

0.342 

0. 1808 

0.0962 

0.0520 

0.0273 

0.01452 

0.00774 

0.00412 

0.00220 

6.  300 

0.  666 

0.337 

0.  1775 

0.0942 

0,0500 

0.0266 

0.01412 

0.00750 

0.00399 

0.00212 

6.400 

0.  656 

0.331 

0. 1742 

0.0922 

0.0488 

0.0259 

0.01372 

0.00727 

0.00385 

0.00204 

0.  646 

0.  326 

0. 1709 

0.0902 

0.0477 

0.0252 

0.01332 

0.00705 

0,00373 

0.001971 

6.  600 

0.636 

0.320 

0. 1675 

0.0882 

0.0465 

0.0245 

0.01294 

0.00683 

0.00360 

0.001900 

6.  VOO 

0.626 

0.314 

0. 1641 

0,0862 

0.0454 

0.0239 

0.01256 

0.00661 

0.00348 

0.001831 
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Values  of  T 


/X  - 


1 

i 

1 

G 

1  ; 

1  . 

9*24 

3.61 

1.756 

0.930 

0.516 

0.294 

0. 1652 

0. 1937 

0.  0685 

0,0247 

i  .  .!00 

5.07 

1.  897 

0.  883 

0.482 

0.238 

0. 1304 

0.0711 

0.0855 

0.0288 

0.00721 

1  .  >v0 

3.50 

1.261 

0.566 

0.277 

0. 1419 

0.0749 

0,0401 

0.01926 

0.01949 

0.00445 

1.  iv-r 

2,  66 

0.929 

0.403 

0. 1905 

0.0944 

0.0481 

0.0249 

0.01272 

0.  01505 

0.00277 

1  .  SCO 

2.  14 

0.726 

0.305 

0, 1398 

0.0670 

0.0331 

0.01659 

0.00834 

0.  01667 

0,001796 

1 . 

1.791 

0.590 

0,241 

0. 1071 

0,0498 

0.0238 

0.01155 

0.00567 

0.00275 

0.001325 

1  .  ,00 

1.535 

0.493 

0.1958 

0.0845 

0.0381 

0.01764 

0,00830 

0.00395 

0,001886 

0.000859 

i  .  ^0i) 

1.341 

0.420 

0. 1626 

0.0683 

0.0299 

0.01343 

0.00612 

0.00283 

0..001311 

0.000605 

\  .  OOU 

1.  189 

0.  364 

0. 1374 

0.0562 

0.0239 

0.01043 

0.00461 

0.00207 

0.000931 

.0.000422 

ooo 

1.067 

0.320 

0.1177 

0.0469 

0.01943 

0,00824 

0.00354 

0.001540 

0.000674 

0.000298 

100 

0.  967 

0.284 

0. 1021 

0.0396 

0.01600 

0. 00660 

0.00275 

0.001166 

0.000497 

0.000214 

:'ou 

0.  883 

0.  254 

0.0.893 

p.0339 

0.01332 

0. 00535 

0.00217 

0.000895 

0. 000371 

0.000156 

.A  ooo 

0.812 

0.229 

0.0789 

0.0292 

0.01120 

0,00438 

0.001732 

0.000696 

0, 000282 

0.000116 

. , .  •  1 0  0 

0.751 

0.208 

0.0701 

b.0254 

0.00949 

0.00362 

0.001396 

0.000548 

0. 000217 

0,000088 

.  sou 

0.698 

0, 1899 

0.0627 

0.0222 

0.00811 

0.00302 

0.001135 

0.000436 

0, 000169 

0.000067 

0.  600 

0.651 

0. 1743 

0,0564 

0.01952 

0.00697 

0.00253 

0.000931 

0.000350 

0,000133 

0.000052 

0.  voo 

0.610 

0. 1606 

0.0510 

0.01727 

0.00603 

0.00214 

0.000769 

0. 000283 

0.000106 

0.000041 

0.  ?.oo 

0.574 

0. 1486 

0,0463 

0.01535 

0.00525 

0.001823 

0.000639 

0.000231 

0. 000085 

0.000033 

0,  000 

0.541 

0. 1380 

0.0422 

0,01371 

0.00459 

0,001560 

0.000535 

0.000190 

0.000069 

0.000027 

3.  000 

0.512 

0. 1285 

0.0386 

0.01230 

0.00403 

0.001341 

0.000450 

0,000158 

0. 000057 

0.000022 

A  i<Vo 

0.485 

0. 1200 

0.0354 

b. 01107 

0.00356 

0. 001159 

0.000381 

0,000132 

0. 000047 

0.000019 

3.  .300 

0.461 

0. 1124 

0.0326 

0.01000 

0.00315 

0,001006 

Q. 000324 

0.000111 

0. 000039 

0.000016 

3.  300 

0,438 

0. 1055 

0.0301 

0.00906 

0. 00280 

0.000877 

0.000277 

0. 000094 

0.000033 

0.000014 

100 

0.418 

0.0992 

0.  0278 

p, 00824 

0.00250 

0.000767 

0.000238 

0.000080 

0. 000028 

0.000012 

A  300 

0.  399 

0.0935 

0.0258 

0.00751 

0.00234 

0. 000674 

0.000205 

0. 000068 

0, 000024 

0.000011 

/.  600 

0.  381 

0.0883 

0.0240 

0.00686 

0.00201 

0.000594 

0.000177 

0,000059 

0. 000021 

0.000010 

..  ,  00 

0.  365 

0.0835 

0.0224 

0.00629 

10.001805 

0.000525 

0.0001541 

0,000051: 

'  0. 000018 

0,000009 

3 .  -HOG 

0. 350  . 

0.0791 

0.0209 

0.00577 

10.001629 

0,000465 

0.000134| 

0.000044!  0,  000016 

10.000008 

,3.  900 

0.  336 

0.0750 

0.01951 

0,00531 

0.001474 

0. 000414 

0.000117 

0.000039i0. 000015 

0.000008 

4.  000 

0.  323 

0.0712 

0.  01827 

0.00489 

0.001337 

0.000369 

0,000102 

0.000034;.  0.000013 

0.000007 

4.  100 

0.310 

0.0677 

0.01714 

0,00452 

0.001215 

0.000330 

0. 000090 

0.000030 

0. 000012 

0.000007 

4  .  200 

0.298 

0.0645 

0.01610 

0.00418 

0.001107 

0.000296 

0.000079 

0.000027 

0.000011 

0.000006 

4.  300 

0.287 

0.0614 

0.01514 

0.00388 

0.001010 

0.000266 

0.000070 

0.000024 

0, 000010 

0.000006 

•t.  iOO 

0.277 

0.0586 

0.01425 

0.00360 

0.000924 

0.000239 

0.000062 

0.000022 

0. 000009 

0.000006 

4.  GOO 

0.267 

0.0559 

0.01344 

0, 00335 

0. 000846, 

0.000216 

0.000055 

0.000020 

0.000009 

0.000006 

4.^00 

0.  258 

0.0534 

0.01269 

0.00312 

0.000777 

0.000195 

0.000049 

0.000018 

0.000008 

0.000006 

4.  ,'00 

0.249 

0.0511 

0.01199 

0.00291 

0.000714 

0.000177 

0.000044 

0.000017 

0.000008 

0.000005 

4.000 

0.240 

0.0489 

0,01134 

0.00271 

0. 000658 

0. 000161 

0.000039 

0.000015 

0. 000007 

0.000005 

.  voo 

0.232 

0.0468 

0,01073 

0.00254 

0.  000607 

0.000146 

0.000035 

0.000014 

0. 000007 

0.000004 

0  ,  (>  0  0 

0.224 

0.0448 

0.01017  \ 

0.00237 

0.000561 

0.000133 

0.000032 

0.000013 

0. 000007 

0.000004 

3.  100 

0.217 

0.0430 

0.00965  j 

0,00223 

0,000519 

0.000122 

0.000028 

0.000012 

0,000007 

0.000004 

4.  400 

0.209 

0.  0412 

0,00915 

0.00209 

0.000481 

0,000111 

0.000026 

0,000012 

0.000006 

0.000004 

6.  300 

0.203 

0.0395 

0.  00869 

0.001961 

0.000446 

0.000102 

0.000023 

0.000011 

0. 000006 

0.000003 

3 . 4  GCf 

0. 1959 

0.0379 

0.00826 

0.001843 

0.000415 

0.000093 

0.000021 

0.000010 

0.000006 

0,000003 

•3.  :';GC 

0, 1895 

0.0364 

0.00786  ! 

0.001743 

0,000386 

0.000086 

0.000019 

0.000010 

0. 000006 

0.000003 

oOO 

0. 1834 

0.0350 

0.00748 

0.001634  i 

0.000359 

0. 000079 

0.000017 

0.000009;  0.000005 

0,000003 

T-  .  F  0  ’  ' 

0. 1774 

0,0336 

0.00712 

0.001540 

0.000335 

0.000073 

0.000016 

0.000009 

0. 000005 

0.000003 

5.  '-'00 

0. 1716 

0.0323 

0.00678  1 

0.001453 

0.000313 

0.000067 

0.000014 

0.000009 

0.000005 

0.000002 

0.  000 

0. 1660 

0.0310 

0.00646 

0.001372  1 

0.000293 

0. 000062 

0.000013 

0.000008 

0.000005 

0.000002 

0.  r;C0 

0. 1605 

0.0298 

0.00616  1 

0,001296  1 

0. 000274 

0.000058 

0.000012 

0,000008 

0.000004 

0.000002 

0.  100 

0. 1552 

0.0286 

0.00587  1 

0.001226  1 

0.000257 

0.000054 

0.000011 

0.000007 

0. 000004 

0.000002 

,  '  >n  '* 

0, 1500 

0.0275 

0.00560 

0.001160  1 

0,000241 

0.000050 

0,000010 

0. 000005 

0. 000003 

0.000002 

0.  .30(' 

0. 1449 

0,0264 

0.00535  ( 

0.001098  1 

0.000226 

0,000046 

0.000009 

0.000005 

0. 000003 

0.000002 

6.  I'-O 

0. 1400 

0.0254 

0.00510 

0.001040  1 

0. 000213 

0.000043 

0.000008 

0.000005 

0. 000003 

0.000002 

(  .  000 

0, 1351 

0.0244 

0.00487  1 

0.000986  1 

0. 000200 

0.000040 

0.000008 

0.000004| 

0,000003 

0.000001 

(■.0  0:j 

0. 1303 

0.0234 

0.00464  ( 

0.000934  1 

0,000188 

0. 000038 

0.000007 

0. 000004 

0. 000003 

0.000001 

0. 1256 

0.0225 

0.00443  ( 

0.000886  1 

0.000177 

0.000035 

0.000007 

0.000004 

0. 000002 

0.000001 
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1  .  iOO 

1 .  200 
]  .  300 
1  .  100 
1  .  rOO 
1 . 600 

1 .  roo 
r.  000 

1 . 900 
2.000 
J.  100 

2.  ^00 
2.  300 

2.400 
2.  500 
2.600 
2.  700 
2.  800 

2,  900 
3. 000 

3.  lf)0 
3.  200 
3.  300 

3. 400 
3.  500 
3.  600 
3.  700 

3. 800 

3.  000 
4.000 
4.100 

4.  ZOO 
4.  300 

4.400 

4.  500 

4. 600 
4.700 

4. 800 
4. 000 
5. 000 

5.  100 
5.  200 
5.  300 

5. 400 
5.  500 
5.  600 

5.  700 
5.S00 

5. 900 
6.000 

6.  100 
6.  200 
6.  300 

6.400 
6.  500 

6.600 
6.  700 


11.24 
7.  17 
5,54 

4.62 
4,02 
3.59 
3.26 
3.00 
^.79 

2.62 
2.47 
2.  34 
2.23 
2.13 
2.04 
1.961 
1.888 
1.821 
1.760 
1.703 
1.650 
1.601 
1.555 
1.511 
1.470 
1.431 
1.394 
1.359 
1.325 
1.293 
1.262 
1.232 
1,204 
1. 176 
1. 149 
1.123 
1.098 
1.073 
1.049 
1.025 
1.002 
0.980 
0.957 
0.936 
0,914 
0,  893 
0,  871 
0.  850 
0.  830 
0,809 
0,788 
0.767 
0.746 
0.725 
0.704 
0.683 
0.661 


8,74 

5,47 

4.15 

3.40 

2,91 

2.57 

2.30 

2.09 

1.923 

1.782 

1.663 

1,560 

1,470 

1.390 

1.320 

1.256 

1.199 

1.147 

1.099 

1.055 

1.014 

0.976 

0.941 

0.908 

0.877 

0.  848 

0.821 

0.795 

0.770 

0.747 

0.724 

0.703 

0.683 

0.663 

0.644 

0.626 

0.609 

0.592 

0.575 

0.560 

0.544 

0.529 

0.515 

0.501 

0.487 

0.473 

0.460 

0,447 

0.434 

0,421 

0.408 

0,396 

0.383 

0.371 

0.359 

0.346 

0.334 


1,780 

1.597 

1.448 

1.325 

1.221 

1.131 

1.054 

0.986 

0.926 

0.872 

0.824 

0.780 

0.740 
0.704 
0.671 
0.640 
0.612 
0.585 
0.561 
0.538 
0.516 
0.496 
0.477 
0.459 
p.442 
0.426 
0.411 
0,396 
0,  383 
0.369 
0.357 
0.345 
0.  333 
0,322 
0.311 
0.301 
0.291 
0.281 
0.272 
0.263 
0,254 
0.246 
0.237 
0.229 
0.221 
0.213 
0.206 
0, 1981 
0. 1906 
0. 1831 
0.1757 


6.64 

4,00 

2.92 

2.32 

1.921 

1.639 

1.427 

1,261 

1,127 

1.017 

0.924 

0,846 

0.778 

0.719 

0.667 

0.621 

0,580 

0,543 

0.510 

0.480 

0.453 

0,428 

0.405 

0.384 
0.365 
0.  347 
0.330 
0.315 
0.300 
0.287 
0.274 
0.262 
0.251 
0.240 
0.230 
0.221 
0.212 
0.203 
0.1951 
0. 1874 
0, 1800 
0, 1730 
0.1662 
0.1598 
0. 1536 
0.1476 
0. 1418 
0.1363 
0.  1309 
0.1258 
0.  1207 
0. 1158 
0.1111 
0.1064 
0.1019 
0.0974 

0.0930 


6.03 
3,56 
2.55 
1.989 
1.620 
1.359 
1.163 
1.012 
0.891 
0.792 
0.710 
0.  640 
0.581 
0.530 
0,486 
0.447 
0.413 
0.382 

0.355 
0.331 
0.  309 
0.289 
0.271 
0,254 
0.239 
0.225 
0.213 
0.201 
0. 1899 
0. 1799 
0. 1705 
0. 1618 
0.1537 
0.1461 
0. 1390 
0. 1323 
0.1261 
0.1202 
0.1146 

0.1093 

0.1043 

0.0996 

0.0951 

0.0909 

0.0868 

0.0830 

0.0793 

0.0757 

0.0724 

0.0691 

0.0660 

0.0630 

0.0601 

0.0573 

0.0545 

0.0519 

0.0493 


5.53 
3,20 
2.26 
1,726 
1.381 
1, 138 
0.958 
0.819 
0.710 
0.621 
0.549 
0.488 
0.437 

0,393 

0.356 

0.323 

0,295 

0.270 

0,248 

0.228 

0.211 

0.1952 

0.1811 

0. 1684 

0. 1570 

0,1465 

0.1370 

0. 1283 

0.1204 

0.1130 

0. 1063 

0.1001 

0.0943 

0.0890 

0,0840 

0.0794 

0.0751 

0.0711 

0.0674 

0.0638 

0,0605 

0.0574 

0.0545 

0.0517 

0.0491 

0.0467 

0.0443 

0,0421 

0.0400 

0.0380 

0.0361 

0.0343 

0.0325 

0.0308 

0.0292 

0.0277 

0.0261 


4,93 
2,  85 
2.00 
1.507 
1. 185 
0.959 
0.793 
0.667 
0.568 
0.490 
0.426 
0.  373 
0.329 
0.292 
0.261 
0.234 
0.211 
0. 1908 

0. 1732 
0.1578 
0. 1442 
0.1321 
0.1214 
0.  1118 
0. 1032 

0.0954 

0.0884 

0.0821 

0.0763 

0.0711 

0.0663 

0.0619 

0,0579 

0.0542 

0.0508 

0.0477 

0.0448 

0.0421 

0,0396 

0.0330 

0.0351 

0.0331 

0,0312 

0,0295 

0.0278 

0.0263 

0^0248 

0.0234 

0.0221 

0.0209 

0.01974 

0.01864 

0.01759 

0.01660 

0.01565 

0.01474 

0.01386 


.659 
.298 
.015 
.811 
.659 
,545 
.456 
.387 
1.331 
0.286 
0,249 
0,218 
0. 1917 
0.1698 
0.1511 
0.1351 

0.1212 
0.1092 
0.0987 
0.0895 
0.0814 
0.0742 
0,0679 
0.0622 
0.0571 
0.0525 
0.0484 
0.0447 
0,0414 
0.0384 
0.0360 
0.0331 
0.0308 
0.0287 
0.0267 
0.0249 
0.0233 
0.0218 
0.0204 
0.01911 
0.01790 
0,01679 
0,01575 
0,01478 
0.01388 
0.01303 
0.01224 
0.01150 
0.01080 
0.01014 
0.00952 
0.00893 
0.00  838 
0.00785 
0.00735 


4.  11 
2.25 
1.617 
1.214 

0.  68r^ 

0,548 

0.445 

0.  367 

0,306 

0.258 

0.219 

0.1879 

0.1622 

0,1410 

0.1232 

0.1083 

0.0957 

0.0849 
0.0756 
0.0676 
0.0607 
0.0546 
0.0493 
0.0446 
0.0405 
0.0369 
0.0336 
0. 0308 
0.0282 
0,0258 
0.0238 
0.0219 
0.0202 
0.01862 
0.01722 
0.01595 
0.01478 
0.01372 
0.01274 
0.01185 
0.01102 
0.01026 
0.00957 
0.00892 
0.00832 
0.00777 
0.00725 
0.00677 
0.00633 
0.00591 
0.00552 
0.00515 
0.00481 
0.00449 
0.00419 
0.00390 


2.45 
.1.429 
1.168  , 
0.935 
0.745 
0.605 
0. 447 
0.363 
0.295 
0.243 
0.201 
0. 1683 
0. 1421 
0. 1210 
0.1037 
0.0895 
0,0777 
0.0678 

0.0594 
0.0524 
0.0463 
0.0411 
0.0366 
0.0328 
0.0294 
0.0264 
0.0238 
0.0215 
0.01952 
0.0  1773 
0.01614 
0.01471 
0.01344 
0.Q1230 
0.01128 
0.01035 
0.00951 
0.00876 
0.00807 
0.00745 
0,00688 
0.00636 
0.00588 
0.00545 
0.00505 
0.00468 
0.00435 
0.00403 
0,00375 
0.00348 
0.00323 
0.00300 
0.00279 
0.00259 
0,00241 
0.00223 
0.00207 
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APPENDIX  A 


FLUID  DYNAMIC  MODEL 

Experience  has  shown  that  useful  results  for  the  behavior  of  solids 
under  high  stress  intensities  can  be  obtained  under  the  assumption  of  stress 
isotropy  and  thermodynamic  equilibrium^.  This  is  equivalent  to  assuming 
that  the  material  behaves  as  an  ideal  compressible  fluid.  Such  behavior 
in  a  blast  absorber  can  be  approximately  insured  by  suitable  selection  of 
material.  In  this  appendix,  we  indicate  the  character  of  the  approximations 
involved  in  a  fluid  dynamic  model  of  a  non-linear  solid. 

Generally,  one  expects  that  a  complete  thermodynamic  description 
of  a  material  would  involve  a  highly  complex  equation  of  state  comprising 
a  relation  among  many  variables  such  as  density,  temperature,  internal 
energy,  stress,  and  deformation  rate.  Alternatively  the  stress  is  regarded 
as  a  function  of  a  number  of  other  state  variables,  thus  defining  a  stress- 
deformation  law  involving  variables  in  addition  to  density  or  strain.  This 
would  be  the  general  "stress- strain"  law,  different  approximations  to 
which  would  be  expected  to  hold  in  various  applications.  We  can,  however, 
discuss  properties  of  the  stress  without  regard  to  a  general  law,  in  order 
to  define  what  is  meant  by  a  fluid  dynamic  approximation. 

In  Cartesian  coordinates  Xj^,  x^,  x^  the  stress  condition  at  a 
point  in  a  body  is  represented  by  the  tensor  ,  a  set  of  nine  numbers 

(i,  j  =  1,  2,  3)  which  are  the  stress  components: 


1  ctj! 

^12 

^13  \ 

a. .  = 

^^21 

^22 

^23 

\  ^31 

^32 

^33  / 
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The  components  are  interpreted  as  tensile  and  shear  forces  per  unit  area 
acting  on  the  sides  of  an  infinitesimal  cubical  element.  Thus  ^ 
force  in  the  ~  direction  per  unit  area  acting  on  a  plane  whose  normal  is 
in  the  positive  Xj  -  direction.  The  stress  tensor  can  be  shown  to  have  the 
symmetry  property  (7..  =CT..;  thus  there  are  in  general  only  six  independent 
stress  components. 

We  can  separate  the  stress  tensor  into  so-called  isotropic  and 
deviatoric  parts  by  writing 


a..  =  -  P  a'--  +  a'..  , 

ij  ij  ij 

where  S--  is  the  Krone  eke  r  delta,  and  the  mean  normal  pressure  P  is 

defined  by  P  =  -  3  .  Here  the  repeated  subscript  indicates  summation 

over  the  values  1,  2,  3  of  the  index,  amd  P  is  therefore  the  negative 

14 

average  normal  stress. 

A  hydrostatic  state  of  stress  in  any  material  is  by  definition  such 
that  all  components  of  the  stress  deviator  vanish;  in  such  a  case 

O'lj  ~  ^ 2.Z  ~  ^33  ~  ^  simply  referred  to  as  the  hydrostatic 

pressure.  We  are  interested  in  materials  which  behave  so  that,  as  the 
magnitude  of  a  general  stress  state  is  increased,  the  magnitude  of  the 
stress  deviator  becomes  and  remains  small  compared  to  the  mean  normal 
pressure.  Then  the  state  of  stress  at  sufficiently  high  stress  magnitudes 
becomes  nearly  hydrostatic.  The  stress  magnitude  is  evaluated  in  terms 
of  the  scalar  invariant  O..  (J..  : 


=  3  P'^  +  a'. .a'..  , 

ij  Ji  ij  Ji 

since  by  definition  =  0.  This  invariant  of  0^^  is  therefore  the  sum 
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of  the  invariants  of  the  isotropic  and  deviatoric  parts  of  For  a 

linear  elastic  material 5  can  be  shown  to  be  proportional  to  the 

J  ^ 


energy  of  distortion  or  change  of 


15 

of  the  material  .  An  approximately 


a*.,  a’..  «  3P  . 
iJ 


isotropic  state  of  stress  is  defined.' 

In  many  applications^  one  is  interested  in  material  behavior  in  a 
state  of  uniaxial  strain  such  as  one « dimensional  fully  confined  compression 
in,  say,  th^  x^  =  directiono  For  example,  the  uniaxial  strain  condition 
is  approached  at  large  distances  from  a  point-source  explosion.  Referred 
to  principal  axes  of  deformation,  the  (symmetric)  strain  tensor  .  is 
such  that  ^  Q  while  all  other  components  0.  Then  by  symmetry 

~  ^33^  ®  i  ^  jo  For  this  case 


Now  it  is  easily  shown  that  twice  the  maximum  shear  stress 

in  the  material  for  this  stress  state.  If  theshear  stress  is  limited  by 

some  maximum  value  Q  ,  a  property  of  the  material,  then  a'. .a'..  =8a^/3 

J  ^ 

for  sufficiently  high  or||«  The  mean  normal  pressure  for  this  case  is 


1 

3 


■I'  Z 


11 


3  ^^22  ■  ^11 


). 


and  P  =  -  Cji  -4a/3  when  the  maximum  shear  stress  is  reached.  There¬ 
fore  P  —  “  sufficiently  high  compressions,  and  the  state  of  stress 

becomes  nearly  hydrosta^tic. 

Unfortunately,  both  0  and  P  are  likely  to  be  rate -dependent 

'■'V 

quantities.  For  easily  crushable  solids  it  is  plausible  to  assume  that 
the  rate  dependence  of  0  is  not  serious.  In  order  to  make  possible  the 
formulation  of  an  ideal  fluid-type  equation  of  state  applicable  at  high  stress 
magnitudes,  it  is  also  necessary  to  assume  that  the  rate  dependence  of  P 


A-3 


Report  No. 
ARF  4132-6 


is  embodied  solely  in  changes  of,  say,  temperature  and  internal  energy, 
and  to  discard  any  explicit  dependence  of  P  on  rate  of  straim  To 
illustrate  this  point,  consider  a  stress -strain- strain  rate  relation  of  the 


following  type : 


a..  -  (a..)  S--  +  2  jxe..  +  X'e,,  S--  +  2|x'€..  (A-l) 

ij  '  ij'  kk  ij  kk'-'ij  ^  ij  '  ' 


where  the  strain  rate  tensor  €..  is  defined. by 

ij 


(Uo  .  +  u.  .), 
^  i.J 


being  the  velocity  vector^  The  subscript  zero  indicates  a  reference 
state  of  rest  at  zero  strain.  Subscripts  following  a  comma  denote  partial 
differentiation  with  respect  to  the  indicated  coordinate,  and  the  overhead 
dot  denotes  total  differentiation  with  respect  to  time  following  the  motion. 
The  natural  strain  ,  so  defined  in  terms  of  its  derivative  following  the 
motion,  is  not  restricted  to  infinitesimal  values.  For  the  moment,  the 
four  quantitites  A  ,  jx’  are  assumed  to  be  constant  properties  of  the 

material.  Equation  (A“l)  above  represents  the  simplest  possible  such 
relation  for  an  isotropic  material,  since  the  strain  and  strain-rate  depen¬ 
dencies  are  linear  and  aAditiveo  For  infinitesimal  strains  in  the  absence 


of  rate  effects,  A  and  jx  would  be  the  Lfame  constants  for  an  isotropic 
elastic  solid.  If  the  material  were  an  ideal  fluid  which  cannot  support  shear 
stress  statically,  A  (S  -  +  2  |x  ~  -  (P  -  P^)  necessarily,  and  A* 

and  jx*  would  be  the  viscosity  coefficients  for  a  Newtonian  fluid.  Equation 
(A-l)  may  be  interpreted  as  a  possible  approximate  relation  describing  a 
material  at  stress  levels  and  loading  rates  somewhat  greater  than  those  at 
which  a  rate -independent  linear  law  would  apply.  The  mean  normal  pressure 


for  this  case  is 
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^>  ]  - 


(X  +  3  F)€j^ 


From  kinematic  (continuity)  considerations  alone,  it  can  be  shown  that 

€  =  -  p  /  p  and  €  =  -  In  (p/p^)*  where  p  is  the  mass  density. 

Equation  (A-2)  thus  represents  a  relation  among  pressure,  density,  and 

rate  of  change  of  density.  Now,  to  discard  the  explicit  rate  of  dependence 

2 

of  P  at  this  point  would  be  equivalent  to  the  assumption  X'  3  ®  » 

however,  P  can  still  be  considered  as  implicitly  rate -dependent  if  X  and 
are  allowed  to  be  functions  of  state  (e,  g.  ,  of  temperature  and  internal 
energy)  rather  than  constants.  This  is  exactly  analogous  to  the  Stokes 
hypothesis,  which  assumes  that  X'  +  fJ-'  =  0  for  an  ideal  fluid;  only 
under  this  condition  can  the  thermodynamic  pressure  in  a  moving  fluid 
be  identified  with  the  thermostatic  pressure,  i.  e.  ,  that  which  obtains  when 
the  same  state  is  reached  quasi -statically.  Assumptions  of  this  general 
nature  are  implicit  in  the  characterization  of  a  solid  at  high  pressures  by 
a  fluid  equation  of  state. 

To  recapitulate,  the  approximations  involved  in  a  fluid  dynamic 
model  of  a  solid  at  high  pressures  are 

1.  The  stress  tensor  is  nearly  isotropic;  that  is,  the  magnitude 
of  the  stress  deviator  is  small  compared  to  the  square  of  the  mean  normal 


pressure. 

2.  States  of  the  material  under  consideration  are  states  of  thermo¬ 
dynamic  equilibrium;  the  thermodynamic  pressure  is  identified  with  the 
pressure  when  the  same  state  is  reached  quasi -statically. 

Returning  to  the  case  of  confined  compression  (uniaxial  strain) 
discussed  earlier,  consider  a  change  of  state  of  the  material  by,  say,  an 
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adiabatic  process.  Assume  that,  at  stress  conditions  such  that 

(J 22  “  ^11  ^  »  strains  are  small  and  the  material  is  linearly  elastic 

and  rate -independent.  Given  a  reference  state  of  low  hydrostatic  stress 

(a..)  =  -  P  O  ,  the  stress- strain  relations  at  low  stress  magnitudes 

'  ij  o  o  y 

are 


a. . 

IJ 


+  p  ^ 


IJ 


\  € 


kk 


+  2|J.€ 


ij 


and  ~  ^11  ~“^P"Po^^Po‘  conditions  when  ^22  "  =  2a 

by  the  subscript  1.  It  is  readily  verified  that  the  isotropic,  longitudinal 
and  transverse  normal  stress  components  at  this  point  are  given  by 


^1  -  ^0=  ■ 


and  the  density  is  given  approximately  by 

Pi  =  Po  ^  IT-)  • 

The  linear  stress -density  relations  up  to  this  point  are  shown  schematically 
in  Fig.  A-L  The  absolute  maximum  normal  stress  at  state  1  thus  differs 
from  the  mean  normal  pressure  by  the  amount  4?f/ 3.  Now  for  higher 
compressions  the  maximum  stress  difference  is  assumed  to  have  at  most 
the  value  2  a  .  In  materials  such  as  metals  it  has  been  conjectured^^  that 
the  rigidity  modulus  actually  decreases;  this  is  certainly  plausible  in  the 
case  of  easily  crushable  porous  solids.  In  such  a  case  the  absolute  maiximum 
normal  stress  and  the  mean  normal  pressure  approach  each  other  at  high 
compressions  as  sketched  in  Fig.  A-1;  in  addition,  there  must  be  a  cusp  in 
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the  -O' 


curve.  In.  normal  shock  impingement  experiments  on  crushable 


11 

solids  in  a  shock  tube^^,  the  axial  normal  stress  component  in  the  solid  is 
equal  to  the  pressure  measured  in  the  gas  upstream  of  the  specimen,  and  a 
cusp  is  indeed  observed  in  the  curve  of  stress-density  states  reached  by 
shock  in  the  solid. 
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Overdensity  Ratio  (P"Pq)/Pq 


Figo  A 


-1  ASSUMED  BEHAVIOR  OF  SOLID  IN  CONFINED 
COMPRESSION 
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APPENDIX  B 


WASTE  HEAT  CALCU1-.ATION 


The  approximate  theory  of  explosion  growth  described  in  this  report 
requires  a  simple  method  of  estimating  the  so-called  waste  heat  in  shock.  To 
do  this,  it  is  assximed  that  we  are  dealing  with  condensed  matter  whose 
thermal  equation  of  state  is  approximated  by  the  form  P  =  f(V),  where  P 
is  absolute  pressure  and  V  is  specific  volume;  a  third  variable,  such  as 
the  absolute  temperature  T.  is  absent.  The  assumption  is  a  reasonable 
one  for  pure  solids  emd  liquids  in  the  absence  of  phase  changes. 

The  constant-pressure  and  constant-volume  specific  heats  (c^  and 
c^)  of  typical  condensed  matter  are  very  nearly  equal.  The  difference  of 
specific  heats  of  any  pure  substance  is  given  by 


3P  \ 

1 

c  -  c  =  T 

P  V  1 

kiv 

\"^l 

If  c  -  c  =0,  then  either  (8P/aT)„  =  0  or  (av/8T)p  =  0;  both  cases  lead 
to  the  result  that  the  equation  of  state  of  the  material  is  of  the  form 
P  =  f(V).  By  use  of  the  identity 


(3P\ 

ap) 

j 

f  av 

Hv  ^  ■ 

j"3Vj 

T 

rar 

\ 

we  can  write 


=  VTG  /3^  . 


where  G  and  ^  are  respectively  the  isothermal  bulk  modulus  and  iso- 

baric  volume  expansivity  of  the  substance.  For  water  at  standard  conditions, 
for  example,  we  have  G  =  2,  1  x  10^^  dynes/ cm  ,  ^  =  2,  1  x  10  K  , 
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V  =  1.  0  cm^/gm,  T  =  20®C  =  293®K,  and  €^  =  1,0  cal/gm-^Kfl  Hence 
c  -  c  =  0,  0065  Gal/gm-°K,  and  the  equation  of  state  of  water  in  the 
absence  of  phase  changes  is  approximately  of  the  form  P  =  f(V). 

An  equation  of  state  of  the  form  P  =  f(V)  implies  that  all  states 
which  the  material  can  attain  lie  on  a  single  curve  in  the  P-V  plane 
(Fig.  B-1),  In  the  present  approximate  theory,  the  Hugoniot  is  assumed 
to  be  given  by  the  form 


P-P 

o 

C 


(B-1) 


where  C,  M-,  and  n  are  empirical  constants  to  be  determined  by  fitting 
to  experimental  shock  compressibility  data^^.  The  initial  pressure  P^ 
was  neglected  in  comparison  to  the  pressures  of  interest.  If  the  equation 
of  state  is  of  the  form  P  =  f(V),  then  Eq.  (B‘^1),  in  addition  to  being  the 
locus  of  states  which  can  be  attained  by  shock  from  the  initial  state 
V^),  is  also  the  equation  of  any  process  involving  pressure  and  volume 
changes.  After  passage  of  a  shock  pressure  P,  a  particle  of  material 
expands  nearly  adiabatically  to  the  initial  pressure  P^,  the  path  of 
expansion  necessarily  coinciding  with  the  Hugoniot  in  the  P-V  plane. 

It  is  pertinent  to  explore  the  implications  of  these  assumptions. 
From  thermodynamics,  we  have  the  reciprocity  relation  for  any  pure 
substance. 


dE 

W 


'T 


-  T  I  i 

-  T  Iw), 


-  p . 


where  E  is  specific  internal  energy.  Hence  for  this  case 


9E  \ 


-  f(V)  . 
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From  the  definition  of  the  constant-volvime  specific  heat 


3E 

"Ft 


=  c 


but  for  this  case  c^  =  c^  =  c.  Since  E  is  a  state  variable  and  therefore 
its  total  differential  is  exact,  we  must  have 


ac' 

FV 


afl 


=  0; 


T  r^iv 

thus  c  is  at  most  a  function  of  temperature  only. 
We  therefore  write  the  total  differential 
dE  =  c(T)  dT  -  f(V)  dV. 

Assuming  c  =  constant,  integration  gives 


E  -  E^  =  c(T  -  T^)  +  F(V), 


where  we  define 


(B-2) 


(B-3) 


F(V)  =  J 

V 


f(V)  dV  , 


and  the  subscript  zero  denotes  conditions  at  the  state  (Fq»  Evidently 

the  internal  energy  consists  of  additive  thermal  energy  and  deformation 
(or  strain)  energy  fractions. 

By  combination  of  the  conservation  laws  for  mass,  momentum, 
and  energy  flux  through  a  shock,  one  obtains 

E  -  Eq  (P  +  Pj,)(V^  -  V).  (B-4) 

Thus  the  specific  internal  energy  jump  in  shock  is  represented  by  the  trape 
zoid  below  the  chord  of  the  Hugoniot  in  Fig.  B-1.  The  fraction  of  specific 
internal  energy  which  is  recovered  mechanically  as  expansion  work  in  the 
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adiabatic  expansion  to  V^)  is  precisely  F(V):  the  lens -shaped  area 

between  the  Hugoniot  and  its  chord  is  therefore  unrecovered  specific 
internal  energy,  equal  to  c(T  -  T^).  This  "waste  heat"  is  evidenced  as  a 
residual  temperature  rise  in  the  material.  Eventually  the  temperature 
gradients  are  equalized  by  heat  transfer  processes  at  constant  pressure, 
and  none  of  the  shock-induced  thermal  energy  is  mechanically  recoverable. 

If  the  equation  of  state  of  the  material  is  of  the  form  P  =  f(V)  and 
the  Hugoniot  is  given  by  Eq.  (B-1),  then  Eq.  (B-1)  is  evidently  also  the 
thermal  equation  of  state  of  the  material.  Equation  (B-3)  is  the  corres¬ 
ponding  caloric  equation  of  state.  Writing  the  first  law  of  thermodynamics 
for  any  reversible  process,  we  have 

TdS  =  dE  +  PdV  =  c(T)  dT  +  [p-f(V)]  dV, 
where  S  is  specific  entropy  and  Eq.  (B-2)  has  been  used.  Thus 

dS  =  c(T)  ^  , 


and  entropy  is  a  function  of  temperature  only*  Assuming  c  =  constant 
and  integrating, 

S  -  =c  In  (^)  .  (B-5) 

o 

We  can  also  write  the  canonical  equation  of  state  E  =  E(V,S)  by 
substitution  into  Eq.  (B-3): 


E  -  E  =  cT 
o  o 


1  S ,  —  S  \ 

- 

exp 

I  ° 

-  1 

1  ) 

- 

+  F(V). 


(B-6) 


The  state  surface  (P-V-T  surface)  for  the  class  of  condensed  matter 
having  a  thermal  equation  of  state  of  the  form  P  =  f(V)  is  therefore  a  ruled 
surface  with  generators  parallel  to  the  T-axis.  Such  a  surface  is  shown 
in  Fig.  B-2.  The  temperature  scale  can  also  be  regarded  as  a  non-linear 
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entropy  scale  in  view  of  Eq.  (B-5)»  If  the  material  is  raised  from  the 

state  (P  ,  V  ,  T  )  to  the  state  (P  ,  V^,  T  )  by  a  shock  process,  the 
o  o  o  s  s  s 

temperature  difference  (T^  ”  is  given  by  QV^/c  ,  where 

denotes  the  lens^shaped  area  between  the  Hug  oniot  and  its  chord.  The 
adiabatic  (isentropic)  expansion  following  shock  is  also  isothermal  in 
view  of  Eq,  As  the  state  point  moves  to  V^,  T^)  in  the 

adiabatic  expansion,  a  unit  mass  of  the  material  performs  on  its  surround- 
ings  an  amount  of  mechanical  work  equal  to  F{V^).  No  work  is  done  in 
the  constant -pressure  (and  hence  constant « volume)  process  through  which 
the  material  returns  to  its  initial  state. 

It  may  be  concluded  that  the  lens-shaped  area  between  the  Hugoniot 
and  its  chord  furnishes  an  estimate  of  the  energy  made  mechanically  un¬ 
available  by  the  shock  process;  the  estimate  is  accurate  to  the  extent  that 
phase  changes  are  absent  and  the  ratio  of  specific  heats  for  the  condensed 
matter  considered  is  close  to  unity. 
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Pressure 


entropy  scale  in  view  of  Eq,  (B-5),  If  the  material  is  raised  from  the 
state  T^)  to  the  state  (I*g»  Vg»  by  a  shock  process,  the 

temperature  difference  (T^  -  T  )  is  given  by  QV  /c  ,  where  QV 
denotes  the  lens-shaped  area  between  the  Hugoniot  and  its  chord.  The 
adiabatic  (isentropic)  expansion  following  shock  is  also  isothermal  in 
view  of  Eq.  (B-5).  As  the  state  point  moves  to  (P  ,  V  ,  T  )  in  the 
adiabatic  expansion,  a  unit  mass  of  the  material  performs  on  its  surround¬ 
ings  an  amount  of  mechanical  work  equal  to  F(V  ).  No  work  is  done  in 

s 

the  constcint-pressure  (and  hence  constant-volume)  process  through  which 
the  material  returns  to  its  initial  state. 

It  may  be  concluded  that  the  lens-shaped  area  between  the  Hugoniot 
and  its  chord  furnishes  an  estimate  of  the  energy  made  mechanically  un¬ 
available  by  the  shock  process;  the  estimate  is  accurate  to  the  extent  that 

absent  and  the  ratio  of  specific  heats  for  the  condensed 
is  close  to  unity. 


phase  changes  are 
matter  considered 
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Pressure 


Specific  Volume  V 


Fig.  B-1  PRESSURE- VOLUME  RELATIONS  IN  CONDENSED  MATTER 
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Fig.  B-2  P- V-T  SURFACE  AND  SHOCK 

COMPRESSION-EXPANSION  CYCLE 
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